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1 .  INTRODUCTION 


This  Is  the  final  report  on  RADC  Postdoctoral  Contract  No.  CCT-SC-0102-442, 
Task  Order  1.  The  work  statement  covering  the  portion  of  the  project  beginning 
May  8,  1979  and  terminating  on  November  1,  1979  is  presented  below: 

The  objective  of  the  project  Is  to  construct  an  analytical  model  of  a 
class  of  RF  intrusion  sensor  systems  which  can  be  used  to  predict  their 
behavior  under  a  wide  range  of  conditions.  There  are  a  variety  of  RF  intrusion 

A 

sensor  systems  now  In  use  or  under  development?”1  throu9h  8-9  5^  this  work 
will  be  focussed  on  a  specific  configuration?"®’®"4,8”5  ,B"7  In  this  system 
a  signal  Is  generated  at  a  point  along  a  leaky  coaxial  cable**”1  throu9h  C-14 
laid  out  circularly  on  the  periphery  of  an  area  to  be  protected. 

A  receiving  antenna  is  mounted  at  or  near  the  center  of  the  area  to  respond 
to  signal  energy  "leaking"  from  the  small  apertures  placed  along  the  cable. 

The  entrance  of  an  intruder  perturbs  the  field  in  the  vicinity  of  the  cable 
and  this  perturbation  is  sensed  by  the  antenna. 

The  work  will  be  divided  into  four  sequential  tasks,  as  follows: 

Task  1:  A  review  of  the  pertinent  literature  to  determine  what  analytical 
work  has  been  done  that  is  applicable  to  the  problem.  Particular  attention 
will  be  focused  on  papers  dealing  with  the  fields  generated  by  leaky  coaxial 
caMes**"1  throu9h  c*14  and  the  perturbation  of  fields  by  various  kinds  of  ob¬ 
jects,  including  humans  and  animals?”1  t^,rou5h  D-35 

Task  2:  Generation  of  an  analytical  model  for  the  specified  configuration, 
using  past  results  uncovered  in  Task  1  whenever  they  are  applicable  and  gener¬ 
ating  new  analysis  if  (as  is  already  evident)  the  past  work  does  not  cover  the 
entire  problem.  The  model  envisioned  involves  the  use  of  theory  already  formu¬ 
lated  on  the  fields  around  a  leaky  coaxial  cable  in  free  space,  generalization 
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of  that  theory  to  include  the  effect  of  the  ground  on  these  fields,  and 
finally  the  perturbation  of  the  fields  due  to  foreign  objects  entering  the 
environment.  Particular  emphasis  will  be  placed  on  the  polarization  character¬ 
istics  of  the  unperturbed  and  perturbed  fields,  with  a  view  toward  the  use  of 
polarization  changes  as  a  means  of  identifying  and  tracking  intruders. 

Task  3:  Construction  of  a  Fortran  program  based  on  the  analytical  results 
obtained  in  Task  2.  This  program  is  expected  to  contain  a  large  number  of  un¬ 
specified  variables.  Among  those  should  be  the  radio  frequency,  constitutive 
parameters  of  the  ground  and  those  of  foreign  objects  perturbing  the  fields, 
size  and  shape  parameters  of  these  objects,  polarization  characteristics  of 
the  fields  and  other  variable  parameters  associated  with  the  cable,  receiving 
antenna  and  the  intruding  object. 

Task  4:  A  parametric  study  based  on  the  computer  program  constructed  in 
Task  3  to  determine  the  effects  of  various  parameters  on  the  expected  behavior 
of  the  system. 

***** 

The  contract  was  extended  beyond  November  1,  1979.  However,  there  was  a 
hiatus  between  that  date  and  March  1980,  when  the  second  portion  of  the 
contract  officially  began.  The  extension  covered  the  period  from  the  above 
date  to  November  1,  1980. 

The  work  statement  for  the  second  portion  of  the  contract  is  presented 
below: 

Tne  objective  of  tne  project  is  to  analyze  a  particular  class  of  RF 
intrusion  systems.  Results  will  be  used  to  predict  their  behavior  under  a 
wide  range  of  conditions.  There  are  a  variety  of  RF  intrusion  systems  now  in 
use  or  under  development,  but  this  work  focuses  on  a  specific  configuration. 

In  this  system  a  signal  is  generated  at  a  point  along  a  leaky  coaxial  cable 


laid  out  circularly  on  the  periphery  of  an  area  to  be  protected. 

A  receiving  antenna  is  mounted  at  or  near  the  center  of  the  area  to 
respond  to  signal  energy  "leaking"  from  the  small  apertures  placed  along  the 
cable.  The  entrance  of  an  intruder  perturbs  the  field  in  the  vicinity  of  the 
cable  and  this  perturbation  is  sensed  by  the  antenna. 

Work  on  this  problem  was  initiated  in  May  1979.  A  comprehensive  review 
of  the  pertinent  literature  was  completed.  An  analytical  model  was  constructed 
which  includes  (a)  mathematical  expressions  for  the  fields  generated  by  the 
cable  slots  in  free  space,  (b)  the  effect  of  ground  reflections  on  the  fields 
from  the  slots,  (c)  a  generalized  mathematical  model  for  the  scattered  fields 
not  accounting  for  ground  reflections  of  these  fields,  and  (d)  a  computer  pro¬ 
gram  for  the  slot-induced  fields  both  with  and  without  ground  reflections. 

The  subject  task  constitutes  an  enlargement  of  the  scope  of  the  work  to 
include  detailed  results  on  the  fields  scattered  from  human  frame  targets. 

The  projected  phases  of  this  task  are: 

Phase  1:  An  analysis  of  the  scattering  from  human  frame  targets,  where 
the  target  is  modelled  as  a  lossy  dielectric  of  ellipsoidal 
or  spheroidal  shape. 

Phase  2:  Incorporation  of  the  scattering  results  of  Phase  1  into  the 

generalized  mathematical  model.  The  scattered  fields  obtained 
will  include  the  effect  of  the  fields  from  the  slots  impinging 

i 

directly  on  the  target  and  the  ground- reflected  slot- induced 
fields  also  incident  on  the  target. 

Phase  3:  Analysis  of  the  effect  of  ground- reflections  on  the  fields 
obtained  in  Phase  2. 

Phase  4:  Development  of  a  computer  program  incorporating  the  results  of 
the  analysis  done  in  Phases  1,  2  and  3  and  containing  the 


previously  developed  program  [see  (d)  above]  as  a  subprogram. 
Phase  5:  Use  of  the  computer  program  developed  in  Phase  4  to  obtain  a 
set  of  numerical  results  showing  the  variation  of  the  signals 


received  at  the  antenna  with  key  system  parameters. 

***** 

The  report  describes  the  very  extensive  mathematical  analysis  undertaken 
to  meet  the  above  objectives.  Numerical  results  are  presented  for  a  limited 
range  of  parameter  values.  The  limitations  are  those  of  time  and  cost.  The 
analytical  model  and  the  resulting  computer  programs  are  sufficiently  general 
to  include  the  capability  of  treating  a  much  wider  range  of  geometries  and 
parameter  variations.  The  outputs  of  the  project  should  be  considered  as  two¬ 
fold.  First,  the  particular  geometries  for  which  computer  results  are  shown 
are  of  direct  importance  in  themselves.  Secondly,  the  computer  program  should 
be  considered  as  an  output  of  the  work,  since  it  could  now  be  used  to  study  a 
wide  range  of  cases  that  may  be  of  practical  interest. 

The  body  of  the  report  is  divided  into  11  sections.  There  are  also  five 
appendices.  Section  2  contains  a  mathematical  description  of  the  basic  model 
of  the  system.  The  generic  configuration  being  analyzed  consists  of:  (1)  a 
set  of  apertures,  small  compared  with  wavelength,  on  which  electric  and  magnetic 
fields  can  be  specified.  These  apertures  are  placed  somewhere  in  the  vicinity 
o*  a  flat  ground  surface;  (2)  an  electromagnetic  scatterer  modelled  to  resemble 
an  "intruder";  (3)  an  antenna,  placed  at  an  arbitrary  point,  which  receives 
signals  resulting  from  the  field  distributions  on  the  apertures.  The  received 
signals  are  influenced  by  the  presence  of  the  ground  and  should  also  be  in¬ 
fluenced  by  the  presence  of  the  intruder.  It  is  clear  that  the  major  objective 
of  the  work  is  to  determine  the  magnitude  and  polarization  of  the  fields  at  the 
receiving  antenna  both  with  and  without  the  presence  of  the  intruder  and  to 
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determine  the  effect  of  the  intruder's  presence  on  the  magnitude  and  polariza¬ 
tion  of  the  fields. 

In  order  to  obtain  numerical  results,  it  is  of  course  necessary  to  specify 
more  tightly  the  set  of  apertures,  the  scatterer  used  to  reDresent  the  intruder 
and  the  location  of  the  antenna.  Hence,  in  the  situation  actually  modelled  on 
the  computer,  the  apertures  are  small  slots  along  a  coaxial  cable  laid  along 
the  ground  surface  in  a  circular  pattern.  The  antenna  is  somewhere  near,  but 
not  necessarily  exactly  at,  the  center  of  this  circular  configuration.  The 
scatterer  is  a  spheroid  somewhere  near  the  cable.  The  spheroid  must  have  uni¬ 
form  constitutive  parameters  but  aside  from  that  constraint  is  not  restricted, 
since  it  is  in  general  a  lossy  dielectric.  The  frequency  range  of  validity  of 
the  analysis  is  roughly  from  50MHz  to  500MHz,  implying  a  wavelength  range 
between  roughly  0.6  and  6  meters.  The  cable  diameter  is  less  than  two  centi¬ 
meters  and  the  length  of  each  slot  on  the  cable  is  less  than  one  centimeter, 
hence  slot  dimensions  are  always  small  compared  with  wavelength.  This  is  a 
feature  that  simplifies  the  analysis  somewhat. 

It  was  not  possible  within  the  time  and  cost  constraints  of  the  project 
to  model  the  cable  fields,  the  scatterer  and  the  effect  of  the  ground  in  a 
rigorous  manner.  "Engineering  approximations"  were  necessary  to  render  the 
analysis  feasible.  Before  delineating  these  approximations,  we  will  first 
describe  the  basic  methodology. 

The  electric  and  magnetic  fields  were  all  expressed  as  superpositions 
of  plane  waves.  This  "plane  wave  spectrum"  type  of  analysis  can  be  viewed  as 
a  process  of  three-dimensional  Fourier  transformation  from  the  space  of  a  set 
of  three  position  coordinates  (x,  y,  z)  to  the  space  of  a  set  of  three  com¬ 
ponents  of  a  wave  propagation  vector  k0  B  =  (kQ  Bx,  kQ  B  ,  kQ  6Z).  Since  kQ, 
the  propagation  constant  of  free  space,  is  a  constant,  we  deal  with  the  vector 
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space  (Bx,  By,  3Z),  subject  to  the  constraint  Bx  +  By  +  Bz  =  1 .  Because 
of  that  constraint,  only  two  of  the  components  By,  Bu  or  B,  (say  By,  Bu)  can 
be  specified  independently  and  the  other  (say  Bz)  can  be  calculated  therefrom. 
Hence,  the  process  degenerates  to  a  two-dimensional  Fourier  transformation  be¬ 
tween  two  position  coordinates  (say  x,  y)  and  two  components  of  B  (say  B„,  B  ). 

-*•  x  y 

Section  II  of  the  report  describes  the  general  mathematical  modelling  of 
the  problem.  The  plane-wave  spectrum  of  the  field  at  the  antenna  consists  of 
four  contributions,  as  follows: 

(1)  the  direct  wave  from  the  slots  in  the  cable  [as  if  the  ground  and 
intruder  (scatterer)  were  not  present] 

(2)  the  ground-reflected  wave  from  the  slots  in  the  absence  of  an  in¬ 
truder  (scatterer) 

(3)  the  field  scattered  by  the  intruder  directly  into  the  antenna 

(4)  the  ground-reflected  field  at  the  antenna  resulting  from  scattering 
by  the  intruder. 

The  plane-wave  spectrum  type  of  modelling  enables  us  to  determine  the 
contributions  (2)  and  (4)  using  standard  theory  of  reflection  of  a  plane-wave 
from  an  infinite  surface.  Another  very  important  reason  for  using  that  type 
of  modelling  is  the  fact  that  the  scattering  program  we  are  using  (the  best 
available,  in  the  writer's  opinion,  for  the  purpose  at  hand)  assumes  a  plane- 
wave  input.  Thus,  by  considering  the  fields  at  every  stage  of  the  process 
as  the  fields  of  a  plane-wave,  we  are  enhancing  the  accuracy  of  the  calcula¬ 
tions  (as  compared  with  approximating  fields  of  spherical  waves  or  near-zone 
waves  or  other  more  complicated  field  patterns  as  plane-waves). 

However,  the  negative  feature  of  this  kind  of  modelling  is  that  it 


necessitates  a  process  of  inverse  two-dimensional  Fourier  transformation  at 
the  last  stage  of  the  calculation  in  order  to  convert  the  field  spectra, 


functions  of  (6y,  8U),  into  functions  of  the  horizontal  position  coordinates 
a  y 

(x,  y).  That  process  is  computer-time  intensive  and,  hence,  the  computations 
require  more  computer  time  than  would  simpler  but  more  approximate  modelling 
schemes  we  could  have  used. 

Section  3  and  Appendix  I  contain  material  on  the  method  variously 
called  the  “Kirchhoff-Huyghens"  or  "Stratton-Chu"  integral  method  of  deter¬ 
mination  of  the  fields  at  an  arbitrary  point  in  space  due  to  fields  on  a 
closed  surface  surrounding  that  point.  This  technique,  widely  used  in  antenna 
theory,  is  our  method  of  calculating  the  fields  from  the  slots  along  the  cable. 
The  method  is  presented  in  general  in  Appendix  I.  In  Appendix  II  the  general 
concept  of  the  plane-wave  spectral  representation  of  fields  is  presented. 
Appendix  III  covers  the  spectral  representation  of  the  fields  from  an  aperture 
as  calculated  using  the  Kirchhoff-Huyghens  integral  method.  Section  3  of  the 
main  body  of  the  report  specializes  the  analysis  in  Appendix  III  to  the  case 
of  the  rectangular  slot.  The  use  of  a  rectangular  slot  shape  for  mathematical 
convenience  is  an  approximation  to  the  actual  shape  of  slots  on  the  cable  of 
central  interest.  However,  since  slot  dimensions  are  small  compared  with 
wavelength,  the  results  are  not  significantly  dependent  on  slot  shape  but  only 
on  slot  area  regardless  of  shape. 

Section  4  specializes  the  analysis  in  Section  3  to  the  specific  case 
of  a  slotted  coaxial  cable  laid  in  a  circular  pattern,  which  is  the  actual 
case  treated  in  our  numerical  work.  The  wave  modes  propagating  in  the  cable 
are  detailed  in  Appendix  IV.  The  analysis  in  Section  4  is  the  basis  of  Sub¬ 
routine  SOURCE,  which  is  our  computer  program  to  evaluate  the  fields  from  the 
cable  as  if  the  cable  were  in  free  space. 

Section  5  treats  the  reflection  of  plane-wave  fields  from  the  ground 
surface.  This  analysis  is  the  basis  for  a  section  in  the  main  program  which 


is  applied  twice,  first  to  the  field  from  the  cable,  and  secondly  to  the  field 
from  the  scatterer. 

Section  6  treats  the  coordinate  transformations  required  to  enter  the 
scattering  program.  The  latter  is  called  “Subroutine  BARBER",  after  Pro¬ 
fessor  Peter  Barber  of  the  Bioengineering  Department  at  the  University  of 
Utah,  who  developed  the  program  and  kindly  gave  us  the  program  cards  for  use 
in  Northeastern's  VAX  computer.  Barber's  program,  discussed  in  Section  7, 
considers  an  incoming  plane-wave  specified  by  a  magnitude  and  polarization  in 
a  "lab  frame"  (Barber's  terminology)  coordinate  system.  The  basic  coordinate 
system  used  in  our  problem,  which  we  call  the  "ground  frame",  is  the  co¬ 
ordinate  system  in  which  our  "Subroutine  SOURCE"  delivers  the  field  components 
from  the  cable  back  into  the  main  program.  In  the  main  program,  the  ground 
reflection  operation  is  performed,  and  the  superposition  of  plane-wave  spectra 
of  direct  and  ground- reflected  source  fields  is  then  evaluated.  The  resultant 
field  components  are  expressed  in  the  ground  frame.  The  analysis  discussed 
in  Section  6  describes  the  transformation  of  field  components  between  Barber's 
lab  frame  and  our  ground  frame,  so  that  the  input  to  Barber's  scattering  pro¬ 
gram  can  be  expressed  in  his  lab  frame.  The  program  to  implement  the  analysis 
in  Section  6  is  called  "Subroutine  BIS"  ("Barber  Input  Subroutine") . 

As  indicated  above,  the  Barber  scattering  program  is  discussed  in  Section 
7.  In  Section  8,  we  present  the  coordinate  transformation  required  to 
transform  the  output  of  the  scattering  program,  i.e.,  the  scattered  field 
components,  from  Barber's  lab  frame,  in  which  these  components  are  expressed 
in  his  program,  back  to  the  ground  frame.  The  program  based  on  this  analysis 
is  called  "Subroutine  BOS"  ("Barber  Output  Subroutine"). 

In  Section  9  we  describe  the  analysis  of  the  effect  of  ground  reflection 
on  the  scattered  field.  Computationally,  the  implementation  of  this  analysis 


is  accomplished  in  the  main  program,  which  receives  from  Subroutine  BARBER 
via  Subroutine  BOS  the  scattered  waves  in  two  directions,  one  being  that 
directed  from  the  center  of  the  scatterer  to  the  antenna  position,  the  other 
being  that  directed  from  the  center  of  the  scatterer  to  a  "ground  reflection 
point",  which,  as  dictated  by  the  law  of  reflection,  sends  the  ground- reflected 
wave  toward  the  antenna.  The  first  of  these  scattered  fields  is  logged  in  the 
main  program  as  that  received  by  the  antenna  [Contribution  (3)]  and  the  second 
is  driven  through  the  ground-reflection  process  ir.  the  main  program  and  the 
result  becomes  Contribution  (4)  at  the  antenna. 

The  justification  for  considering  scattering  in  only  these  two  directions 
as  opposed  to  all  directions  (which  would  have  required  prohibitively  large 
computer  time  and  would  have  required  still  another  double  Fourier  trans¬ 
formation),  is  the  fact  that  the  scatterer  is  in  the  far-zone  of  the  antenna 
and  subtends  a  very  small  angle  at  the  antenna,  so  that  the  antenna  sees  it 
as  very  nearly  a  point  source.  Hence,  it  appears  as  a  plane-wave  at  the 

antenna.  Barber's  scattering  program  produces  a  plane-wave  in  any  given  di- 

jk  r 

rection,  which  when  multiplied  by  the  Green's  function  e  /r,  can  be  viewed 
by  an  observer  at  the  antenna  as  a  spherical  wave  from  that  direction,  the 
curvature  of  whose  phase  front  Is  negligible,  i.e.,  equivalent  to  a  plane-wave 
weighted  by  -p  .  The  ground-reflected  scattered  wave,  to  be  treated  rigorously, 
would  require  that  scattered  plane-waves  in  all  directions  be  reflected  from 
the  ground  and  that  the  wave  seen  at  the  antenna  is  the  superposition  of  all 
of  these  waves.  However,  the  theory  of  the  plane-wave  spectral  representation 
of  fields^”1 ,A'2  shows  that,  if  the  receiving  point  is  sufficiently  far  away,  the 
only  significant  contribution  is  that  which  obeys  the  law  of  reflection.  Hence, 
according  to  this  approximation.  In  view  of  the  large  distance  between  the 
antenna  and  the  scatterer,  and  the  fact  that  they  are  each  in  the  other's  "far 
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zone",  it  is  sufficient  to  consider  only  one  ground- reflected  plane-wave 
from  a  single  ground-reflection  point.  Invoking  this  approximation  saves  an 
enormous  amount  of  computer  time  since  it  saves  us  from  another  round  of  double 
Fourier  transformations. 

Section  10  contains  the  analytical  basis  for  the  way  in  which  the  overall 
computation  must  be  done  to  provide  for  the  inverse  Fourier  transformation  of 
the  composite  fields.  This  is  done  in  the  main  program,  with  the  aid  of  two 
small  subroutines  labeled  C0MPUTEC0  and  COMPUTEEF.  The  first  of  these  sub¬ 
routines  performs  a  numerical  integration  of  each  component  of  the  spectral 
field  vector  on  the  angle  [=  tan"1  (Sy/gx)]  from  0  to  2ir.  The  second  performs 
another  numerical  integration  on  (=  Jb*  +  )  from  0  to  1 .  This  double 

integration  procedure  is  done  in  computing  the  inverse  Fourier  transform  of  the 
field  components  at  the  antenna  in  the  absence  of  the  scatterer. 

It  was  attempted  to  carry  out  this  same  procedure  in  the  case  where  the 
scatterer  is  present.  In  this  case,  the  inverse  Fourier  transform  must  be 
computed  at  the  position  of  the  center  of  the  scatterer,  rather  than  the  an¬ 
tenna  position.  It  was  found  that,  with  the  scatterer  present,  computer  time 
becomes  prohibitively  large  when  the  above  procedure  is  used.  It  was  decided 
to  make  use  of  the  stationary  phase  principle  to  evaluate  (approximately)  the 
integral  on  <j>&  analytically.  A  justification  for  the  use  of  this  method  is 
presented  in  Section  10. 

Section  11  contains  numerical  results  and  conclusions  therefrom. 
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2.  MATHEMATICAL  MODELLING  OF  RF  INTRUSION  SENSOR  SYSTEM 


The  generic  system  that  is  the  subject  of  this  study  is  shown  in  Figure  2.1. 
The  point  0  on  the  diagram  is  the  origin  of  coordinates. 


Scatterer 


Consider  a  collection  of  N  sources  numbered  with  an  index  1  ranging  from 
1  through  N.  The  source  numbered  l  is  centered  at  a  position  designated 
by  a  vector  r^  originating  at  0  and  terminating  on  P^.  It  occupies  a  small 
volume  around  the  point  P^,  but  this  volume  is  assumed  to  subtend  a  very 
small  solid  angle  when  viewed  from  the  position  of  a  receiving  point  and  hence 
appears  to  be  a  point  source  at  the  receiving  point.  There  may  be  a  number 
of  receiving  points;  hence  they  are  numbered  with  an  index  m,  ranging  from  1 
through  M.  The  mth  receiving  point  Pm  Is  designated  by  a  vector  r,  originating 
at  0  and  terminating  on  P  . 

The  scatterer  is  centered  at  a  point  P0$,  designated  by  a  vector  rQS 
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originating  at  0  and  terminating  on  PQ$. 

The  point  P  is  an  arbitrary  point  In  space,  designated  by  a  vector  r 
emanating  from  0  and  terminating  on  P. 

We  denote  the  electric  field  vector  at  the  point  P  due  to  excitation  of  the 
t**1  source  by  a  three-element  column  vector  C^^)]  whose  elements  are  the 
rectangular  components  of  E  In  the  basic  (x,y,z)  coordinate  system.  This 
system,  as  shown  in  Figure  2.';,  has  its  z-axls  in  the  vertical  direction  (up¬ 
ward)  and  its  (x,y)  plane  along  the  ground  surface.  Thus,  the  electric  field 
at  the  arbitrary  point  P  due  to  excitation  at  source  l  is: 


"tx 


(r) 


£E,  (r)] 
-► 


Ety  ^ 
Etz  ^ 


(2.1) 


Using  the  plane-wave  spectral  representation  for  the  electric  field 
(Appendix  II,  Eq.  II-3),  we  have 


A-l .A-2 


"  I! 


2  .  _Jk!h'£  ,.-Jk|S2|Zrf 


d  & e 


J*k|B2!  z 


(2.2) 


where  r  *  o  zz;  p  «  xx  +  yv 


f7*1*  ■ 


1  “h  ■  2®x  *  ley  • 


-  and  +  refer  to  the  cases  S2  ■  - J 0Z J  (downward  propagation)  and  S2  ■  +  |b2( 
(upward  propagation),  respectively. 
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where 


ll+.x^ 

St+,y^ 


|t+,z^ 


The  magnetic  field  vector  at  Pm  due  to  excitation  at  source  l 


tyy]  ■ 


r  O  Jk6h'Pm  -Jk } 6, | Zm 
1  d2  8h  e  *”  (e  "t£.(&)]CEe.(|h)] 


j k I 0_ I  2 

*=  "Wyityiy]) 

where  (from  the  Maxwell  equations  H  *  ^ — 7  x  E) 

>  JUVq  -*■ 

o  ±!szl  sy 

?|S2I  0  -Bx 

-6y  6x  °_ 

The  fields  E.  and  H9  have  four  contributions  as  follows: 

t|,(r)]  .  [E<a><r)]  *  CE<b>(r)]  ♦  [|<c»(rl]  *  [E<d>(r>] 
tH.(r)]  •  [H(a)(r)]  *  CH<b,(r)]  *  tH<c>(r)3  *  [H<d,(r» 

Xj  “►  %  ^  **  ~  > 


given  by 

(2.3) 

(2.4) 

(2.5a) 

(2.5b) 


where  the  superscripts  a,  b,  c,  d  correspond  to  the  following: 


(a)  Direct  wave  from  source  i  to  point  P  in  infinite  free  space  (i.e.. 

In  absence  of  ground  or  scat terer) 

(b)  Reflected  wave  from  ground  due  to  excitation  of  source  l  in  absence 
of  scatterer,  as  observed  at  point  P 

(c)  Direct  wave  from  scatterer  (as  if  the  scatterer  were  in  infinite 
free  space)  where  the  fields  incident  on  the  scatterer  consist  of 
the  superposition  of  the  direct  wave  from  P^  to  the  scatterer  and 
the  ground-reflected  wave  at  the  scatterer  resulting  from  excitation 
of  source  i 

(d)  Ground-reflected  wave  at  P  in  response  to  the  scattered  wave  fields. 

It  will  be  shown  in  Section  5  (see  5.18a)  that  the  plane  wave  spectrum 
of  the  reflected  wave  field  (superscript  (b))  is  related  to  that  of  the  in¬ 
cident  wave  field  (superscript  (a))  by  the  expression 

■  CRE(8h)][|1.(6h)]  (2.6) 

where  [R^( 0^) D  is  given  by  Eq.  (5.18a)  and  Eq.  (5.19).  The  receiving  point 
is  always  above  ground-level;  hence  the  subscripts  and  the  incident  wave  at 
the  ground  reflecting  point  always  propagate  downward;  +  and  -  are  used  on  the 
reflected  and  Incident  waves,  respectively. 

The  plane-wave  spectrum  of  the  reflected  wave  as  given  by  Eq.  (2.6)  cor¬ 
responds  to  the  wave  with  superscript  (b)  while  the  spectrum  of  the  incident 
wave  corresponds  to  that  with  superscript  (a).  Using  these  superscripts  in 
their  appropriate  places  in  Eq.  (2.2),  with  the  aid  of  Eqs.  (2.5a,b)  and  (2.6), 
we  obtain  the  following  results: 
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Case  A:  Presence  of  ground  not  accounted  for;  scatterer  not  present.  Fields 
at  point  P  are  [see  Eqs.  (2.2)  and  (2.5a)] 


a 


ft  o  Jkeh*P  -jk | B_ | 2  /.x 

JJ  d  |h  e  e  <2-7a> 


if  source  is  above  P 


,  Jk&'P  jk|62|z  („ 


“  lh» 


where 


if  source  is  below  P. 


(fy,)  a  spectrum  of  electric  field  due  to  source  «.  in  infinite 
free  space. 


and  where 


Bz  *  +|BZ|  for  subscript  + 


|  Bz 1  for  subscript  - 


Also,  from  Eqs.  (2.3),  (2.4)  and  (2.5b) 


CS»(I):lca«  A  '  d  ®h  e 


o  jkBu-P  - jk | S_ | z  ,  x 

2  a  -*■"  -*•  o  2  ro  fa  nrc^aJ 


CP-(^h)]t4-  (!h)]  (2'7b) 


if  source  is  above  P 


l\ “2  a 


2  ■  _JkV£  .jk|S2|z 


if  source  is  below  P. 


Case  B:  Presence  of  ground  1^  accounted  for;  scatterer  not  present.  The 
fields  are  (see  Eqs.  (2.2,  (2.3),  (2.4),  (2.5a,b)  and  (2.6)). 


:[a)(r)]  +  Cf[b)(r)] 

(2.8a) 

|ia)(r)]  +  [H[b)(r)3 

(2.8b) 

where  [E^]  and  [H^ab  are  given  in  general  by  Eqs.  (2.7a)  and  (2.7b), 
respectively,  with  the  aid  of  Eq.  (2.4),  and  where  (from  Eqs.  (2.2),  (2.3), 
(2.4),  (2.5a,b)  and  (2.6)), 


<:>}««  b 


oo 


jksh*p  jk I e  1 2 

e  +  e  z 


(2.9a) 


[s!ib)<r>  W  b 


d2  & e 


e  2  [Meh)]cnE(8hn[E<!)(!hn 

(2.9b) 


where  [P+(|h)]  Is  given  by  Eq.  (2.4)  and  the  matrix  [1^(6^)]  is  given  in 
Section  5  [Eqs.  (5.18a)  and  (5.19)]. 


Case  C:  Presence  of  ground  js_  accounted  for  in  determination  of  the  fields 
incident  on  the  scatterer  but  not  in  evaluating  the  scattered  field 
(i.e.,  ground-reflection  of  scattered  wave  is  neglected).  In  this 
case,  the  fields  at  point  P  are:  (see  Eqs.  (2.4  through  2.8b)) 

C  •  ‘  CEib)(r)]  ♦  [E<C>(X)]  (2.10a) 
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i 


(2.10b) 


[HA(jr)]caSe  c  *  ^a)y  ]  +  CH^b)(r)J  +  [H<c)(r)] 

where  CE[a^(r)],  [Hja^(r)],  [E^(r)]  and  [H^(r)]  are  given  by  Eqs.  (2.7a), 
(2.7b),  (2.8a),  (2.9a)  and  (2.8b),  (2.9b),  respectively,  aided  by  Eqs.  (2.2) 
through  (2.6).  The  field  vectors  [E^(r)]  and  [H^(r)l  will  be  discussed 
in  what  follows. 

Referring  to  Figure  2.1,  we  can  express  the  scattered  field  in  the  form 


[|<c)(r)]  ■  \\  d2  6h  *  t|Lc>(|h) 


-jk|8zlz 


Jkl6  jz 

+  e'c'(8  )  e  2  } 
+  h+  lSh;  e  ' 


(2.11) 


where 

„)  ■  cs‘-)(6h)Hii!Wb,(8l,)j 

I'fdft)  *  tsit)(eh)][li!)+(b,(s„)] 

*  Csi*,(8|l!Kli:H<b)(8h)] 


and  where 


and 


are  3x3  matrices,  and 


C|(aK(b)(Sh)]  .  ♦  [E<b>(5h)] 


We  will  not  discuss  the  details  of  this  constituent  of  the  field  at  this 


point. 


3.  SPECTRAL  FIELDS  FROM  A  RECTANGULAR  SLOT 


We  will  approximate  the  slots  in  the  cable  as  rectangular.  Since  the 
slots  are  always  small  relative  to  wavelength,  the  exact  shape  of  the  slot 
will  not  be  critical.  This  supposition  is  borne  out  by  the  approximations  that 
can  be  justified  when  actual  parameter  values  are  assigned,  as  will  become 
evident  later  in  this  report. 

Let  us  assign  a  set  of  coordinates  (x^,  y^,  z^)  to  the  £th  slot,  together 
with  a  "length"  L^  in  the  x^  direction  and  a  "width"  in  the  y^  direction  and 
a  slot  center  located  at  a  point  (x^Q,  y^,  z^o).  We  note  that  the  ix^  *  yp 
plane  is  in  the  plane  of  the  slot  and  hence  =  0.  We  have  not  specified  the 
origin  of  this  coordinate  system. 

Referring  to  Appendix  III,  Eqs.  (III.6a,b)  and  (III.7a,b),  points  on  the 
slot  are  designated  by  a  vector 


r' 


r1 

-U 


(3.1) 


where  r^Q  is  the  vector  representing  the  center  of  the  slots. 

The  vector  referred  to  in  Eqs.  (III.6a,b)  and  (III.7a,b)  emanates  at  the 
origin  of  the  basic  (x,  y,  z)  coordinate  system.  Therefore,  it  would  be  cor¬ 
rect  to  designate  the  vector  r'  as  follows: 


r'  = 
Uo 


XX 


to 


Iho  +  zz 


to 


(3.2) 


where  (xi0,  ylQ,  z^)  are  the  slot  center  coordinates  in  the  basic  (x,  y,  z) 
system. 

The  vector  ar',  on  the  other  hand,  which  refers  to  the  displacement  of  a 
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point  on  he  slot  surface  from  the  slot  center,  should  be  expressed  In  the 
slot  coordinates  as 


2i 


+ 


(3.3) 


/s 

where  (x^,  y^,  zp  are  the  unit  base  vectors  in  the  slot  coordinate  system  and 
Ax^  and  Ay^  are  the  displacements  from  the  slot  center  in  the  and  y^  di¬ 
rections,  respectively.  There  is  no  displacement  in  the  z^  direction,  so  it 
is  clear  that  Az^  s  0. 

We  will  now  make  an  additional  assumption  which  will  apply  in  every  case 
considered,  namely  that  the  most  general  spatial  dependence  of  key  terms  in 
the  vector  field  components 


n'xE,  n'  •  E,  n1  x  H  and  n*  •  H 


in  Eqs.  (III. 6a, b)  and  (III. 7a, b)  is  exponential  in  and  y^.  Moreover,  the 
exponential  dependence  is  the  same  for  all  four  of  these  vectors.  Mathematically, 
we  can  express  this  with  the  statements: 


*  E  ■  i  6 


, .  yko(«nxi+Vi> 


(3.4a) 


- 1  *  i  cbM(zi> e 
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(3.4b) 
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(3.4c) 
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(3.4d) 


H  =  /  C 
‘  n 
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,(zj) 


.jko,ci.x«t'ny;) 


( £  \ 

whore  it  is  indicated  that  the  vectors  and  Cfna  and  tin?  scalars  Cb|)?  and 

^dnJi  are  ™  9enera^  functions  of  z^  and  Cn  and  r>f)  are  in  general  complex. 

We  should  carefully  note  here  the  fact  that  the  coordinates  (x^,  z‘  ) 

are  the  coordinotcs  of  the  slot  center  in  the  (x^,  y£,  z,! )  system.  The  origin 

of  this  system  will  generally  be  placed  at  the  launching  point  for  the  energy 

that  drives  the  source  and  this  origin  has  no  relationship  to  the  origin  of 

the  basic  (x,  y,  z)  coordinate  system.  Failure  to  take  note  of  this  point 

could  lead  to  confusion  in  the  formulation  of  the  equations  to  follow. 

Substituting  Eqs.  (3.4a,b,c,d)  into  Eqs.  (III. 6a, b)  and  (III. 7a, b),  we 

would  obtain  [where  6-  =  Ch  +  z|fl_|  and  (3;;,,  3V.  are  the  (x‘  y*  z') 

-y  >n  -  >  z  Xj  yj  z £  y.  x.  y. 

components  of  f’,+ .  Note  that  the  components  of  3^  in  the  basic  x,  y,  z  system 


are  Bv,  3W,  +  IflJ  and  the  components  of  the  slot  center  in  that  system  are 

xto’  yto’  Z£o’^A1so  note  that  we  have  not  indicated  explicitly  the  spatial 
dependence  of  the  modified  spectral  fields  in  what  follows.] 

tMSh>],  •  “ft—  l  ^k°<6,,VV£0il82'2t0Vl|(a),)],)n  (3.5a) 
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_  (j  •o.Ka.VW 
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(3.7b) 


The  form:.  (3.7a,b)  will  be  used  in  Section  4  as  the  basis  for  development 
of  the  field  expressions  to  be  used  in  our  final  results. 
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4.  SPECTRAL  FIELDS  FROM  SLOTS  IN  A 
COAXIAL  CABLE  LAID  CIRCULARLY 


We  consider  the  source  configuration  to  consist  of  a  leaky  coaxial  cable 
laid  In  a  circular  pattern,  with  the  radius  of  curvature  or  the  configuration 
sufficiently  large  to  justify  considering  the  propagation  down  the  cable  to 
be  equivalent  to  that  along  a  straight  cable.  The  geometry  is  shown  in  Figure  4.1 
The  local  cylindrical  coordinates  used  within  the  cable  in  the  vicinity 
of  the  slot  are  designated  as  (rp  #p  zp  with  corresponding  rectangular 
coordinates  (x£,  y£,  zp.  The  inner  and  outer  conductor  radii  are  denoted  by 
a  and  b  respectively.  The  finite  thickness  of  the  outer  conductor  is  neglected 
in  designating  the  radial  coordinate  of  the  slot;  hence,  r£  =  b  in  all  cases. 

The  radius  of  the  cable  configuration  is  denoted  by  r.  In  all  cases  to 
be  considered,  it  will  be  true  that 


k  a  «  1 
o 


k.  b  «  1 
o 


R  »  a 


(4.1a) 

(4.1b) 

(4.1c) 


R  »  b 


(4. Id) 


kQ  i?  »  1 


(4.1e) 


To  obtain  needed  relationships  between  coordinate  systems,  we  consider 


Figure  4.2a,b,c,d. 


Figure  4-1.  Geometry  of  cable 


Slot  #t 


Coordinates  of  point  P 
In  system  =  (x$,  y$,  z$) 

x2  *  z3  cos  ♦*,  -  ^3  sin  ♦ic 
y2  *  x3 

z2  *  z3  s1n  ♦to  +  y3  c°s  ♦Jo 
x3  *  y2 

y3  *  ”x2  sin  ♦to  +  22  cos  ♦t, 
z3  a  x2  cos  ♦to  +  z2  sin  ♦to 


(e) 


Figure  4-2.  Coordinate  Systems  (cont'd) 


From  the  equations  In  Figure  4.2a,b,c,d  the  point  P  has  (x,  y,  2) 
coordinates  given  by 

x  *  *  cos  *l0  +  x\  -  R  cos  <pl0  +  x£  cos  <pZ0  -  y'2  sin  ^ 

=  *  cos  ♦to  +  <z3  cos  ♦Ho  "  y3  sin  ♦Ho)  cos  ho  •  x3  sin  ho 

3  *  cos  ho  +  C(b  +  ZV  cos  ♦Ho  '  y4  sin  ho]  cos  ♦to  *  x4  sin  ♦to 

(4.6a) 

y  3  R  s1n  *£o  +  y\  ■  *  sin  *£o  +  x£  sin  *£o  +  y£  cos  4£q 
3  *  sin  +  (z$  cos  qQ  -  y‘  sin  ^Q)  sin  *£q  +  x^  cos  *£o 

3  R  sin  ♦to  +  [(b  +  z4}  cos  ho  *  y4  sin  ♦to3  sin  ♦to  +  x4  cos  ♦to 

(4.6b) 

z  =  b  +  z-j  »  b  +  z'2  =  b  +  (z^  sin  <t>^Q  +  yj  cos  <f>"Q) 

=  b  +  [(b  +  z^)  sin  +  y^  cos  (4.6c) 

If  P  is  at  the  center  of  Slot  H,  then 

x4  a  *4  =  z4  3  0  (4.7) 

Substituting  Eq.  (4.7)  into  Eqs.  (4.6a,b,c),  we  have  for  the  (x,  y,  z) 
coordinates  at  the  tth  slot  center  [with  the  aid  of  the  condition  (4. Id)] 


**0  3  l*  +  b  cos  ^o}  cos  ho  3  R  cos  *£0 


(4.8a) 


y„A  «  (R  +  b  cos  )  sin  *  r  sin 


(4.8b) 


2to  *  b(1  +  sin  V 


(4.8c) 


Note  that  [see  Eqs.  (3.1),  (3.2),  (3.3)],  If  Ar,  *  x  ax.  +  y  Ay.  +  z  az0. 


r  *  *  6X  AX£  +  6y  **£  ±  l®zl  *Z£ 


a  Ax*(-Bx  sin  yiQ  +  By  cos  »‘Q) 


*  Ay;[-(BX  cos  <PiQ  +  By  sin  ^Q)  sin  ^  +  |BZ|  cos  ♦“„] 


(4.9) 


where  It  follows  from  Eq.  (4.8)  that 


ex*£  3  Sx*£  3  A  s1n  *£o  +  By  cos  *£0 


(4.10a) 


=  "(Bx  cos  *£0  +  6y  sin  *£o}  sin  *£o  i  l6zl  cos  ho  (4J0b) 


3  ±  J1  "  -  (&*)2 


(4.10c) 


The  exponential  factors  In  Eqs.  (3.5a,b)  can  be  written  [with  the  aid  of 
Eqs.  (4.8a,b,c)]  in  the  form 


i 


jt.n-i 


(4.11) 


g- j  kQ ( Sxx4o+0yyio^l ®z I 2  to ^ 

-  i  k0Sh  ( R*  bcos$£0 )  cos  ( <(>6 ) + j  kQ  I  e2 1  b(  1  +s  i  n<|>  "Q ) 


where 


6h  *  sy  1  *e  ■  tan'’  (  rx  > 

The  slot  center  coordinate  x^Q  appearing  in  Eqs.  (3.7a,b)  is  the  displace¬ 
ment  along  the  cable  of  the  slot  center  from  the  launching  point  for  wave 

energy  propagating  down  the  cable  [as  indicated  earlier,  this  point  will  be 
designated  as  the  origin  of  the  coordinate  system  defined  with  respect  to  the 
cable  itself,  i.e.,  the  (r",  0",  z")  system].  It  is  easy  to  deduce  (see 
Figure  4.2)  that  (if  subscript  to  indicates  slot  center  for  all  coordinates) 


'to  to 


to 


(4.12) 


The  other  slot  center  coordinate  appearing  in  Eqs.  (3.7a,b),  denoted  by 
y^Q,  is  the  linear  coordinate  corresponding  to  the  angular  position  of  the 
slot  center  around  the  preiphery  of  the  cable.  It  is  evident  from  Figure  4.2 
that  this  coordinate  is  given  by 


ylo a  So 

The  approximations  used  to  justify  the  modelling  of  the  slots  as  rec¬ 
tangular  (ignoring  the  curvature  of  the  slots  both  along  the  cable  and 
around  its  periphery)  and  thereby  justify  the  use  of  Eqs.  (4.12)  and  (4.13) 
are: 


4-6 


(1)  The  radius  R  is  so  large  and  the  slot  length  so  small  that  in  the 

vicinity  of  the  slot  we  can  Ignore  the  curvature  of  the  slot  along  the 

cable,  i.e.,  in  the  direction  of  wave  propagation. 

(2)  The  slot  width  is  small  enough  to  neglect  the  curvature  of  the 
slot  around  the  cable  periphery. 

***** 

The  forms  (3.4a,b,c,d)  used  to  represent  fields  are  valid  in  the  case  of 
propagation  along  a  coaxial  cable,  since  fields  can  be  represented  as  a  series 
of  propagating  modes,  where  the  ntn  mode  has  a  factor 


Fn(z\  ♦")  -  e 


e 


(4.14) 


where  k2n  is  the  complex  propagation  constant  of  mode  #n  and  n  is  the  mode 
integer. 

To  relate  Eq.  (4.14)  to  the  forms  (3.4a,b,c,d)  we  set  z"  and  <i>"  to  z"Q 
and  <(>£0  in  Eg.  (4.14)  and  then,  with  the  aid  of  Eq.  (4.12) 


where 


znR 


Jo. 


=  k 


znR 


^zn1 


znR 


k  =  — 
znR  vn 

vn  *  phase  velocity  of  nth  mode  in  meters/second 
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u  :  .it Lrnjj  1  .inn  o'  n  mode  in  ncpers/Kiotor 


?ni\  i; 


kznl  -an  k. 


from  Eqs.  (4.13),  (3.4a,h,c,d),  and  (4.14) 


i' 

'n  ”  F  b 
o 


(4.16) 


Ths  variable  z^  in  Eqs.  (3.4n,b,c,d)  is  the  radial  coordinate  r";  lienco,  we 

can  desi qna te  C  ,  C,  . . .  etc.  in  Eqs.  (3./ia,b,c,d)  as  functions  of  the  value 
.van  -vDn 

of  r'  at  the  Slut,  which  is  equal  to  b,  the  outer  radius  of  the  cable.  Thus, 
the  forms  of  Eqs.  (3.7a,b)  for  the  specialization  in  this  section  are  [using 
Eqs.  (4.11)  through  (4.16)]: 

-  ,  „  O  kf  T 

*  -4, . l e 

■*.ik0IS2|b(lnin,}J0)  -jk0(k2„»)4t0  in^o 


k  L  . 

O  C.  r«.t 


K  w.  ,  + 

l  u  r  1 


sine  [  2  +  kzn)]  sinc  [  ~T"  ( B 7 ' e.  ‘  (T‘b  }1 


(-20  £^>(b)  4  If  C<''(l.) 


(4.17a) 


o  n 


i  n 


H,  B, )1-  ■  -  >.u  -  ■  «• 

>  '  >1:  a  't  ,  '■ 

n 


+  ik  lb(l  i-sina"  )  -j!i  (k  •  );.  jtr" 

‘  o'  /.'  v  'to'  o'  zn  ’  1  . n  J  1  v.o 

...  O  0  -i' 


...  sine  [\-(^,£,kzn):  sine  lV"  k'b 


<Y0  Wb>  +  [S+  cdnJb)  -  (./  x  ccn:(b);i>  (4 

...  f  + 

whore  in  Eqs.  (4.1/a,b)  -.r,^  and  p^.^  are  given  by  Eos.  (4.1  Ha  »L' ) . 

From  Appendix  IV,  we  know  that  an  arbitrary  mod 3  will  have  the  goner; 
forms  [where  subscripts  (r",  q",  z")  denote  components  along  the  indicated 
directions  and  (f",  </' ,  z")  are  the  unit  base  vectors  in  the  double-primed 

■t  -  >  > 

cable  coordinate  system) 


Sane  ■  il  -  -fe 


w  *  4:;)<b> 


Sena  ■ 3».  -  s  Hz»’(b> 


Cdn<,  -  "r»’<b> 


where,  in  general. 


-In)  "jkokznz" 

E(b,  r.  /")  F(n,(b)  e  0  zn  •  e 


'(b,  <f>" ,  z")  -  'li(b)  <■ 


-jk  k,z"  .in-?" 


» *• 


and  where 


0  for  a  TE  mode 

T^?}  ■  0  for  a  TM  mode 
and 

E^  =  =0  for  a  TEM  mode 

It  is  easily  deduced  from  the  cable  geometry  that 

S  3  ^cos  ho  cos  ^o}  +  2(s1‘n  ho  cos  *£o}  +  |(sin  V 
♦t  *  *(-cos  ho  sin  ho}  +  i(~sin  ho  sin  ^o5  +  !(cos  ♦to) 

It  *  -5  S1‘n  ho  +  1 cos  ho 

Substituting  Eqs.  (4.19a,b,c)  into  Eqs.  (4.18a)  and  (4.18c),  the 

components  of  the  vectors  C,B  and  C  „  are: 

-►an  ♦cn 

Can£x  3  'sin  ^oTi")(b)  +  cos  *£o  Sin  ♦to  ^")(b> 

Can£y  3  cos  **o  ^")(b)  +  s1n  ♦to  sin  *£o  *T'){b) 

Can£z  3  -cos  ho 

Ccn£x  =  -sin  ho  ^")(b)  +  cos  ho  sin  ho  4")(b) 


(4.19a) 

(4.19b) 

(4.19c) 

x,  y,  z 

(4.20a) 

(4.20b) 

(4.20c) 

(4.21a) 
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w  -  *«,  i!"  <■„  *«»  Ko 


(4.2lc) 


From  Eqs.  (4.20a,b,c)  and  (4.213,6,0),  the  bracketed  expressions  in  Eqs. 
(4.1/a,b)  are  given  by 


sf  *K  Wb>  *  [f  cbns,(h)  -  (f  X  Can>(b)]l)  r* 


/N  "  +  S'  "  >  /v  + 

=  .?  rEx  +  y  FEy  +  *  FEz 


(4.22a) 


*  “(Yo  Wb)  +  t  Cdn^b>  '  tf  x  Wb)4 


A  ~4-  /\  ~4I  A  ~  -f- 

i  %  +  Z  Fly  *  ?  Flz 


(4.2?b) 


where 


FEx  =(Zo[sin  H1”)(b)  -  cos  ho  Sin  Hz")(b)] 


+  (3X  E^(b)  +  |8Z|  [cos  4>£0  E^(b)  +  sin  <*)?0  sin  <frj0 ‘^(b)] 


£.o . £0  r 


+  $  [cos  q"0  E'j^(b)])  f> 


(4.22ax) 1 


Fty  =(-Z0[c0S  V.  ^'')(b)  +  sin  *to  Sln 


+  ey  E^}(b)  -  (>x[cos  <y;0  E^(b)] 


-  'W-'"  -  cos  <V.o  sin  ho  LZU  {ir 


Atcos  ^0II^)(1))J+  |Pzl^!.°(i.) 


Uyf  ;ii1  ^f,o  E([")(b)  "  tos  Ho  sin  ho  4"}<W 


'  r'x[cor’  Ho  Ei")(b)  +  sin  ho  Sin  ho  Ez")(b)]) 
:(-Yo[sin  ho  E|"){b)  “  cos  ho  sin  *£o  Ez")(b)] 


+  Px  H^(b) 


+  |P.Z|[C05  (J)Jl0  fiJiJ?)(b)  +  sin  <j);o  sin  Hj,!?;(b)] 


i  n  A  **  ll(n)( 


k"  u(n) 


+  ®y[cos  ho  Hz 


"}(b  )])f4 


'o^cos  ho  +  sin  ho  sin  ho  E["^b^ 


►  Py  Hp!.^(b)  -  Px[cox  <j>£0  H^(h)] 


t  I  1  [sin  B  jj)  -  cos  4>£o  sin  ^H^fb)] 

=  (-Yo[cos  Ez")(b)]  ±  l3z[H^}(b) 


3y[sin  <f'£o  H^(b)  -  cos  <j>Pn  sin  H^"'(b)] 


lo  lo  " z' 


'  Bx[cos  ^0  HI"^b)  +  sin  ^£0  S1’n  *£0 


fr  Oil  Eqs .  (/1.17c:,b),  (4  a  ,x,y  ,7) '  and  (4.22i>x,y,/) '  we  have  t!.,-  fol¬ 
io:-:  in «  expression::  for  the  x,  y,  z  components  of  the  plano-vu-'  spectra  of 
the  fields  from  Slot  Hi  with  Node  £n  propagating  alone  the  cable: 


y  U  ko,?  T  •iW"'"’!/*#-*/ 

Z  "  7:i  " l-  6 

n 


[Sin!(?h»  ‘ 


TjgsjikOni^y  -:k0(K2n^)*to  j«*Jo 


Sine  f  (s±  t  {  >]  Sine  f  -°y'-  <«,■*  -  r  L  » 

J  0 


|<*!  r<;>(b) 

*  W 

a]3  tz„ 

21 

22 

23 

31 

32 

33 

41 

42 

43 

51 

52 

53 

61 

62 

63 

+  R  ^ 
“ll 

nj!?)(b)  ■ 

1-  cj^  H^}(b) 

+  M  M 

+  b13  i.z„ 

21 

22 

23 

31 

32 

33 

41 

42 

43 

51 

52 

53 

61 

62 

63 

- 

b41 

11 

IS) 

X 

(4.23a ,b) 


a(^)  =  jj(^)  -  p 

ft21  “51  Py 


If  <j>£0  Is  not  the  same  for  all  slots  and  If  the  slots  are  not  uniformly 
spaced  along  the  cable,  then  the  procedure  for  computation  is  the  programming 
of  Eqs.  (4.23a,b)  as  given  above  for  each  value  of  the  index  £,  followed  by 
accumulation  of  the  sum  over  l  from  l  *  1  to  £  *  N^,  where  is  the  total 
number  of  slots  along  the  cable.  If  this  technique  is  to  be  used,  then  taking 
advantage  of  various  approximations  that  apply,  such  as  Eqs.  (4.1a,b,c,d,e) , 
to  simplify  the  expressions  is  not  very  productive,  since  the  effect  of  these 
simplifications  on  computer  time  is  negligible.  A  program  has  been  written  in 
accordance  with  the  above  procedure  to  allow  for  the  possibility  that  we  might 
wish  to  evaluate  the  general  case  where  the  slots  are  not  all  of  the  same 
dimensions  and  not  all  uniformly  spaced  both  along  and  around  the  cable.  In  a 
practical  sense,  this  computation  might  involve  a  prohibitve  cost  in  computer 
time  because  N^,  the  number  of  slots,  might  number  in  the  thousands. 

For  those  reasons,  we  find  it  necessary  to  search  for  justifiable  simpli¬ 
fications  in  Eqs.  (4.23a,b).  Fortunately,  in  the  problem  of  imnediate  interest, 
the  cable's  slots  are  all  placed  at  the  same  angle  around  the  cable  and  are 
uniformly  spaced  along  the  cable.  Thus,  for  this  case 

<f>£0  *  *o  f°r  va^ues  *  (4.24) 

♦no  •  *1.  *<»-’>*  »0  •  *10  *  <4-25> 

where  Is  the  position  of  the  first  slot  center,  A  is  the  angular 
spacing  between  slot  centers,  and  since  the  slots  are  distributed  uniformly 


-  ?v  . 


along  the  out  i  re  cable,  fh  4  v 


1  ~  L  for  .i'll  9 


W  W  for  dll  9 


(4.2,) 

(4.27) 


It  is  shown  in  Appendix  V  that  over  the  rang?  of  parameter  values  of 
direct  interest  in  this  study 


k  L. 

:inc  [  — ?-  (px^  +  kzn)I  -•  1  for  all  £,  n 


(4. 28) 


and 


k  W,  ,  n  W 

sine  t  -  V  -  TV  »  *  (  -ar  > 


for  all  J l,  n 


(4.23) 


Using  Eqs.  (4.24)  through  (4.29)  in  Eqs.  (4.23a,b)  ive  obtain,  after 
summing  over  all  slots  (i.e.,  I  ~  1  to  Ns), 


[ll!n)(3j] 

L  +  \>h  x 


(j  k  r  +jkj8Jb(l+sin;.")  jn<i>" 


4ir 


O'  e  '"O'-Z 


•  c 


o  .  ,  nW  \ 

sme  (  2b  )  f 


-Jko(kzn")U0 


4-16 


1 


•  ** 


V 


where  (with  explicit  indication  of  argument  6^  deleted) 


[Einh  =  [8  E^  +  8  cos  4"  E^]  S<"> 

L  +  Jx  LPx  r"  yo  2  J  o 

Jd>lo 

+  ^-2—  [+|SZ|  E^}  +  j|Bz|  sin  ♦“  l! 

-  z0  sin  t*Q  H^]  sjn) 

-J*lo 

*  ■S-? —  Cil8r]  ■$’  *  |SZ|  sin  »jt<! 

-  Z0  sin  ^  H^]  s[n) 

[E(n)]y  =  [8y  E*?*  -  Bx  cos  <D“  E^}]  S<n) 

J*lo  t 

+  C+  j|8zl  ^  +  |6Z|  sin  **0t\ 


+  jZ0  sin  ♦-  H^]  s|n) 


-j<0 


+  C±  JlBJ  E<">  *  |6J  sin  ♦-  tM  .  Z  «<") 


'Z'  (j) 


o  z”  -o  r 


-  jZ0  sin  ♦"  h'IJ]  sfn) 


C+18,1  E^l?1  ♦  Z„  cos  ♦;  «>:>]  s<"> 


+  S-Jl  C-(ex  -  jsv)  eJI)  +  j(sx  -  j6v)  s|n) 


y'  <t> 


y,  -2« 


-j*lo 


+  ^“7—  [-(0*  +  JBy)  “  j(S*  +  j6J  ^  sfn) 


y'  z’ 


C0X  Hr"}  +  6y  cos  ♦;  Hi")]  Son) 


j* 


10  U(n)?j|07l  sin  ♦“  H<2>  ♦  jYn  E<!) 


•o  <r 


+  Yo  sin  4>”  E^}]  S<n) 


-J* 

♦  S—  um  i  j i b2i  sin  ♦;  "i"1  -  jy0 


+  Yq  sin  E^]  s[n) 


tK  H<!?>  -  6,  cos  ♦-  H<!)}]  $in) 


1 T -  W|BZI  Hr}  *  le2l  sin  ♦;  TT<!?>  +  Yo  t<!?) 

JY0  s1n  ♦;  ^"}3  s<n) 


C±J ( e2 1  hJ2)  ?  |6Z|  sin  ♦;  H^}  +  Yo  t*!J} 


+  jY0  sin  ♦-  T<?>]  S_(n) 


(4.30e) 


L"in\  3  t±j e2 1  Uj.!?)  -  yq  cos  *“  £<!?>]  s£n) 


*  ~  1 0 

1  [-(6x  -  «y)  ♦  J(sx  -  J6  )  H<:>]  S<"> 


f-(sx  +  V  -  Mx  *  ity)  h'I1]  sf">  C-30^' 


where 


2ir(JM) 


s(")  ,  f  e’jkoVcost<,,1o'^~lir“"^:i  ‘‘,ko(kzn*,t2Tr(JM)/NsJ 


(4.31a) 


2it(JM ) 


$(")  =  I  e’Jk°0h/?COS[01o'K^S - *8]e'J^(k2n*)t2*(*‘1)/N$] 


...  e 


j[2*U-l)/M 


(4.31b) 


a 


ye  3  0 

-1 


-j[^(f-D/r:,3  , 


(4.31c) 


To  evaluate  Eqs.  (4.31a,b,c)  in  cases  where  the  effect  of  the  factors  f’ 

can  be  approximated  as  unity,  we  can  invoke  the  well  known  Bessel  function  rela- 
A-4 

tionships 


gjxs  in 


jxcos(e-p 


w  . 

I  JD(x)  e3p 

p  =  -C*a  H 


(4.32a) 


and 


Jp(-x)  =  (-l)p  Jp(x) 


(4.32b) 


which,  when  substituted  into  Eqs.  (4.31a,b,c)  yields  after  interchanging  the 
order  of  the  summations. 


-(n)' 


:(n) 


-(n) 


I  H)P  Vko  3h  e  1  3  1 

P  =  -'-o  ^ 


S  jp[2ir(l-l)/Ns]-j[2nU-l  )/M,] 

e  •  fko<kZn  *>] 


[ko<kzn  *)  ♦  1] 


[k0(kz„  n)  -  1] 


(4.33) 
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•  V 


We  now  invoke  the  well-known  formula  for  the  truncated  geometric  series 


N-  Nc-1  ,  nS 

f  /-1  =  l  zVa  VrV  z*  1 

£*1  v=o 


Ns  if  2  =  1 


(4.34) 


Applying  Eq.  (4.34)  to  Eq.  (4.33)  and  noting  that  [as  indicated  below 
Eq.  (4.15)] 


k  -  kznR  "  jan  .  r  * 

kzn  ^znR  "  Jan 


(4.35) 


we  obtain 


pJL<-’>pVko‘b*>* 


Jp(*lo‘*8'7) 


’c(np)' 

o 

c(np) 


c(np) 


(4.36a) 


where 


c(np) 


j2n[p-kk^]  -2TO„S 

1  -  e  •  e 


1  -  e 


j(27T/N.)[p-kP 


1f  an  *  0 


*  N$  ^  P  *  ko  *<2nR  *  +  2irm,  where  m  is  an  integer  and  on  *  0 
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;<nP>  , 


J2it[p+l-k  k  pff]  -2™„S 

l  -  e  •  e _ 


1  -  e 


j(2ir/Ns)[p+l-k0kznR/?]  -2Tr<y?/N$ 


if  “n  *  0 


3  Nc  if  p  *  k^  k  D  /?  -  1  +  2irm  and  a  3  0 
bo  znR  n 


;(np) 


JMp-l-k  k  „»]  -2 ™n« 

1  -  e  *  e 


1  -  e 


j(2Tr/Ns)[p-l-k0kznRfl]  -2TOni?/Ns 


if  an  *  0 


3  Nc  if  p  3  k  k  Q  R  +  1  +  2m  and  a  3  0 
bo  znR  n 


Using  as  a  basis  the  approximations  justified  in  Appendix  V 
ill  assume  that  kQ  kznR  R,  which 
be  approximated  by  an  integer,  i.e.. 


we  will  assume  that  kQ  kznR  R,  which  is  a  very  large  positive  number,  can 


kQ  kznR  R  =»  q,  where  q  is  a  positive  integer  greater  than  20  (4.36b) 


Also  justified  in  Appendix  V  is  the  approximation 


-  2Tra_i?/N  - 
e  n  S  *  1 


(4.36c) 


and  the  fact  that 


an  f  0  in  all  cases  of  interest 


(4.36d) 
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At  this  point,  we  will  test  the  equations  (4.17a,b)  from  which  Eqs. 
(4.30a,b,c,d,e,f)  are  derived,  to  determine  whether  they  satisfy  the  Maxwell 
equations  in  their  present  forms.  As  noted  in  Appendix  III,  If  Eqs.  (4.16a,b) 
are  consistent  with  Eq.  (III. 6a),  that  is  sufficient  to  confirm  this  con¬ 
sistency.  Performing  the  operation  (Y  6^  x)  on  as  given  by  Eq.  (4.22a), 
we  obtain  (with  the  aid  of  the  vector  identity  ax(bxc)  =  b(a,c)  -  c(a,b),  and 

the  fact  that  by  definition  Y„  =  1  and  of  course  gix3±-30,  B*  •  =  1): 

0  0 


'of 


-  fV 0  f  '  +  ‘Wb)] 


(4.37a) 


Performing  the  operation  (-ZQ  S*-  x  . ..)  on  F^  as  given  by  Eq.  (4.22b),  we 
obtain 


-Z-  ^  *  S  *  !e  +  ■  C,„„(b))  -  Chn,(b)] 


cn£' 


'bn£v 


(4.37b) 


Performing  the  operation  (8^  •  . ..)  on  F>  and  Ft:,  we  obtain 


"Vf  •  ScnH<b»  ’  W>> 


f  '  $  *  Yot!±  '  W<b”  *  Wb> 


(4.37c) 

(4.37d) 


It  follows  from  Eqs.  (4.17a,b),  (III .6a,b,c,d)  and  (4.37a,b,c,d)  that 
the  conditions  required  for  consistency  with  the  Maxwell  equations  are 


z-Cb±  ‘  Wb)]  ■  cb».<b>  * 0 


bn£' 


(4.38a) 
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and 


Vf  •  *  Wb>  •  0  (4'38l>> 

For  the  particular  case  of  Interest,  Eqs.  (4.18a,b,c,d)  translate  Eqs.  (4.38a,b) 
into  the  form 

Zo  H^)«)(b)(B±  •  z")  -  Z0  •  $")  -  E^b)  =  0  (4.39a) 

Yo  ^(bKS*  •  z£)  -  Yq  E^S^b)^  •  |")  +  H^(b)  =  0  (4.39b) 

For  TEM  modes,  where  E^,  "E^,  all  vanish,  Eq.  (4.39a)  be¬ 

comes 


(bn<™> 


■  Z0CB<J>(b)]<TEM>(f  •  IJ) 


and  Eq.  (4.39b)  is  an  identity 


(4.40) 


For  TM  modes  In  the  special  case  where  the  boundary  conditions  require 
that  E^(b)  =  E^(b)  =  H^(b)  =  0  and  the  TM  nature  of  the  modes  requires 
that  =  0,  again  Eq.  (4.39a)  becomes 

EEjJ?,(b)](™)  *  Z0CH^)(b)3(™)(|±-  •  !")  (4.41) 

and  Eq.  (4.39b)  is  again  an  identity. 

For  TE  modes,  in  the  special  case  where  the  boundary  conditions 
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require  that 


^(b)  -  H^}(b)  «  0 

and  the  TE  nature  of  the  modes  requires  that  E^  *  0,  Eq.  (4.39a)  becomes 

ClJ?)(b)]TE  =*  ZQ[Hjlf)(b)]TE(fft  •  I”)  -  ZQ[H^)(b)]TE(g±  •  f)  (4.42) 
and  Eq.  (4.39b)  is  an  identity. 

If  0*  were  along  z",  this  would  imply  a  plane-wave  propagating  along  the 
cable. 

At  this  point  we  will  make  an  adjustment  in  the  expressions  (4.17a,b)  and 
the  subsequent  relationships  derived  therefrom.  In  effect,  we  are  forcing 
Eqs.  (4.17a,b)  into  consistency  with  the  Maxwell  equations  in  £-space  (i.e., 
those  of  a  plane-wave  with  propagation  vector  £) .  Another  way  of  looking  at 
this  is  to  consider  it  as  a  matching  of  the  fields  just  inside  the  slots  with 
the  fields  just  outside  the  slots.  The  latter  must  obey  the  conditions  (4.38a,b) 
which  in  turn  must  be  equivalent  to  the  vanishing  of  the  last  set  of  terms  in 
Eq.  (III. 7),  since  both  of  these  conditions  imply  consistency  with  the  £-space 
Maxwell  equations.  Such  consistency  is  required  of  the  plane-wave  spectrum  of 
a  field.  The  original  assumption  that  the  fields  on  the  slots  can  be  approxi¬ 
mated  by  those  of  the  wave  propagating  along  the  cable  as  if  the  slots  were 
not  present  is  being  replaced  by  the  assumption  that  there  is  a  small  per¬ 
turbation  due  to  the  slots.  This  perturbation  changes  the  fields  on  the  slots 
just  enough  to  force  the  plane-wave  spectrum  of  the  fields  at  an  arbitrary 
point  in  space  to  obey  the  Maxwell  equations  for  a  plane  wave  with  propagation 
vector  |.  The  reason  the  original  plane-wave  spectra  (4.17a,b)  did  not  show 
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x 


k 


consistency  with  these  equations  was  because  of  our  insistence  i  he  t  the  fields 


or,  the  slot  be  exactly  equal  to  those  propagu Liny  down  a  slotless  table  at  tin 
slit  position.  Doth  of  these  sets  of  plane-wave  spectra  (4.17c‘.,b)  and  those 
to  ba  given  in  what  follows  arc  rough  approximations  to  the  truth.  It  is  as¬ 
serted  that  the  latter  are  better  approximations  than  the  forcer. 

The  indicated  adjustment  in  the  original  forms  (4.17a,b)  consists  of  re¬ 
definition  of  fr  and  Fn  in  Eos.  (4.22a)  and  (4.22b)  respectively,  as  follows: 


o  icnt(b>  -  t  *  W’)  +  t  Ijg  ■  JcntC*”)^ 
£aot<b>  -  f  *  Ecn*<b>  -  f  Vf  '  Wb”)f± 


(4.43a) 


(4.43b) 


Applying  (Y  fr  x  ...)  to  as  given  in  Eq.  (4.43a)  yields 


+ 

0"  x 
-► 


;h 


(4.44a) 


Applying  (-Z  0*  x  . . . )  to  Ft;  as  given  in  Eq.  (4.43b)  yields 

0  ->  -vn 

-zofyS"Et  <4-<<b> 

Applying  (p+‘  •  ...)  to  F>  and  F^  yields 


FT+  •  Fr  =  0  (4.44c) 

■>  >L 

0*  •  FjJ  =  0  (4.44d) 

The  equations  (4.44a,b,c,d)  establish  the  consistency  of  the  forms  of  ^(3,  ) 


4-X,! 


*■(3.)  impl  ied  by  Eqs.  (4.43a,b)  with  the  Maxwell  equations  ( II 1 .6a,b,c,d) . 
>n 

The  quantities  replacing  F^  ....  F^z  in  Eqs.  (4.22a,b,c,d ,e,f ) 1  are: 

"(4")[0]  *  cos  *i0] 

*  e'">11Y  *Jo  *  iP'Z  Sin  ':!.0  Si"  "Jo5  * 

*  Sjf’Ro'O  '  sx>  sin  *u  4  8x  ry  cos  VJ 

*  H^’rz^-sin  *J0((1  -  S*)  cos  ho  -  Bx  Py  sin  *w) 

+  nx|n2|  COS  *J0)])f1  (4.45a) 

(F|y)1"  sin  <(0] 

*  E^"! f-Bx  cos  i';o  + 1 6Z |  cos  0,o  sin  4"0)  +  H^tOl 

*  -  By)  “4  *so  *  6>:  By  4-in  «;„») 

+  ^.’[-Zofsin  <.;o(0  -  Bp  sin  t)0  -  Bx  Py  cos 
1  |Bz|8y  cos  f*  (4.45b) 

(>t/n  -(4")t°]  4  4;)['(6x  cos  *10  4  6y  sin  ♦to” 

*  e^Csin  4»;0(-Bx  sin  ♦  By  cos  ♦„)] 
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+  V  [0.1  +  ll:..,C+Z0|fiz|(-!'ix  sin  0>o  >  0y  cos 

+  "l-^V-lPzl2  cos  *Io  -  IM  sin  cos  ho 


+  Ry  sin  V)])f± 


(4.45c) 


=(4")[0]  +  El")|:“Yo{(1  -  3X>  sin  ^0  +  3x  By  cos  W1 
+  f2")[Y0{si'n  *£0((1  '  $  cos  *£o  -  Bx  By  si"  V 
*  ex I P2 1  cos  *lo}]  +  “J^ro]  +  cos  ^0] 


+  lll")|:By  cos  ^£0  ±  l0zl  s1n  «£0  sii1  ^o])f± 


(4.45d) 


'(4"’[0]  *  4">[yol(]  *  /'  r-os  ho  *  \  ey  sin  »20,] 


+  E<;>[V0{sin  (,Jo({l  -  6y)  sin  -  Sx  ey  cos  *lo) 
*  Ry I ez I  cos  *;„>]  +  SjS’tO]  >  ii^CilS,!  sin  *[0] 


*  "2">['ex  cos  ho  i  IM  cos  ho  sin  ho1'/ 


(4.45e) 


-(e^}[0]  +  C^C+ISJ  Y0(-Sx  sin  0;q  +  3y  cos  4^)] 

+  cos  *"0  ±  |6Z|  sin  o;'0(3x  cos 

+  By  sin  ho™  +  *$^[0]  +  cos  «}n  +  K  sin  4.J] 


4”  L  'Kx  —  v£oy 


+  \n  +  i'v  cos  v£o):ij  f- 


(4.46:) 


The  case  actually  programed  for  the  computer  vw  >  that  wherein  4."  is 
constant  an.:  the  slots  arc  spaced  equally  along  the  cable,  resulting  in  the 
conditions  (4.1:4)  and  (4.25).  lor  this  case,  the  expressions  analogous  to 
figs.  (4.3.7a,b,c,d,e,f) '  [obtained  frc:r.  tqs.  (4.45a,b,c,d,c, f)]  are  as  fol¬ 
lows  (in  abbreviated  notation): 


[i|n)]x  =  sW[(lin))0]x  +  -J--  sJn)[(L-in))+:ix  -<  e -2-U’  s5n)[(tir,))_] 


'Jr 


10 


y 

2 


y 

2 


(4.46a) 


j  0  ■»  , j  -  J^IQ 

[ii<n)]x  »  s'nl[(ii<">)0]x  ♦ s  r  s'nlL(fi(n))+]x  *  s  -2-  s(n>[(fi'"!)J 


y 

2 


y 

2 


y 

2 


(4.46b) 


where 


[(E^)  ]  = 

LV  +  /oJx 


y 

z 


/  u( fl)  7 

2"  Z0 

+  F{n) 
2" 

■  fly~ 

A 

Ib2I2  _ 

0 

:os  <k" 
1 0 


A")  v 
V  Yo 

Xlf-zf 

+  h<S> 

'  By" 

?By!B2l 

-Bx 

-  Ib2I?  . 

0 

cos  4>" 
0 


(4.47a) 


(4.47b) 
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*  Hr’  Zo 

1 

C-4. 

X^ 

J 

4.  n(n)  2 

V  0 

-A*x 

-\ 

•jAy 

_+j|62|(0x+jBy)_ 

_±|Bzl(Bx+j8y)_ 

(4.47e) 


[(H^n))Jx  =  /  H^n  } 

'±|azl 

+ 

±j ! Sz I 

y  ) 

z  \ 

±j|3z| 

+  |BZI 

_-(8x+j  By)_ 

_-j(Bx+jBy)_ 

+  E^  Yn 
o 

-jA* 

+  V 

V  0 

A: 

A 

jA 

y 

y 

_ijlSzl(Bx+j8y)_ 

_+|8z|(6x+j8y)_ 

(4.47f ) 


Specializations  of  Eqs.  (4.46a,b,c,d,eff)  to  the  case  of  TEM,  TE  and  TM 
modes  are  presented  below.  These  are  based  on  the  forms  of  the  fields  for 
these  modes  developed  in  Appendix  IV. 
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TEH.  TM  Modes 


These  fields  are  given  by  Eqs.  (4.46a,b),  where 


(TEM) 

>  T"1 

yOJX 


y 

z 


0 


(TEH) 

)  T"' 

'+Jx 


y 

z 


z0  r. 


ii|0zl(8x-j6y) 


)Jx 


(TEH) 

(TMJ 


y 

z 


jA; 


L+j|8zl(Sx+j&y)J 


(TEM) 

(TM  ) 

)  1  ° 
'+ Jx 


y 

z 


Yw  t<5>(b) 


!lBzl 

♦J|BZ| 


-(Bx-J6y)j 


(TEM) 

)  i1™"* 
-Jx 


y 

z 


Yw  ^,?}(b) 


±le2( 
±J  I e2 1 


L-(sx+jsy) 


(4.48a) 


(4.48b) 


(4.48c) 


( 4 . 48d ) 


(4.48e) 


r 


^’(b)  ■  uji’lb)]™  .  for  TEH  nod. 

[from  Eq.  (IV.4a)]f  where  VQ  *  voltage  between  inner  boundary  at  r" 
outer  boundary  at  r"  ■  b. 


TM. 


tjj^b)  -  [^(b)]  n 


A™n  jk 


zn 


£n<ko  »> 


cn 


for  TE„  mode 
n 

[from  Eq.  (IV. lib)]  where  all  quantities  in  Eq.  (4.48g)  are  defined 
Appendix  IV. 


AX  *  (1  -  B‘)  *  J8X  By 


Ay  *  (1  -  By)  *  JBX  By 


Yw  *  wave  admittance 


For  TEM  mode 


Y  *  f~  \  e  *  complex  dielectric  constant 
w  "  vo  ca  of  cable  material 


For  TMn  mode 


TM  u  £ra 

y  s  y  n  «  ,  — 
w  w 


(4.48f) 

a  and 

(4.48g) 


(4.48h) 

( 4 . 48i ) 

( 4 . 48  j ) 

(4.48k) 

(4.48A) 
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For  TEn  modes: 


r(n)  (TE  )  ,  . 

>0lx  •  "iS’Cb)  Z0[-?6X|8 

y 

|6,|2 


C(H+n))03'TEn)  ■  H<5>(b» 


y 

z 


"By 


r(n)  (TE  ) 

«E+  U  0  " 


x 

y 

z 


-jAv 


-AJ 


+j|0z|(8x-jSy)_ 


&$><» 


y 

z 


JAV 


L+j|62|(Sx+j6y). 


Y^S^b) 


r(n)  (TE  ) 

C(H+  )+]x  " 


'x 

y 

z 


+  |6Z{ 
*Jl0zl 


■(Bx-jBy)j 


\n  ^(b)  -  JH<2J 


(4.49a) 

(4.49b) 
-  jH^(b)  sin  ^ 

(4.49c) 
+  jH^(b)  sin  ^ 

(4.49d) 

(b)  sin  ♦-) 


(4.49e) 


r(n)  (TE  ) 

Ct»+  )JX 


y 

z 


+  p. 


+3 


Y^n  E'Jj^b)  +  jH-^(b)  sin  ^ 


(4.49f ) 


TL'_ 


E^}(b) 


CE<ff)(b)]1E" 


jn  k 


zn 


TE. 


Y  ”n  N’(k  b)  k *  b 
w  n  cn  cn 


‘n  <kcn  b' 


(4.49g) 


TE 


TE 


h<j>(6)  •  rH^oor "  ■  f:r^3T  i.  <k,„  b) 


cn 


(4.49h) 


Y 


<T<g 


(4.491 ) 


Other  quantities  used  in  Eqs.  (4.49a-i)  are  defined  in  Eq.  (4.48h,i,j)  or  in 
Appendix  IV. 

The  forms  (4.48 a-£)  and  (4.49a-i)  are  those  actually  programed  and 

hence  on  which  our  numerical  results  are  based.  However,  in  the  computations  the 

(n) 

sums  over  the  cable  slots,  S^'  1  are  not  treated  through  Eqs.  (4.33)  through 

0 “ 

(4.36d).  These  expressions  would  be  strictly  valid  only  if  f*  were  approximately 
unity.  In  the  general  case  where  fA  is  accounted  for,  the  sums  are  replaced 
by  integrals  and  a  different  technique  is  used  for  their  approximate  evaluation. 
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4  V 


First,  wo  take  advantage  of  the  fact  that,  in  (4.31-a,b,c),  parameter 
values  in  all  cases  of  interest  are  such  that  [see  (4.36b)]: 


kA-R^I  ~  1 


(4.50) 


Equation  (4.50)  is  used  to  justify  the  approximation 

c(n)  ..  <-(n)..  «.(n) 
ao  5+  *  - 

From  (4.31-a,b,c)  and  (4.50),  where 

9  9  +  ki - 


(4.51) 


]  a 


(4.52) 


we  have 


s<">  -  s<B)-  s<">  ’■i-fvr 


(4.53) 


where  the  dependence  of  f4  on  <j>‘  is  explicitly  indicated. 
Noting  that 


kzn“  kznR  '  1  k 


(4.54) 


where  kznR  is  the  real  part  of  k2n  and  »n  is  the  attenuation  along  the 
cable,  we  can  express  (4.53)  gensrically  in  the  form 


s=  Sa  +  Sb 


where 


'a 

b 


1 

-3 


J  J 


i 


d9'ga  (9')e 
b 


(4.55) 


and  where 
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;(■:■')  ■  - 


9a(t  )  ■ 


rW  ♦  %  c“-(;  ' 

e  1  + - j 

R  J 


gb(v,}  = 


(  1-) 

V 


R  -  x  (o  cos  $  -  rv cos  ;>')  +  yjp  sin  $ -;(;sin  <*>')  +  2(7) 


R  =  +  p‘  +  ~  2  cos ( <!> '  -<f>  )  1 

jfR  6h[p  cos(* -«>,.)  -  ;/v  cos  (<p‘  <j> ^  )]  -  |^j  z] 
yVfi+  f>?  +  z?  ‘  2-$p  cos(i>'  -(TT1 
Since  ga(<|>')  and  gb(V)  are  both  real,  and  since 

k  jl  »  1 


(4.5G) 


for  all  cases  of  interest,  the  integrals  Sfl  and  Sfa  are  both  of  a  type  that 
can  be  evaluated  approximately  by  an  asymptotic  form  given  in  Reference  A. 5, 
pages  276-278.  The  method  is  based  on  a  form  of  the  principal  of  stationary 
phase,  but  the  case  involved  is  that  wherein  there  is  no  stationary  point, 
i.e. , 


/*• 


-77  W)  t  0  at  all  points 

fly 


Since 


"  \  Sln  <♦'  -+s>  -  kznR 


(4.57) 


(4.58) 


and  since  kznR  >  1  in  the  cases  under  consideration,  and  | sin  ($'-$  )  |  <  1, 


it  follows  that  (4.57)  must  hold.  Since  (4.55)  also  holds  in  ail  cases 
of  interest,  the  method  is  applicable. 


The  asymptotic  solution  is  [Reference  A. 5,  page  278,  Eq.  (6.57)]: 


(4.59) 


where  upper  and  lower  limits  4.'.  and  V  for  S  and  S.  are  dependent  on 

U  L  d  D 

the  forms  of  g  ($')  and  gbU')  respectively. 

R 

Another  term  was  included  in  the  computations  to  account  for  tha 
term  present  at  points  very  close  to  the  cable.  This  term  was  also 
evaluated  by  the  above  methoo 


5.  PLANE-WAVE  SPECTRAL  REPRESENTATION  OF  GROUNO-REFLECTED  FIELDS 


We  can  do  the  ground  reflection  problem  for  a  particular  spectral  compo¬ 
nent  (l.e.,  particular  values  of  8  and  8 )  and  having  accounted  for  the 

a  y 

boundary  conditions  for  each  such  component,  we  can  then  Integrate  over  all 
8  ,  8  .  This,  of  course,  reduces  the  problem  to  the  classical  case  of  plane- 

A  Jr 

A-l 

wave  reflection  from  an  Infinite  boundary  between  two  media. 

We  will  now  outline  the  method  of  analysis  which  was  carried  out  to  obtain 
the  (x,y,z)  components  of  the  reflected  wave  fields  and  to  construct  a  matrix 
relating  the  reflected  fields  to  the  Incident  fields.  Since  the  literature 
is  replete  with  solutions  of  this  problem  there  is  certainly  nothing  new  or 
original  about  this  portion  of  the  analysis. 

Consider  incident  and  reflected  wave  electric  and  magnetic  field  compo¬ 
nents  for  a  particular  value  of  (8  ,  8  )  where  incident  fields  are  denoted 

x  y 

by  (Ej ,  Hj )  and  reflected  fields  by  (f^.,  i^.).  We  must  also  include  the  fields 
of  the  wave  transmitted  into  the  ground,  denoted  by  (Et,  Ht). 

For  fixed  values  of  (Sy,  SJ  (see  Appendix  II)  the  horizontal  components 
of  the  electric  field  above  ground  (In  free  space)  are: 


+  E_ 


rx 


y  y 

z  z 


jk(Sxx+6vy)  ,  -jk|8_|z  ,  jk|6  | 

6  *  *  ♦  Erx  6  3 

y  y 

Z  Z 


(5.1a) 


(where  It  Is  noted  that  the  Incident  and  reflected  waves  propagate  downward 
and  upward  respectively). 

The  horizontal  components  of  the  electric  field  In  the  ground  (z  <  0)  Is 
(note  downward  propagation  only) 


jk(8xx+8^)  -jkr 
e  *  e 


(5.1b) 
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TEST-1 


where 


Y2  ■  /"2  -  («x  +  *y) 

and  v  is  the  complex  refractive  index  of  the  earth,  given  by 

v  *  vR  +  jvj  (5.2) 

where 


v 


R 


and  where  eR  =  relative  permittivity  of  ground,  tQ  *  permittivity  of  free  space 
-12  -7 

*  8.854(10  )  farads/meter,  uQ  *  4  (10  )  henries/meter  *  Magnetic  permeability 
of  free  space,  assumed  to  be  also  that  of  the  earth,  o  *  conductivity  of  ground 
in  mhos/meter.  The  horizontal  magnetic  field  components  are: 


z  >  0 


/C  jk(Bxx+Byy)  -jk|B |z  .  ,  -  . 

Hix  +  Hrx  *ypQ  e  ^8y  Eiz  +  ^z^  Eiy^ 


jk|B.|z 

+  e  (By  Erz  -  |BZI  Ery)D 


(5.3a) 
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ft~  jk(e  x+6  v)  -jkje  | z 

Hiy  +  HW^e  ^c‘e  ('ezlEix +  6x  ei2> 


Jk | S  1 2 


+  e  ^8z^Erx  "  8x  Erz^ 


(5.3b) 


z  <  0 


/e 7  jk($  x+8  v)-jky z 

Htx  *//e  X  ^  *  <‘y  *tz  +  *z  V 


(5.3c) 


jk(Sxx+Byy)-jky2 


('Yz  Etx  *  8x  Etz^ 


(5.3d) 


To  determine  the  vertical  components  of  the  electric  field,  we  invoke  the 


equations 


7  '  li  a  Jk<8x  Eix  +  6y  E1y  *  I 6z I  Eiz^  *  0 


V  *  Sta  jk(0x  Etx  +  8y  Ety  "  Yz  Etz^  ® 


(5.4a) 


(5.4b) 


Solving  Eqs.  (5.4a,b)  for  Eiz’  Erz’  Etz  1n  terms  of  E1x>  Erx>  Etx  and  equating 

y  y  y 

the  tangential  components  of  both  electric  and  magnetic  fields  at  z  3  0,  we 
arrive  at  the  set  of  equations: 


(|8z!yz)  EfX  ■  (Ibz|y2)  Etx  3  -(|sz|y2)  Eix 


(|ez|Yz)  Ery  -  (|6zt  z)  Ety  3  -(|BZ|YZ)  E1y 


(5.5b) 


<8x8yV  Kx  *  «'  -  V  Sry  *  <8x8yl8zl>Etx  *  ^ -8x>l«zl>  Ety 


(SxSyY2)  t^x  +  (0  ■  Bx)  Yz^  Eiy 


(5.5c) 


«1  -  By)  y2)  8rx  *  <Wy\>  gry  *  '  6y))Bz[)  E(x  *  (6x6y|62|)  E(y 


((1  -6  y)yz)  Ejx  +  (8x8yYz^  Eiy 


(5.5d) 


Solving  Eqs.  (5.5a,b,c,d)  by  Kramer's  rule  yields  the  horizontal  compo¬ 
nents  of  the  incident  wave,  as  follows: 


E  *  2(  ^  ) 

Erx+  2(  y,  +  |8_|  ‘ 


1 


Z T  'BZ'  •  (1  -  bz  +~|e~|Y2>  lL"x " 2 


il»l  -  7  (1  -  &l  +  |BZI Yz)]  Eix 


+  l-8x8y^  Eiy-^ 


(5.6a) 


Y,  -  |6,|  i 

Srv+  =  2(  -z  -|5Z-]  )  - g-1 -  {[6x6y]  E. 

ry+  Yz+I8zl  (1  -  S2  +  |$z|yz)  Xy  lx 


+  0  ’  4  +  !6zlYz)]  v> 


(5.6b) 


where  the  minus  sign  on  Eix_  indicates  downward  propagation  of  the  incident 

y 

wave  and  the  plus  sign  on  lrx+  Indicated  upward  propagation.  Invoking  Eq. 

y 

(5.4a)  for  both  incident  and  reflected  wave  fields  we  obtain  the  vertical 

component  of  E  .  from  Eq.  (5.6a,b) 

T 
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We  can  easily  determine  the  magnetic  field  components  from  (5.6a,b,c): 


Hrx+  MJu0  ^By  Erz+  -  l6zl  Ery+^ 


fe  y  —  1 6  j  i  „ 

hr  (  Tgz i  ) - r - {(-2|e,|8x8v)  i. 

jUo  Yz  *  lBz>  (1  -  6‘  +  |82|y2)  2  xy  ix’ 


+  (-2|b2|CbJ-^(i  -  4+  !s2Iy2)3)  E1y. 


+  <-6y(l  -  82  -  |82|y2))  Eiz.} 


(5.7a) 


Hry+  J  vQ  ^Sz^rx+  "  ®x  Erz+^ 


Fo  /  YZ  ^ 

/m0 


z1  (1  -  S2  +  |8z|y2) 


•  ((2|bz|[62  -  '7  (1  -  6*  *  |e2|yz)])  E1x  *  (2|Bz|BxBy)  i,y. 


+  (-exO  -  s2  -  |82|y2))  E1z.} 


(5.7b) 


Hrz+  * JuQ  ^®x  E ry+  "  ®y  Erx+^ 


«  S  (  II — 


le 


TeTT  )(8y  Eix-‘  3x  Eiy-} 
(5.7c) 


To  compare  Eqs.  (5.6a,b,c)  and  (5.7a,b,c)  with  well-known  results  on 
the  reflection  coefficients  of  plane  waves  Incident  on  an  Interface  between 


Consider  first  the  case  of  horizontal  polarization,  implying  that 


where  y2 


Jv2  -  sin2  6g 


The  results  (5.10a,b,c)  are  consistent  with  well-known  results  on  re¬ 
flection  coefficients  with  horizontal  polarization  [e.g.,  Stratton,  p.  493]. 

From  Eqs.  (S.9a,b,c)  and  the  Maxwell  equations,  the  Incident  wave  mag¬ 
netic  field  components  are 


(implying  that  E1y 


-  0 


(5.11a) 


(5.11b) 


(5.11c) 


(Implying,  as  does  Eq.  (5.11a),  that  E^y_  ■  H^  )  where  H^_  Is  the  complex 
amplitude  of  the  magnetic  field  vector  of  the  incident  plane  wave. 

From  Eqs.  (5.9a,b,c),  (5.11a,b,c)  and  (5.7a,b,c) 


\  -  I  cos  80| 
\  +  I  cos  0g| 


)|cos  eg| 


(5.12a) 


Unct-  o 

fi1- 


(5.12b) 
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(5.14c) 


1 


From  Eqs.  (5.9a,b)  and  the  Maxwell  equations 


Hix_  *  0  (5.15a) 

V  C-f0lleix-  •  6x  !1l->  ■  St-tcos2  »6  *  s1"2  V 


(implying  that 


-  0 


(5.15b) 


(5.15c) 


Substitution  of  Eqs.  (5.13a,b,c)  and  (5.15a,b,c)  Into  Eqs.  (5.7a,b,c) 
yields 


0 


!Vyi 


Yz  -  |cos  eB|  sin2  e6  -  |cos  86!yz 

Yz  +  lcos  V  sin2  9g  +  |cos  0g|yz 


,  v2lcos  9gl  ~  Yz 
v2 1 cos  e  g |  +  yz 


H 


1- 


0 


(5.16a) 


(5.16b) 
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The  results  (5.16a,b,c)  are  consistent  with  well-known  results  on  re¬ 
flection  coefficient  with  vertical  polarization  [e.g.,  Stratton,  p.  494]. 

It  is  sometimes  convenient  to  write  the  results  (5.6a,b,c)  and  (5.7a,b,c) 
in  terms  of  the  angles  (0g,  #g),  i.e.. 


yz  -  [cos  eBi 


"rx+ 


Yz  + 


COS  8 


) 


81  [sin  eg  +  | cos  0g | yz3 


{E.x_[sin2  8g(cos  24>g) 


-  Icos  ©bIy23  +  E.y_  -  [sin2  8g  sin  2$g] 


(5.17a) 


Yz  -  icos  eg| 


Ery+  =  ^  Y,  +  | COS  07T  '  r,.->2 


8‘  [sin  0g  +  (cos  9g|Yz] 


(Eix_[sin2  9g(sin  2<j>g)] 


+  Ety.Cstrr  8a(-cos  2<J>fi)  -  |cos  9r|y7]) 


e"z- 


(5.17b) 


where  Y, 


=  j/v2  -  sin2  0. 


y2  -  |cos  eJ  sin2  0q  -  | cos  e  |y. 
Erzt  ■  -sin  «„(  o-V  |m.  /l  )(  - - SLl- 


8  +  l“s  9bT  "  sin2  e3  *  |cos  es|y! 


)  Eiz.  (5.17c) 


where  yz  *  ^v2  -  sin2 


5#  Yz  '  icos  V 


rx+  >o  (  Yz  +  'cos  0SI  }  (sin2  0fi  +  Icos  9g|Yz) 


{(-  |cos  0g|  sin2  0g  sin  2$g)  Eix_  +  (-|cos  9g|[sin2  8g(-cos  2<pg) 


:_2 


Icos  9g|Yzl)  Eiy_  +  (-sin  0g  sin  4>g[sin2  eg  -  |cos  0g|Yz])  Ei2) 

(5.17d) 
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W  *  -I'gtBy  -  Bx  ■  I’z^  ■  -|S2|6e22 


^13  ‘  ‘6yf6h  "  ‘  "By  ^E33 

«H21  ■  lBzltBx-8y-  !BzM  *  !Bzl  *E11 

®H22  *  2!Bz'Sx6y  *  6Z'^E12 

®H23  =  ”Bx^Bh  ‘  lBzlYz^  ’  "Bx  ®E33 
*H31  "  By^Bh  *  lBz'Yz] 

*mz  ■  -BxtBh  *  lBzM 


"tM*0 


(5.19) 

(cont'd) 
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6.  FIELDS  INCIDENT  ON  THE  SCATTERER-COORDINATE  TRANSFORMATION 
AT  INPUT  TO  SCATTERING  PROCESS 


The  scattering  theory  that  we  are  using  confines  Itself  to  the  scattering 
of  a  plane-wave  by  a  body.  The  spectral  fields  as  given  by  the  superposition 
of  Eqs.  (4.16a,b)  (the  fields  associated  with  the  direct  wave  from  the  cable 
slots)  and  Eqs.  (5.18a,b)  (the  fields  associated  with  the  ground-reflected 
wave  from  the  slots)  are  the  fields  of  a  plane-wave  propagating  in  the  di¬ 
rection  of  the  wave  vector  S+.  Denoting  these  plane-wave  fields  by  E.(3+) 

-+*  I  -H 

and  H^(g+).  the  next  phase  of  the  solution  to  our  problem  is  to  evaluate  the 
scattered  fields  when  the  incident  fields  are  E. (8+)  and  H.(8+). 

Before  we  can  do  the  scattering  problem,  we  must  perform  certain  coordi¬ 
nate  transformations. 

The  Barber  scattering  program,  which  will  be  discussed  in  Section  7, 
contains  two  coordinate  systems.  One  is  the  "lab  frame"  and  the  other  is  the 
"body  frame."  The  diagram  of  Figure  6.1,  due  to  Dr.  Peter  Barber,  illustrates 
the  two  systems. 

First,  note  that  Figure  6.1(a)  and  (b)  which  illustrate  the  "lab  frame" 
(whose  coordinate  axes  are  denoted  by  x^,  yL,  zL)  and  the  direction  of  the 
incident  wave  with  respect  to  the  lab  frame.  ,  The  incident  wave  travels  in  the 

j 

+z.  direction  in  the  lab  frame  (hence,  8+  ly'ln  the  +z.  direction).  The 

L  -*  L 

electric  field  vector  of  the  incident  wave,  denoted  by  E,. (8+),  lies  in  the 

-►1  -► 

(xL  -  yL)  plane  and  is  at  an  angle  with  respect  to  the  x^  axis  in  a  clock¬ 
wise  direction,  looking  along  the  axis. 

The  "body  frame"  [Figure  6.1(c)]  is  characterized  by  a  set  of  rectangular 
coordinates  (xg,  yg,  Zg).  The  scatterer  is  pictured  as  a  spheroid  with  its 
axis  of  symmetry  along  the  +Zg  axis.  (This  is  only  for  graphical  convenience; 


the  theory  does  not  require  that  the  scatterer  be  a  spheroid  oriented  as 
shown.)  The  direction  of  the  +Zg  axis  In  the  lab  frame  corresponds  to  two 
spherical  coordinate  angles,  a  polar  angle  Op  and  an  azimuthal  angle  4>p,  as 
shown  In  Figure  6.1(a). 

Our  first  task  is  to  develop  a  transformation  between  the  body  frame  and 
lab  frame  coordinates. 

A  simple  way  to  do  this  makes  use  of  the  diagrams  of  Figure  6.2(a),  (b), 
(c),  and  the  accompanying  equations  (6.1a,b,c),  (6.2a,b,c)  and  (6.3a,b,c) 
appearing  in  the  figure. 

From  Eqs.  (6.1a,b,c),  (6.2a,b,c)  and  (6.3a,b,c),  the  body  frame  coordi¬ 
nates  of  a  point  (xg,  yg,  Zg)  In  terms  of  the  lab  frame  coordinates  of  that 
point  (xL,  yL,  zL),  expressed  in  vector-matrix  form,  are: 

[rg]  *  [Mbl]0j_3  (6-4a) 

and  the  inverse  form 

(6-4b) 


where 


*L  *» 


Point  P  has  (x^  ,yL‘  ,zL' )  coordinates 
XjJ  *xLcos*p  +  yLsin<>p  (6.1. a) 

yl  “•\sin*p  +  -Ucos^p  (6.1.b) 

zL’  *zL  (S.l.c) 


Coordinate  system  (xL  yL.zL) 

Point  P.  has  coordinates  (xL,yL,zL) 
Zg  Is  unit  vector  In  direction 


Rotate  around  z^  axis  through  angle 

♦p,  such  that  In  new  system 

(*L  »yL'.zL').  Zg  lies  1n(xL'-zL')  plane 


Figure  6-2.  Transformation  between  lab  and  body  frame 


r  r 


Point  P  has  (xL",  yL",  z“)  coordinates: 


xJ«xL‘cos0p  -  2L’sinep 
*lm*l 


(6. 2. a) 
(6.2.b) 


•  xL‘s1n0p  +  zL'c°s0p  (6.2.c) 


V  ZB^ 


Rotate  around  yL'  axis  through  angle 
0p,  such  that  +zL"  and  +zg  directions 
are  parallel . 


Rotate  around  z^'  axis  through  angle 

counter  clockwise  looking  -zg  direction. 
Point  P  has  Xg.yg.2g  (body  frame 
coordinates) 

XB*  V005^  +  (6. 3. a) 

yg»  xL"sin^p  +  yL"cos<j»F  (6.3.b) 


W 


(6.3.C) 


Figure  6-2.  Transformation  between  lab  and  body  frame  (cont'd.) 


6-4-a 


i 


(cos  0„  COS  A  COS 

P  P  r 

(cos  6p  sin  (jip  cos  $F 

(-sin  0p  cos  <pF) 

-  sin  ♦  sin  <j>F) 

+  cos  $p  sin  4»p) 

(-cos  e„  cos  a  sin  <pc 

P  p  r 

(-cos  6p  sin  4>p  sin  $p 

(sin  ep  sin  $p) 

-  sin  $p  cos  <|>F) 

+  cos  $p  cos  4>p) 

(sin  ep  cos  «j)p) 

(sin  sin  <f>  ) 

P  P 

(cos  ep) 

(6.4c)' 

(cos  6p  cos  4>p  cos  <J)p 

(-cos  0p  cos  <J>p  sin  0F 

(sin  6p  cos  <j>p) 

-  sin  4>p  sin  <j>F) 

-  sin  <pp  cos  $p) 

(cos  9p  sin  <(ip  cos  tpF 

(-cos  0p  sin  4>p  sin  <f>p 

(sin  6p  sin  <j>p) 

+  cos  4>p  sin  $p) 

+  cos  <^p  cos  <t> p) 

(-sin  0  cos  $c) 

(sin  0_  sin 

(cos  9_) 

We  can  use  Eqs.  (6.4a,b)  to  obtain  a  vector-matrix  relationship  between 
the  unit  base  vectors  in  the  two  coordinate  systems.  Note  that  the  elements 
of  [rg]  and  [r^]  are  coordinates  of  an  arbitrary  point  in  the  body  and  lab 
coordinate  systems  respectively.  The  unit  base  vector  along  the  Xg  direction 
(for  example),  denoted  by  Xg,  terminates  at  a  point  P  which  has  coordinates 
(1,  0,  0)  in  the  body  frame.  In  this  case,  Eq.  (6.46)  would  read 


XL 

(MLB^11 

^MLB )l 2 

^  MLB  ^13 

1 

yL 

= 

(Mlb)21 

(mlb)22 

(MLB}23 

0 

ZL 

(MLB*31 

^MLB^32 

*MLB*33 

0 

wm  «. 

■» 

^ehi 


(mlb)2i 


(MLBJ31 


(6.5a) 


By  the  same  reasoning  (for  another  example)  the  unit  vector  in  the  yL  direction, 
denoted  by  y^,  terminates  at  a  point  P  with  coordinates  (0,  1,  0)  in  the  lab 
frame.  Thus,  Eq.  (6.4a)  in  this  case  would  read 


*B 

<Wll 

(MBL)12 

^Wl3 

0 

CsJ 

r— 

^B 

* 

(mBl)21 

^L^22 

^MBL^23 

1 

* 

(mbl)22 

(6.5b) 

ZB 

^  MBL  ^  31 

^MBL^33 

^MBL^33 

0 

(MBL}32 

From  Eqs.  (6.5a,b)  and  the  fact  [evident  from  Eqs.  (6.4c,d)']  that 
C'Wjk  “  ^B^kJ*  H  fol1ows  that 


h  *  (XL  *L  +  yL  &  +  ZL  iL}xB=l  *  ^ MLB 5 1 1  *1  +  (MLB>21 


W 


+  (MLBJ31  h. 
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4 


KMBL>i1*,Wl2(,Wl3^ 


-J  ' 

a 


& 


ih. 


(6.6a) 


and 


i.  *  (x»  a  *  j'b  h  *  h  !b>Vv«  •  (v>,2  «b  *  <n,L>22  a 

V1 


^L^32  §B  *  ^\b^21  ^\b^22^\b^23^ 


l& 


(6.6b) 


Arguments  like  those  leading  to  Eqs.  (6.5a,b)  and  (6.6a,b)  applied  to 
all  the  unit  vectors  in  both  systems  yield  the  following  two  vector-matrix 
equations: 


[Ug]  ■ 


(6.7a) 


and  its  Inverse  form 


(6.7b) 


where  [u^]  and  [tig]  are  "vectors  of  unit  vectors"  of  the  form 
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(6.7a)' 


(6.7b) * 


Cu,  3 


2l 


4 


’2b 

& 


We  now  wish  to  obtain  (1)  a  transformation  between  the  lab  frame  (xL,  y^, 
zL)  and  a  coordinate  system  to  be  called  the  "ground  frame"  (xQ,  yQ,  zQ),  ob¬ 
tained  from  the  basic  (x,  y,  z)  coordinate  system  by  a  translation  such  that 
the  origin  of  the  ground  frame  is  at  the  center  of  the  scatterer;  and  (2)  a 
transformation  between  the  ground  frame  and  the  body  frame  (Xg,  yg,  zg).  To 
accomplish  this,  we  first  note  that  Figures  6.3  and  6.4  contain  diagrams  which 
are  both  the  equivalent  of  those  in  Figure  6.2  except  for  changes  in  coordinate 
names. 

From  Figures  6.3  and  6.4  and  accompanyi ng  sets  of  equations  (6.6a,b,c) 
through  (6.13a,b,c),  by  analogy  with  Eqs.  (6.4(a,b)  and  (6.5a,b),  we  obtain 
the  following  sets  of  vector-matrix  transformation  equations,  analogous  to 
Eqs.  (6.4a,b)  and  (6.5a,b): 


C^l  *  C"b0K;0] 

(6.14a) 

[r0]  • 

(6.14b) 

■  ["80%] 

(6.15a) 

t",!  ■  CmobH=b3 

(6.15b) 

\ 

) 

/ 
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Analogous  to  Fig.  6. 2. a,  xL-xQ 

zrzo 
V*  b 
V6b 


Rotate  around  z0  axis  through  angle 


so  Zg  lies  in  (x0*  -yj) 

plane: 

x0‘  -x0cos^b  +  y0sin«|,b 

(6. 8. a) 

yQ*  ^x0sin4,b  +  y0cos<(,b 

(6 . 8 . b ) 

z0  “  z0 

(6.8.c) 

Analogous  to  Figure  6.2. 

b 

and  Eqs.  (6.1.a,b,c) 

Rotate  around  y’  axis  through  angle 
®b  50  zo"H  ZB 


Figure  6-3.  Transformations  between  ground 
frame  and  body  frame 


(c.2) 


V  "xocoseb*zos1n0b 


(6. 9. a) 


V  *yo’ 

V  *xo‘  s1n6b  +  zo'  coseb 

Analogous  to  Fig.  (6.2.c) 
and  Eqs.  (6.2.a,b,c) 

Rotate  around  Zg"  axis  through  angle 
CCW  looking  in  direction. 


XB"X0  cos*H  +  VSln*H 

(6. 10. a) 

yB*-xo,,sinV  yo"  cos*h 

(6. lO.b) 

ZB  *zo" 

(6.10-c) 

Analogous  to  Fig.  (6.2.d) 

and  Eqs.  (6.3.a,b,c) 

J 

(6.9. b) 
(6.9.C) 


Figure  6-3. 


Transformations  between  ground  frame 
and  body  frame  (cont'd.) 


6-9-a 


Analogous  to  Figure  6.2. a 

V  x0  XB"*’XL 

yCy0 

V  Z0  zb""zL 
V9s 

V 

(Note:  8  is  along  +  z^  direction) 


y 


"•'tate  around  z  axis  through  angle 


so  zjies  in  (x^-Zg)  plane 


x2'  *xocos0g +yQsin0B  (6. 11. a) 

y2'« -xQsin0g +y0cos<>s  (6.11  -b) 

-zQ  (6.11.c) 

Analogous  to  Fig.  (6.2. b)  - 

and  Eqs.  (6.1 .a,b,c) 


Figure  6-4.  Transformation  between  ground  frame  and  lab  frame 


fill  -  w 
W  •  t«oL3CXL3 
C^]  •  [\0][S0] 
C=0]  ■  [moL%] 

where 


(6.16a) 

(6.16b) 

(6.17a) 

(6.17b) 


(6.17c)1 


where  xs,  y$,  z$  are  the  basic  (x,  y,  z)  coordinates  of  the  scatterer 
center. 


(6.17d) ' 
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(cos  eb  cos  $b  cos  4>h 

(cos  efa  sin  0b  cos 

(-sin  0b  cos  <j>H) 

-  sin  $b  sin 

+  cos  <f»b  sin  $H) 

(-cos  eb  cos  $b  sin  <f>H 

(-cos  6b  sin  <(>b  sin 

(sin  0b  sin  4>H) 

-  sin  <(>b  cos  <*>H) 

+  cos  ^  cos  4>h) 

(sin  9b  cos  *b) 

(sin  6b  sin  $b) 

(cos  eb) 

(6.15a)' 

(cos  eb  cos  0b  cos  4>h 

(-cos  0b  cos  <frb  sin  <t>H 

(sin  eb  cos  (»»b) 

-sin  <|>b  sin  <J»H) 

-sin  <|)b  cos  <(>H) 

(cos  0b  sin  4>b  cos 

(-cos  6b  sin  <(»b  sin  <&H 

(sin  0 b  sin  <t>b) 

+  cos  $b  sin  $H) 

+  cos  4>b  cos  $H) 

(-sin  0b  cos  <j>H) 

(sin  0b  sin  $H) 

.O 

<x> 

1/1 

o 

o 

(cos  eg  cos  <pQ  cos  (cos  0g  sin  $g  cos  (-sin  eg  cos  $j) 


sin  <j>g  sin  <|>j) 


+  cos  $g  sin  ^j) 


(-cos  eg  cos  $g  sin  (-cos  0g  sin  $g  sin  (sin  0g  sin  <|ij) 


sin  $.g  cos 


+  cos  $.g  cos  4>j) 


(sin  0g  cos  <j>g)  (sin  0g  sin  <|>g)  (cos  0b) 


( 6 . 1 7a )  ‘ 


(cos  0g  cos  $g  cos  (-cos  0g  cos  sin  (sin  9g  cos  <j>g 


sin  *g  sin  4>j) 


-  sin  4»g  cos  $j) 


(cos  0g  sin  <frg  cos  4>j  (-cos  0g  sin  $g  sin  <pj  (sin  0g  sin  <f>g) 


+  cos  <|>g  sin  <pj)  +  cos  $g  cos  4>j ) 


(-sin  0g  cos  4>j )  (sin  0g  sin  ifrj) 


(cos  0g) 


(6.17b)1 
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From  the  viewpoint  of  programing  these  calculations,  we  must  distinguish 
between  the  variables  that  are  specified  at  the  beginning  of  the  program  ("given" 
variables)  and  those  that  are  calculated  from  the  given  variables  ("derived" 
variables).  Of  those  variables  appearing  In  [Mg03  or  (M^L  0b  and  #b  (the 
angles  specifying  the  orientation  of  the  scattering  body)  are  given,  while  <j»H 
Is  derived.  Of  the  variables  appearing  in  [M^0]  or  [M^].  0^  and  <j>g  (direction 
angles  of  6)  are  given,  while  Is  derived.  Finally,  of  the  variables  In 
[MgL]  or  [Myj],  8p,  $p  and  <j>F  are  all  derived. 

To  relate  the  matrixes  [Mg03  (or  [M^]),  and  [MQ^]  (or  CMj_0D ) >  we  form 
the  matrix  equations  [from  Eqs.  (6.7a,b),  (6.15a,b)  and  (6.17a,b)]r 

CSfi]  -  1%0K&  -  3  (6.18a) 

or 

tV  *  9  ^  "  tMLB^  (6-18b) 

Either  of  the  matrix  equations 

[V  ■  Wb0xm0L]  (6-'9a) 

[derived  from  Eq.  (6.18a)],  or 

Dfc]  ■  [\0]["oe] 

[derived  from  Eq.  (6.18b)]  constitutes  a  set  of  nine  algebraic  equations  (not  all 
Independent)  which  contains  the  variables  8p,  <j>p,  $F,  eb,  <pb,  e&,  <frg,  $H  and  <j»j. 
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Scattering 


Scattered 


All  of  the  machinery  for  performing  the  coordinate  transformations  needed 


to  Implement  the  scattering  program  is  contained  in  the  matrices  [Mg0J  (or 

[MqB]),  [M^]  (or  [MoL]),  and  the  matrix  equations  (6.19a,b).  Fortunately, 

we  only  require  limited  knowledge  of  the  variables  contained  in  those 

equations.  First,  we  must  solve  (in  terms  of  known  variables)  for  the  variable 

4>j  contained  in  [M^Q]  or  [MqL],  since  the  other  variables  in  those  matrices, 

da  and  <f>0,  are  known.  This  is  the  purpose  of  the  development  to  follow. 
p  p 

The  Barber  scattering  program  is  constructed  in  such  a  way  that  the 
scattered  wave  is  evaluated  in  the  (xL,  -  z^)  plane.  There  are  two  plane-wave 
scattering  processes  to  be  considered  (see  Figure  6.5).  We  will  call  these 
processes  A  and  6. 

Process  G  is  direct  scattering  from  the  body  to  the  antenna.  The  wave 
propagation  vector,  for  this  process,  denoted  by  IsA  and  totally  independent 
of  8*,  is  along  the  vector  r^  (Figure  6.5a)  from  the  center  of  the  scatterer 
to  the  antenna. 

Process  G  is  the  scattering  from  the  body  to  a  point  G  on  the  ground 
surface  followed  by  ground  reflection  toward  the  antenna  in  accordance  with 
the  law  of  reflection.  The  wave  vector  for  this  process,  denoted  by  !sg  and 
again  independent  of  8*,  is  parallel  to  the  vector  r$G  shown  in  Figure  6.5a. 

For  the  computation  of  the  scattered  fields  using  the  Barber  scattering 
program,  both  rsA  and  r$G  must  lie  in  the  (xL  -  zL)  plane.  Hence,  these 
vectors  must  be  perpendicular  to  the  y^  coordinate  direction.  The  conditions 
to  be  met  are  [with  the  aid  of  the  second  equation  in  the  matrix  relationship 
(6.17b)] 

IsA  *  (xA  ‘  xs)C-cos  9S  cos  *s  sin(0.j)A  -  sin  4>g  cos(^)A] 
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♦  (yA  -  y$)C'cos  sin  sin($j)A  +  cos  ♦g  cos  (*j)a3 


+  (ZA  ■  zS^sin  e0  s1n(*A^  *  0  (6.21a) 


Us  ’  i.  a  (xG  *  xS)C’cos  es  C0S  *B  sin(*A  "  s1n  H  ^A3 


+  A  ‘  ys^-C0S  9g s1n  ♦b  s1n(*A  +  cos  608  (♦A^ 


+  (zG  -  z$)[s1n  eg  s1n(^j)Q3  «  0 


(6.21b) 


where  ($j)^  and  Uj)q  are  values  of  $j  for  processes  A  and  G,  respectively, 
and  where,  from  Eqs.  (f  20e,g,f) 

y  -  XA  ZS  *  "S  2A 

s  zA  *  zs 

„  *A  zs  *  H  h 

H  h  *  H 


H 


0 


^A  *  X<XA  *  XS}  +  +  ^(2A  -  2S> 


Isg  “  x(xg  *  xs}  +  to G  '  yS}  *  2A  *  2$> 


Solution  of  Eq.  (6.21a)  for  cos(<f>j)A  and  s1n($j)A  or  that  of  Eq.  (6.21b) 


[pSG  COS  0g  COS(4»g  -  4>GS>  +  ZS  sin  90]  +  tp 


pSA*J(xA-xS)  ♦^A*JfS,‘ 


,  _  t  -1  ,  yA  '  , 

♦sa  tan  {  x'.  -“x;  ) 


pSG*J(xG-xS)  +(*G'yS>‘ 


-1  ,  yG  "  yS 


4>CJS  *  tan  { 


XG  '  XS 


XA  ZS  +  XS  ZA 

G  “W 


,  s2a 

G  2a  ♦  !s 


At  this  point  we  will  invoke  three  of  the  equations  in  the  matrix  equation 
system  (6.19a),  as  follows: 


<  W>jk 3  DhAk 


(6.23) 


(j  »  3,  k  •  1) 


cos  <pj(- Ceos  6b  sin  8g  -  sin  eb  cos  9g  cos($b  -  $g)]) 

+  sin  $j(s1n  8b  sin($b  -  $g))  *  sin  ep  cos  <f>p  (6.23a)’ 


(j  -  3,  k  *  2) 


cos  $j(s1n  8b  sin(<pb  -  $g) 


+  sin  4»j(-[sin  9b  cos  0g  cos(4>b  -  $g)  -  cos  0&  sin  9g]) 


sin  6„  sin  a 
P  P 


(6.23b) 


(j  -  3,  k  *  3) 

sin  9b  sin  0g  cos($b  -  $g)  +  cos  6b  cos  0g  *  cos  0p  (6.23c)’ 

Note  that  the  left  hand  side  of  Eq.  (6.23c)'  consists  entirely  of  given 
(and  therefore  known)  variables.  This  equation  can  be  used  to  evaluate  0p, 
thus  from  Eq.  (6.23c)', 

8p  *  cos-1[sin  eb  sin  0g  cos(d»b  -  $  )  +  cos  efa  cos  eg]  (6.24) 

From  Eqs.  (6.23a)'  or  (6.23b)'  we  can  solve  for  ($p)A  and  (<bp)G  (the 
azimuthal  angle  of  the  body  axis  in  the  lab  frame  for  scattering  processes 
A  and  G,  respectively)  in  terms  of  known  variables,  as  follows  (if  0  f  0; 

P 

otherwise  (<j>p)A  and  ($p)G  are  arbitrary): 

(♦p>a  =  cos'1  ( sTirr {_cos  (*A  Ccos  9b sin  9b 

p 

-  sin  efa  cos  0g  cos  (<J>b  -  4>g)]  +  sin  (4>j)A(sin  0b  sin  (<f>b  -  4>g))>] 

*  lequ,'",  s,"‘’ 1 A  kos  (s,n  9>  s,n  (*>  ‘  *•” 

+  sin  (<p j )A  [cos  0b  sin  0g  -  sin  0b  cos  0g  cos  (<f>b  -  #g)]}) 

(6.25a) 

where  cos  (<j>j)a  and  sin  Uj)A  are  given  by  Eqs.  (6.22a)  and  (6.22b),  respectively. 


6-20 


and  sin  0p  Is  obtained  from  Eq.  (6.23c)'. 

(♦p>G  "  C0S_1  (  THTe^  {'cos  (^}G  Ccos  9b  s1n  9& 

-  sin  eb  cos  eg  cos  (<frb  -  4>g)] 


+  sin  (0j)G  (sin  @b  sin  (d>b  -  d>6) }) 


*  (equivalently)  sin'1  (  jf-j-  { cos  (<j>j)G  (sin  6b  sin  (d>b  -  *g)) 

p 

+  sin  Uj)G  [cos  0b  sin  0g  -  sin  eb  cos  0g  cos  Ub  -  <j>g)]}) 

(6.25b) 

where  cos  (4>j)g  and  sin  (<l>j)g  are  given  by  Eqs.  (6.22c)  and  (6.22d),  respectively, 
and  sin  0p  Is  obtainable  from  Eq.  (6.23c)'. 

We  must  now  obtain  a  transformation  between  the  ground- frame  coordinates 
of  the  electric  field  plane-wave  spectrum  incident  on  the  scatterer  (l.e.,  the 
field  spectrum  directly  from  the  cable  superposed  on  the  ground-reflected 
field  spectrum  from  the  cable.  Both  of  these  are  expressed  in  ground-frame 
coordinates  (xQ,  yQ,  zQ)and  the  lab  frame  coordinates  of  these  same  field 
spectra  (i.e.,  the  x^  and  yL  components,  since  the  zL  component  must  be  zero 
for  a  plane-wave  propagating  in  the  zL  direction).  To  this  end,  we  write, 
where  E^  Is  the  Incident  field  (spectrum)  vector 


Si 


Eix. 


+  &  Eiy, 


+  Si<0) 


x  E 


lx 


+  l  E1y  +  *  Eiz 


(6.26) 


where 


iiL1 


E.  *  x.  component  of  E. 

1 L  -►! 

jL  *  y,  component  of  E. 
1yL  L  -V 


il\°h 


cix 

Siy 

L1* 


E^x  =  x  component  of 

E.  *  y  component  of  E. 
iy 

E.  *  z  component  of  E. 


and  where  [u,  ]  and  [u  1  are  defined  in  Eqs.  (6.7a)'  and  (6.17d)'  respectively. 

-+-L  -+0 

From  Eqs.  (6.26)  and  (6.18a,b), 


[|(,0)]Tt«oL][jL]  •  cl<1L,:TcuLi  (6.27) 

which  implies  [using  the  fact  that  [MqL]T  =  [MLq]  evident  from  the  definitions 
below  Eq.  (6.17b)]  after  taking  the  transpose  of  both  sides  of  Eq.  (6.26),  that 

[E<L)]  *  [MoL3T[|S0)3  -  [MLo][|S°)]  (6.28) 

in  longhand  notation  Eq.  (6.28)  takes  the  form 
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1 


\  Ccos  «6  CM  »J  -  cos  ®B  s1"  *8  s1"  ♦P 


+  Eiy{coTT  ^sin  *8  Cos  +  C0S  96  C0S  *8  s1n  (6.30a) 

$ 


if  efl  t 


8r1' 


StyL  *  gfx  {-  m h;  Ccos  *8  sin  *J  +  c“  e8  s1n  *6  cos  *0]) 


+  E.  {-  [sin  <t>$  sin  <j>j  -  cos  8g  cos  <f>B  cos  <t>j]}  (6.30b) 


if  0Q  / 


8  r  1  ‘ 


If  9g  3  then  Eq.  (6.28c)  implies  that  (where  subscripts  R  and  I  indicate  real 
and  imaginary  parts  respectively) 


<fiA 

^1y*R 


(flxil 

(Eiy>I 


■tan 


(6.31) 


From  Eqs.  (6.28a,b)'  and  (6.31),  If  6g  3  ^  ,  we  can  infer  the  following 

relationships  for  real  and  imaginary  parts  of  E.  and  E.  : 

1XL  1yL 


<£ix>« 

L  I 


-ixR 

I 


C- 


sin  <pj 

Sin  <) )r 


Eur  cos 


(6.32a) 


if  <J>g  t  0,  ir,  6g  »  £  • 
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R. 


A  further  step  is  required  to  construct  the  inputs  to  the  Barber  scattering 
program. 

The  amplitude  and  the  angle  4^. ,  the  polarization  angle  of  the  incident 
field  in  the  lab  frame,  are  required  as  inputs  to  the  scattering  program.  The 
amplitude,  independent  of  polarization,  is  easily  obtained  without  any  coordinate 
transformations,  i.e.  [with  the  aid  of  Eqs.  (6.29  and  (6.32)], 


li  ■  /Ef  '  'E?  *  /|E,XIS  *  |E(y!Z  + 


J{|Eix|*(!  +  tan2  eg  cos2  (f>g)  +  1 E .y 1 2(1  +  tan2  0g  sin2  <j>g) 


,.  +  2Re(E1x  Eiy  cos  $g  sin  $g  tan'1  0g)}  (if  eg  f  |  ) 
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'sin  <j>c 


If  6g  »  £  ,  f  O.ff 


♦  £*!!♦  ie. 


1f  \  •  I  •  *6  *  5  •  f 


(6.33) 


The  polarization  angle  assuming  that  the  phase  angles  of  E^x  and 


E.  are  the  same  (linear  polarization)  is  given  by 
iyL 

±  |E*  I 

<p.  *  tan"1  (  ) 


(6.34) 


where  E.  and  E.  are  given  by  Eqs.  (6.28a,b,c)‘  or  alternatively  by  Eqs. 

1XL  1yL 

(6.30a,b)  if  eg  f  ir/2  or  by  Eqs.  (6.32a,b)  if  8g  *  tt/2,  f  0,tt  or  by  Eqs. 
(6.32c,d)  if  eg  *  tt/2,  $g  t  ir/2,  3ir/2. 

If  the  phase  angles  of  and  are  not  the  same,  that  fact  would 
imply  that  the  wave  incident  on  the  scatterer  is  elliptically  polarized.  We 
could  exclude  that  possibility  by  assuming  that,  if  a  particular  source  should 
generate  a  superposition  of  elliptically  polarized  plane  waves,  then  for  each 
such  wave  the  rate  of  rotation  of  the  plane  of  polarization  is  sufficiently 
slow  that,  throughout  the  scatterer  region,  where  the  incident  wave  can  have 
an  influence  on  the  scattered  wave,  the  rotation  of  the  plane  of  polarization 
can  be  considered  approximately  constant. 

Once  the  fields  from  the  source  have  been  found,  we  can  determine  to 
what  extent  the  relative  phase  angles  between  x,  y  and  z  components  constitute 
a  problem.  Since  free  space  is  a  reciprocal  medium,  it  would  seem  that  it 
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k 


should  not  be  a  problem  at  all.  Any  plane  waves  propagating  from  the  slots 
should  be  linearly  polarized. 

To  ascertain  that  this  Is  indeed  the  case,  we  can  examine  the  expressions 
for  x,  y  and  z  components  of  fields  In  Eqs.  (4.46a,b),  (4.47c,d,e,f)  or  the 
specialized  forms  of  these  equations  (4.48b,c,d,e)  or  (4.49a,b,c,d,e,f) . 

To  obtain  the  field  components  we  must  add  the  quantities  with  subscript  + 
to  those  with  subscript  -  [i.e.,  Eqs.  (4.47c)  +  (4.47e),  (4.47d)  +  (4.47f), 
(4.48b)  +  (4.48c),  (4.48d)  +  (4.48e)].  The  results  of  the  first  set  of  these 
summations  are  as  follows  [where  we  note  that  =  (S+n^)*]: 


Eqs.  (4.47c)  and  (4.47e): 


X 

y 

±|8z|Re[sjn)  eJ*10] 

±|e2)Im[s|n)  e^10] 

i|8z!l«i[s<n) 

+  |8z(Re[s|n)  eJ4>1°] 

«AX  S<"> 

-Re[Aysjn)  e^10] 

-Re[Ax  s|n)  e^10] 

Im[Ay  S{+n)  e^10] 

(All  real  numbers) 

(All  real  numbers) 

±|6z|Re[s|n)  e^10] 

±|Bz|Im[sjn>  e^10] 

±|6z|lm[sin)  e^10] 

±|Bz|Re[s|n)  e^10] 

-Im[Ax  $jn)  eJ>1°] 

Re [Ay  s[n)  eJ>1°] 

Re[Ax  S<">  e^10] 

-Im[Ay  s[n)  e^10] 

(All  real  numbers) 

(All  real  numbers) 

From  the  chart  above,  it  Is  clear  that  all  three  components,  x,  y,  and  z. 


have  the  same  phase.  Each  component  associated  with  one  of  the  slot  fields, 

e.g.,  Hi?)  or  ,  has  a  complex  factor  common  to  all  three  components  and 
z 

another  complex  factor  that  is  different  for  the  three  components.  It  is  the 
latter  set  of  factors  that  the  chart  refers  to,  and  when  the  plus  and  minus 
terms  are  added,  the  sums  of  those  factors  are  all  real  numbers.  Hence,  the 
plane-wave  fields  from  the  slots  are  all  plane-polarized.  The  ground-reflection 
will  not  change  this,  since  the  ground  is  also  a  reciprocal  medium. 

The  same  arguments  apply  to  the  summations  [(4.48b)  +  (4.48c)], 

[(4.48d)  +  (4.48e)],  [(4.49c)  +  (4.49d)]  and  [(4.49e)  +  (4.49f)],  since  these 
are  only  special  cases  of  those  summations  discussed  above. 

Note  that  the  corresponding  factors  in  the  field  components  with  subscript  o 
[as  evidenced  by  Eqs.  (4.47a,b),  (4.48a)  and  (4.49a,b)]  are  all  real  numbers. 

Me  will  now  summarize  the  steps  required  to  implement  the  calculations 
described  in  this  section  on  the  computer. 

A.  Feed  in  the  variables  xA,  yA,  zA,  x$,  y$,  z$,  ®&,  <j>g,  6b,  <t»b  and 

(Eix,  E^,  Ejz),  the  components  of  the  spectral  field  of  the  plane  wave  incident 
on  the  scatterer,  from  the  main  program  to  a  subprogram  to  be  called  Subroutine 
BIS  ("Barber  JjnPut  Subroutine") . 

B.  Compute  Xg  and  yG  from  Eqs.  (6.20e,f)  (repeated  and  renumbered  for 
conveni ence): 


XA  ZS  +  XS  ZA 
XG  iA+*s 


(6.35a) 


„  yA  ZS  *  ^S  ZA 
*6  ’  ZA  +  ZS 


(6.35b) 
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C.  Using  Xq  and  from  Step  B,  compute  the  following  quantities  [see 
the  definitions  below  Eq.  (6.22d)]: 


PSA  *J<X A  -  XS}‘  +  <*A  * 


♦s*  ■ *•"  <  rrr  > 


cos  Uj)6 


_ PS6  cos  e6  cos  “  *6S>  *  ZS  s1n  e8 _ 

^[pSG  cos  9g  cos  U  -  $05)  +  zs  Sin  0g]2  +  [p$G  sin  -  *G$)]2 

(6.37c) 

~PSG  sin  ^6  ~  ^GS ^ _ 

S  "  *GSJ  +  ZS  sin  e6]Z  +  [pSG  sin  "  ^GS^2 

(6.37d) 

where 


sin 


ft* 


COS 


0g  cos 


E. 


PSA  -  J<XA  -  XS)Z  +  <*A  -  *S)5 


$Sa  *  tan"  ( 


-1 ,  yA  -  >S 


XA  -  XS 


) 


From  the  results  of  A,  B,  C,  and  D,  compute 


(4>  )A  from  Eqs.  (6.25a,b) 
p  G 


repeated  and  renumbered): 

(0p)A  *  cos  1  (  SJ  {-cos  (4»j)aCcos  eb  sin  0g  -  sin  eb  cos  0g  cos  ($b  -  $g)] 
+  sin  (4>j)A(sin  @b  sin  (*b  -  4>g))» 

»  (equivalently)  sin'1  (  {cos  (*j)A(sin  9b  sin  (*b  -  *)) 


+  sin  ($j)A[cos  0b  sin  0g  -  sin  0fa  cos  0g  cos  (*b  -  <j>g)]})  (6.38a) 
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'^tgafs#*****-*-  v--  --- 

JL i 

_ : _ _ _ I _ 


where  cos  ($j)A  and  sin  (4>j)a  are  given  by  Eqs.  (6.22a)  and  (6.22b)  respectively, 
and  sin  0p  Is  obtainable  from  Eq.  (6.23c)'. 

(<0p)s  *  cos'1  (  (-cos  Uj)g[cos  eb  sin  8g  -  sin  9b  cos  0g  cos  ($b  -  *g)] 


+  sin  (^j)G  (sin  8fa  sin  ($b  -  4>g))» 

*  (equivalently)  sin  ’’  Itttc  {cos  (4>j)G(s1n  0b  sin  (d>b  -  <j>g)) 

P 

+  sin  (4>j)q[cos  9b  sin  0g  -  sin  0fa  cos  0g  cos  ($b  -  <*>g)]})  (6.38b) 

where  cos  (<t>j)G  and  sin  (<pj)G  are  given  by  Eqs.  (6.22c)  and  (6.22d)  respectively, 
and  sin  0p  Is  obtainable  from  Eq.  (6.23c)'. 

F.  Compute  [from  Eq.  (6.33)  ,  repeated  and  renumbered]  the  amplitude  of 
the  plane-wave  incident  on  the  scatterer: 


G.  Compute  [from  Eq.  (6.34,  repeated  and  renumbered]  the  phase  of  the 
Incident  wave  for  both  scattering  process  (antenna-directed  and  ground- 
directed): 
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i 


(6.40) 


♦<*>■<«  -  [an’’  (  >]<A>-<6> 

|£<xj 

where  the  arguments  are  chosen  in  accordance  with  the  "lab  system" 
for  each  scattering  process. 

H.  Transfer  the  results  0p,  ($p)A,  ($p)s,  (<^)/«  and 

the  Barber  scattering  program. 


def i ned 


(4>1  )G  into 


6-32 


7.  THE  BARBER  SCATTERING  PROGRAM 


It  was  decided  early  in  this  project  to  use,  if  possible,  a  scattering 
program  that  was  already  developed  as  opposed  to  generating  a  new  one.  This 
seemed  sensible  in  view  of  the  time  limitations  of  this  project  and  the  ex¬ 
tensive  research  that  has  been  done  on  electromagnetic  scattering  during  the 
past  30  years.  It  might  have  required  virtually  all  the  project  time  and 
effort  available  to  develop  a  suitable  computer  program  to  treat  electro¬ 
magnetic  scattering  from  an  object  designed  to  simulate  a  human  frame  target. 

If  such  a  program  had  already  been  developed  and  was  available,  it  would  seem 
that  we  should  use  it. 

Hence,  early  in  the  project  an  extensive  literature  search  was  done  on 
electromagnetic  scattering  at  radio  frequencies.  The  bibliography  resulting 
from  that  search  is  included  in  the  list  of  references  at  the  end  of  this 
report?-1  through  0-35 

In  the  process  of  conducting  the  literature  search,  it  was  found  that 
virtually  all  the  scattering  theory  done  by  previous  workers  had  severe  limi¬ 
tations  when  one  considers  scattering  from  a  live  human  body.  The  theory  can 
only  be  done  rigorously  for  uniform  or  layered  spheres,  uniform  or  layered 
infinite  circular  cylinders  (not  finite  cylinders),  and  with  considerably  more 
difficulty  for  uniform  or  layered  ellipsoids  or  uniform  or  layered  cylinders 
with  elliptical  cross-section.  When  scatterer  dimensions  are  very  small  com¬ 
pared  with  wavelength  (or  more  precisely,  when  <<  1,  where  d  is  the  largest 
scatterer  dimension  and  A  is  the  wavelength  In  the  scatterer  material)  whatever 
series  of  special  functions  are  to  be  used  (e.g.,  spherical  or  cylindrical 

A-1 

vector  wave  functions  in  the  cases  of  spheres  and  cylinders,  respectively) 
require  fewer  terms  for  a  given  level  of  accuracy.  In  Implementing  scattering 
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problems  on  a  computer,  even  for  these  simple  tractable  shapes,  it  is  highly 
desirable  that  the  scatterers  be  (a)  of  uniform  constitutive  parameters,  as 
opposed  to  being  layered,  which  although  theoretically  tractable  always  in¬ 
volves  more  computer  time;  (b)  either  perfectly  conducting  (complex  dielectric 
constant  a  pure  imaginary)  or  perfectly  dielectric  (complex  dielectric  constant 
purely  real),  also  resulting  In  a  smaller  expenditure  of  computing  time;  (c) 
very  small  or  very  large  dimension  compared  with  wavelength,  allowing  the  use 
of  fewer  terms  in  the  appropriate  series  in  the  former  case  (as  remarked  above) 
or  the  use  of  physical  optics  approximations  in  the  latter  case.  Either  of 
these  two  extremes  simplifies  the  computations  and  thereby  saves  computer  time. 

Unfortunately,  none  of  the  simplifications  above  are  necessarily  applicable 

to  human  frame  targets  within  the  frequency  range  of  'interest.  Wavelength  in 

free  space  (XQ)  for  our  frequency  range  of  50  to  500MHz  ranges  from  0.6  to  6 

meters.  Considering  2  meters  as  d,  the  largest  dimension  of  a  human  frame 

target,  the  parameter  (  )  ranges  from  about  2  at  50MHz  to  about  20  at 

Ao 

500MHz.  If  we  consider  ,  where  X  is  the  wavelength  in  the  medium  of  which 
the  scatterer  is  composed,  we  must  correct  these  numbers  by  a  multiplicative 
factor  of 


ecs  being  the  permittivity  of  the  scatterer,  and  o$  is  the  conductivity  of 
the  scatterer.  Suppose  we  use  values  ec$  a  20  and  a$  =  1  mho/meter,  which 
were  close  to  the  values  used  in  our  computer  program,  based  on  Reference 
D-34. 

The  correction  factor  would  range  from  about  5  at  500MHz  to  about  20  at 
50MHz.  Thus  the  parameter  ,  referred  to  the  scatterer  material  rather 


than  free  space.  Is  somewhere  between  roughly  40  at  50MHz  and  roughly  100  at 
500MHz.  These  numbers  are  well  above  unity  but  not  necessarily  large  enough 
to  allow  the  use  of  physical  optics  approximations.  They  are  certainly  not 
small  enough  to  allow  the  use  of  highly  truncated  series  of  harmonic  functions 
and,  in  fact,  require  large  numbers  of  terms  for  convergence  of  these  series. 
Also,  human  targets  cannot  be  modelled  as  perfect  conductors  or  perfect  di¬ 
electrics  at  these  frequencies,  which  again  brings  the  problem  into  the  domain 
of  more  computer-intensive  scattering  problems.  A  sphere  or  circular  cylinder 
is  an  extremely  crude  model  of  a  human  body.  An  ellipsoid  would  be  a  much 
better  model,  but  again  it  would  be  desirable  from  a  computer-time  viewpoint 
to  have  some  symmetry  in  the  scatterer;  nence,  a  good  compromise  between 
analytical  simplicity  and  realism  is  a  spheroid,  particularly  a  prolate  spheroid. 

As  a  part  of  the  literature  search,  the  bioengineering  literature  on 
absorption  and  scattering  of  electromagnetic  waves  by  biological  objects  was 
perused?'1  throu^h  0-35  Bioengineers  are  primarily  interested  in  absorption 
rather  than  scattering,  but  the  boundary  value  problem  they  must  attack  to 
determine  the  absorption  also  contains  the  mathematical  machinery  to  determine 
the  scattering  from  the  body.  Hence,  a  particular  piece  of  work  done  by  Pro¬ 
fessor  Peter  Barber  of  the  Bioengineering  Department  at  the  University  of  Utah 

in  Salt  Lake  City  appeared  to  be  highly  applicable  to  the  problem  of  interest 
D-25 

on  this  project. 

Barber  had  authored  or  co-authored  a  number  of  papers  on  tnis  work  and 

had  developed  a  very  extensive  computer  program  to  implement  the  analytical 
sol utions?'1 ^ *^_1 ^ 9 ,D-22 ,0-24 ,0-29 ,D»30 ,D-32  ,D-34  The  prQgram  ha(J  been 

used  on  a  number  of  problems  of  Interest  to  bioengineers  over  a  period  of  about 
five  years  and  had  been  thoroughly  "debugged."  It  had  been  used  by  others  on 
different  computer  facilities  and  had  produced  reliable  results.  It  seemed 
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that  tha  best  procedure  for  us  would  be  to  ask  Barber  for  his  program  cards 
and  adapt  the  program  to  Northeastern's  VAX.  All  of  this  was  done,  as  Indicated 
In  the  quarterly  status  reports  on  this  project,  and  Professor  Barber's  scat¬ 
tering  program  was  adapted  to  the  VAX  and  implemented  as  a  subroutine  in  our 
overall  program.  It  is  called  "Subroutine  BARBER." 

The  details  of  the  theory  behind  this  scattering  program  are  given  in 
Barber's  papers,  particularly  References  D-13,  D-17,  D-19,  D-22,  D-29. 

It  is  based  on  an  integral  equation  technique  called  the  "Extended  Boundary 
Condition  Method"  (EBCM) . 

0-4 

Following  a  development  due  to  Waterman,  the  incident  and  scattered 
electric  fields  at  position  r  =»  (r,  0,  $)  are  expanded  in  a  series  of  spherical 

-a* 

vector  wave  functions  as  follows  (using  some  of  Barber's  notation): 

E*(r)  «  l  0  [a  M](kr)  +  b  (kr)]  (7.1a) 

V3 1 

where  superscripts  i  and  s  refer  to  "incident"  and  "reflected",  respectively, 

where  M1 ,  ft1 ,  M3,  N3  are  vector  wave  functions,  superscripts  1  and  3  referring 

v  v  v  v 

to  particular  classifications  of  these  wave  functions  (e.g.,  3  refers  to  radi¬ 
ation  fields,  i.e.,  Hankel  function  expansions),  k  being  the  propagation 
constant  in  the  ambient  medium,  v  being  a  combined  index  incorporating  spheri¬ 
cal  harmonic  indices,  and  0^  being  a  normalization  constant.  The  constants  ay 
and  by  are  assumed  known  for  a  specified  incident  wave.  The  task  of  solving 
for  the  scattered  wave  field  is  that  of  evaluating  the  unknown  constants  f 

and  g  . 

3v 
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The  method  Is  valid  for  any  homogeneous,  linear,  isotropic  scatterer, 
and  has  also  been  extended  to  include  layered  scatterers  each  of  whose  layers 
is  homogeneous,  linear  and  isotropic.  The  technique  lends  itself  most  readily 
to  scattering  from  spherical  objects  but  can  be  used  to  treat  objects  of 
arbitrary  shape. 

The  scattering  object  Is  completely  surrounded  by  a  sphere  whose  diameter 
is  the  largest  dimension  of  the  object  itself.  If  the  fields  on  the  surface 
of  the  sphere  are  known,  then  the  Kirchhoff-Huyghens  integral  equation  [see 
Appendix  I,  Eqs.  (I.3a,b)  or  {I.6a,b)]  can  be  used  to  determine  the  fields 
outside  the  sphere.  To  relate  the  fields  on  that  sphere  (a  knowledge  of  which 
is  tantamount  to  solution  of  the  scattering  problem)  to  the  boundary  conditions 
on  the  bounding  surface  of  the  nonspherical  scattering  object  itself,  an 
analytic  continuation  process  is  used.  Through  the  Kirchhoff-Huyghens  integral 
equation,  the  expansions  in  Eqs.  (7.1a,b)  and  the  application  of  the  boundary 
conditions  on  the  scatterer's  surface,  a  set  of  matrix  equations  are  developed 
from  which  the  coefficients  fy  and  gv>  and  hence  the  scattered  fields,  may  be 
determi ned . 

Further  details  on  the  basic  theory  behind  the  EBCM  is  explained  quite 
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thoroughly  in  Waterman's  1971  paper,  in  which  he  also  references  a  great  deal 
of  previous  work  on  this  and  related  techniques  for  solving  scattering  problems. 
Barber's  principle  contribution  was  to  apply  that  theory  to  the  development  of 
computer  programs  that  can  handle  a  wide  range  of  difficult  scattering  problems 
and  the  application  of  the  theory,  via  those  programs,  to  problems  involving 
biological  scatterers.  The  particular  scatterer  to  which  we  have  applied  his 
program  is  a  homogeneous,  linear,  isotropic  prolate  spheroid  designed  to  simu¬ 
late  a  human  being. 

It  should  be  mentioned  that  another  analytical  method  of  treating  scattering 


i 


from  a  human  frame  target,  together  with  a  computer  program  to  Implement  the 
method,  was  developed  by  Barber  and  one  of  his  Ph.D.  students,  Mark  Hagmann. 
This  work  is  presented  in  Hagmann ‘s  Ph.D.  thesis.  Reference  D-31  on  our 
reference  list. 

Professor  Barber  gave  the  writer  a  copy  of  this  dissertation  for  possible 
use  on  the  project.  The  method  treats  a  human  body  as  an  electromagnetic  scat- 
terer  in  a  much  more  accurate  way  than  does  the  theory  that  we  actually  used. 
Head,  torso  and  limbs  are  each  modelled  and  scattering  processes  from  all  of 
these  parts  and  interactions  between  them  are  considered.  The  unfortunate 
aspect  of  this  work  from  our  point  of  view  is  that  it  is  an  almost  completely 
numerical  technique  and  is  much  more  computer- intensive  than  the  method  we  have 
actually  used.  An  attempt  to  use  this  technique  would  probably  have  exhausted 
our  computation  resources  before  we  could  have  obtained  any  significant  results 
For  this  reason,  we  decided  upon  the  simpler,  more  analytical  method,  wherein 
the  "human"  scatterer  is  modelled  in  a  somewhat  crude  manner,  but  the  amount 
of  computer  time  required  to  determine  a  point  is  large  but  still  reasonable. 


8.  SCATTERED  FIELDS-COORDINATE  TRANSFORMATION 
AT  OUTPUT  OF  SCATTERING  PROCESS 

The  scattered  field,  as  produced  by  the  Barber  scattering  program,  is 
measured  along  the  (xL  -  y^)  plane  in  the  lab  frame.  This  is  illustrated 
in  Figure  8.1.  The  scattering  angle  is  the  polar  angle  of  85,  the  wave 
propagation  vector  for  the  scattered  wave,  measured  in  the  lab  frame.  (Note 
again  that  85  is  entirely  unrelated  to  8*). 

The  Barber  scattering  program  yields  as  its  output  the  "vertically 
polarized"  and  "horizontally  polarized"  components  of  the  plane-wave  field 
scattered  in  the  e$  direction.  The  former,  denoted  by  Esv,  is  the  component 
normal  to  the  (xL  -  zL)  plane,  and  the  latter,  denoted  by  ESH,  is  the  compo¬ 
nent  parallel  to  the  (xL  -  yL)  plane.  From  E^v  and  ESH,  we  can  determine  the 
lab  frame  components  of  the  scattered  field.  We  must  then  find  the  ground 
frame  components  of  the  field.  The  procedure  for  accomplishing  this  will  be 
described  below. 

From  Eqs.  (6.26)  and  (6.18a,b)  applied  to  in  lieu  of  E^ 

-  cgL)]Tc\0][ut)]  (8.1) 

where 


11 

ESX 

,  CE<1)]. 

V 

[e'o)]  . 

ESy 

5* 

Esz 
■  • 

Sk 

Taking  the  transpose  of  both  sides  of  Eq.  (3.1)  and  noting  again  that 
[MoL3T  3  [MLq],  we  obtain 
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(8.2) 


Ce|0)]  *  CmoL][E<l)] 

In  longhand  [v.om  the  matrix  definitions  (6.17a ,b) ' ] 

E$x  =  (cos  es  cos  cos  ♦j  -  sin  6g  sin  ♦j)  Ej 


(cos  6g  cos  <j)g  sin  +  sin  (jig  cos  4>j ) 


+  (sin  0g  cos  $g)  E 
(8.2a)' 


E$y  *  (cos  0g  sin  (tg  cos  <j>j  +  cos  j>g  sin  <j»j) 


-  (cos  0g  sin  <t>g  sin  d>j  -  cos  j)g  cos  <j>j) 


ESyL  +  <s1"  9b 


sin  <j»g)  E 
(8.2b)1 


hz  =  -{sin  6S  cos  ♦j)  ESxl  +  (sin  0B  sin  ESyL  +  (cos  V  ESzl 

(8.2c) ‘ 

The  angle  0SL  (see  Figure  8.1)  is  the  polar  angle  of  the  vector  6s  as 
measured  in  the  lab  frame.  It  is  evident  from  Figure  (8.1)  that 


and  hence,  from  Eqs.  (8.2a,b,c)‘  and  (8.3a,b,c) 


[(cos  0g  cos  $g  cos  (<j>j)A-s1n  $g  sin  ($j)A)  cos  ®sj_ 

5  6 


w  JA 

(sin  0g  cos  $g)  sin  8^  ] 


•[cos  6g  cos  $g  sin(d>j)A  +  sin  $g  cos  (^j)a^ESv 

S  G 


(8.4a) 


[(cos  0g  sin  <t>g  cos  Uj)A  +  cos  $g  sin  (<t>j)A)  cos  0^J 

G  G 

(A)  JA) 

(sin  8g  sin  $g)  sin  0^1  E^ 


[cos  0g  sin  4>g  sin  (<|>j)A  -  cos  <|>g  cos  ($j)A]  E^ 

G  G 


(8.4b) 


[-(sin  6g  cos  (<>j)A)  cos  9^' 
G 


(A)  (A) 

(cos  0g)  sin  6^] 


(A) 

+  [sin  0g  sin  ($j)A] 

G 


(8.4c) 


where  subscripts  and  superscripts  (A)  and  (G)  in  Eq.  (8.4a,b,c)  correspond  to 
A  and  G  scattering  processes,  respectively,  and  where  cos  Uj)A»  sin  (4>j)a, 
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cos  ($j)q  and  sin  ($j)q  are  obtained  from  Eqs.  (6.22a,b,c  and  d)  respectively. 

Confining  attention  In  what  follows  to  the  A  scattering  process,  it  remains 
to  evaluate  cos  0^  and  sin  0^.  Since  8^  (parallel  to  r^)  is  in  the 
(xL  -  z^)  plane,  the  lab-frame  azimuthal  angle  <j>SL  «  0  or  i.  It  follows  from 
the  transformation  between  spherical  and  rectangular  coordinates  that  (since 
the  antenna  lives  In  the  xL  -  z^)  plane) 


cos  0 


sin  0, 


,(A) 

.3la 

SI 

rSA 

,( A ) 

.  ^xlaI 

SL 

rSA 

(8.5a) 


(8.5b) 


where 


rSA  *  -  XS)2  +  <yA  -  ys>2  +  (ZA  *  ZS^2 

and  where  xLA,  y^A,  zLA  are  the  lab- frame  coordinates  of  the  antenna.  (Noting 
that  the  lab  frame  has  its  origin  at  the  scatterer  center,  it  follows  that 
XLS  *  yLS  =  ZIS  =  °*  A1so*  since  the  antenna  Is  in  the  (xL  -  zL)  plane,  we 
know  that  yLA  =  0. ) 

From  Eq.  (6.16a)  applied  at  the  antenna  position 

XLA  =  [cos  0g  cos  $g  cos  ($j)A  -  sin  sin  Uj)A](xA  -  x$) 

+  [cos  6g  sin  $g  cos  (*j)A  +  cos  sin  (<frj)A](yA  -  ys) 

-  [sin  0g  cos  (4>j)a3(za  -  z$)  (8.6a) 
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yLA  "  "^C0S  06  C0S  *6  Sln  +  Sin  C0S  ^J^(XA  -  XS^ 

-  [cos  eg  sin  (j»g  sin  ($j)A  -  cos  $0  cos  Uj)A](yA  -  y$) 

+  [sin  8g  sin  (4>j)A](zA  -  z$)  ■  0  (8.6b) 

zLA  *  (sin  eg  cos  <1>e)(xA  -  x$)  +  (sin  eg  sin  <*>g)(yA  -  y$) 

+  (cos  0g)(zz  -  z$)  (8.6c) 


We  will  now  summarize  the  procedures  required  to  transform  the  output  of 

the  Barber  scattering  program  into  the  basic  (x,  y,  z)  coordinate  system. 

~(A1  ~(A1 

(A)  The  complex  field  components  and  Egv'  are  extracted  from  the 
output  of  the  Barber  program.  These  are  respectively  the  "horizontally"  and 
"vertically"  polarized  scattered  fields  in  the  direction  of  the  antenna  in 
response  to  the  incident  plane-wave  with  propagation  vector  k^.  These  field 
components  are  fed  from  the  Barber  program  into  a  subprogram  which  we  call 
Subroutine  BOS  ("jlarber  Output  Subprogram") .  The  variables  xA,  yA,  zA>  x^, 
y$»  z$,  4>e  and  eg  are  read  into  the  BOS  from  the  main  program 

(B)  The  following  computations  are  made  [introduced  below  Eq.  (6.22d)]: 

^A  *  /<* A  -  XS)2  +  <yA  ■  yS)Z  {8-7a) 


0SA  -  tan 


yA  ~  yS 
XA  '  x 


(8.7b) 
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i 


l 


and  s' 
for  c 


[psA  cos  0g  cos  (d>0  -  -  (zft  -  z$)  sin  0g]4  +  [p 


mmm 


yLA  “  '^C0S  ®S  cos  *g  s1n  +  Sin  C0S  ^A^XA  "  x$) 


-  [cos  e6  sin  ♦g  sin  (<frj)A  -  cos  $g  cos  ($j)A](yA  -  y$) 


+  [sin  0g  sin  (<i>j)A3< ZA  -  z$)  *  0 
zLA  -  (sin  8g  cos  4>g)(xA  -  x$)  +  (sin  0g  sin  <j>g)(yA  -  y$) 


+  (cos  eg)(zA  -  z$) 


(8.9b) 


(8.9c) 


(E)  From  the  results  of  (A),  (B),  (C),  and  (D),  cos  0^  and  sin  0^ 
are  computed  from  Eqs.  (8.5a,b)  (repeated  and  renumbered  here): 


cos  0  W  =  iA 
bL  rSA 


sin  =  l^LAi 

s "  9sl 


(8.10a) 


(8.10b) 


(F)  From  the  results  of  (A),  (B),  (C),  (D)  and  (E),  compute  E^,  E^ 


and  Eg\);  from  Eqs.  (8.4a,b,c)  (repeated  and  renumbered  here): 


*  [(cos  @g  cos  0g  cos  (d>j)A  -  sin  $g  sin  Uj)A)  cos  ej 


-  (sin  9g  cos  $g)  sin  e^]  E^ 


-  [cos  0g  cos  0g  sin  (<pj)A  +  sin  0g  cos  (4>j)a3  Egy  (8.11a) 


-  *) 


*«u*- 


l 


ESy  *  ^cos  08  sin  ♦?  cos  +  cos  *3  sin  V  cos  eSL 

-  (sin  0g  sin  <j»&)  sin  ©jj^] 

-  [cos  0g  sin  $g  sin  (<*iJ)A  -  cos  $g  cos  (<J>j)A]  EjJ*  (8 

^Sz  *  £"(sin  eg  cos  <4>j)a)  cos  0SL^  ‘  ^cos  eg^  sin  4i]l  ^ 

G 

+  [sin  0g  sin  (<pj)AJ  E^  (8 


(G)  Return  the  results  of  (F)  to  the  main  program. 


9.  EFFECT  OF  GROUND  REFLECTIONS 
ON  THE  SCATTERED  FIELD 


In  this  section,  we  will  evaluate  the  effect  of  ground  reflections  on 
the  scattered  field  at  the  antenna.  We  will  also  discuss  the  procedures  for 
implementing  these  calculations  on  the  computer. 

The  calculations  to  be  performed  to  obtain  this  term  in  the  total  field 
can  be  divided  into  two  phases,  as  follows: 

Phase  1:  Calculation  of  the  field  of  the  plane-wave  scattered  from  the 
body  toward  the  ground  reflection  point  G. 

Phase  2:  Calculation  of  the  plane-wave  reflected  from  G  and  propagating 
toward  the  antenna  A. 

Phase  1  of  these  procedures  begins  with  Eqs.  (8.4a,b,c)  with  superscripts 
and  subscripts  G  rather  than  A.  This  step  is  followed  by  the  equivalent  of 
Eqs.  (8.5a,b)  with  G  substituted  for  A,  i.e.. 


cos  9 


(G)  = 

SL  rSG 


(9.1a) 


sin  9 


(G)  = 

SL 


LG1 


$G 


We  note  that,  in  this  case 


(9.1b) 


rSG  *  '  XS)2  +  •  yS)2  +  (zG  •  zS)Z 


(9.2) 


and  [from  Eqs.  (6.20e,f,g)] 


G  ■  Ground  reflection  point 
A  ■  Antenna 

Figure  9-1.  Ground-reflected  Scattered  Field 


It  will  be  recognized  that  rSG  as  given  by  Eq.  (9.3)  is,  except  for  the 
factor  2s/(zs  +  zA),  the  distance  from  the  "image11  of  the  scatterer  center,  lo¬ 
cated  at  (x$,  y$,  -zs),  to  the  antenna.  The  factor  zs/(zs  +  zA)  implies,  as 
supported  by  our  intuition,  that  if  z$  =  zA,  then  rSG  is  one-half  of  the 
image-to-antenna  separation  distance.  If  zs  »  zA,  then  r^G  the  image-to- 
antenna  separation  distance,  and  if  z^  <<  zA,  then  r$G  is  the  separation 
distance  multiplied  by  the  ratio  of  scatterer  center  height  to  antenna  height. 

The  lab  frame  coordinates  of  the  ground-reflection  point  G,  (xLG,  y^G,  z^q)i 
are  calculated  from  equations  analogous  to  Eqs.  (8.6a,b,c),  where  ( d>j )q  replaces 
and  coordinates  (xG,  yQ,  zQ)  replace  (xA,  yA,  zA). 

We  will  now  summarize  the  steps  required  to  implement  these  calculations 
on  the  computer.  The  steps  A  through  G  below  are  analogous  to  signs  A  through 
G  in  Section  8. 


A.  The  complex  field  components  E^  and  E^  are  extracted  from  the 


output  of  the  Barber  program.  These  components  are  analogous  to  E 
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r(A) 


SH 


l 


A1.. 


~(h) 

and  £jV  ,  where  the  scattering  process  G  (scatterer  to  ground  reflection 
point)  replaces  the  process  A  (scatterer  to  antenna).  These  components 
are  fed  from  the  Barber  program  Into  subroutine  BOS  (see  Step  A  in  the 
summary  in  Section  8).  It  was  indicated  in  Section  8  that  vari¬ 
ables  x^,  y^,  zA,  Xj,  ys,  Zj,  dg  and  0g  are  read  into  BOS  form  the  main 
program.  Also  v,  the  complex  refractive  index  of  the  ground,  is  read 
into  BOS  from  the  main  program. 

8.  The  following  computations  are  made:  [analogous  to  Eqs.  (8.7a,b,c) 

where  xs,  yQ,  zQ  replace  xft,  yft,  z^],  [with  the  aid  of  Eqs.  (9.2a,b,c)'] 

0s;(  (9.4a) 


PSG 


^  >  ’  »SA 


(9.4b) 


[Eq.  (9.3)  repeated  and  renumbered] 


SG 


(  -^TiJ  )  -  *s)2  ♦  (yA  *  *s>2  +  (ZA*  ZS)2  (9-4c) 


C.  Using  the  results  of  B,  cos  Uj)q  and  sin  Uj)g  are  computed  from 
Eqs.  (6.22c,d)  (repeated  and  renumbered  below):  1 


cos  Uj)Q 


PSG  C0S  6S  C0S  ^6  "  ^GS^  +  H  sin  ®3 


/[pSG  cos  93  cos  Us  -  +  sin  egD2  +  Cpsg  sin  Ug  -  <tGS)]2 

(9.5a) 
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sin  Uj)G 


_ _ ~PSG  sin  ~  _ - 

/[PSG  C0S  6S  cos  "  *GS}  +  ZS  sin  eB]2+[pSG  s1n  (*B  "  ♦gS^ 

(9.5b) 

where 

PSG  ’  J(xG  -  XS)2  +  (yG  -  y$>2 


4>2q  =  tan 


-1 


yG-ys 

XG  -  xs 


)  a  tan'^ 


,  yA  -  ys 

*A  '  *S 


) 


D.  Using  the  results  of  B  and  C,  xLG,  yLG  and  z^G  are  computed  from 
equations  equivalent  to  Eqs.  (8.6a,b,c)  where  xG>  yQ,  zG  and  Uj)G 
replace  x^,  yA,  z^  and  (<j>j)^,  respectively  (As  is  the  case  with  the 
antenna,  the  ground-reflection  point  G  is  in  the  (x^  -  zL)  plane; 
hence,  we  know  that  y^G  3  0  and  we  don’t  need  Eq.  (9.6b)  in  the  com¬ 
putations). 


xLG  *  [cos  0S  cos  cos  -  sin  <j>g  sin  (<J>j)g3( xQ  -  x$) 

+  [cos  9g  sin  <*>g  cos  (4>j )Q  +  cos  4>g  sin  Uj)G](yG  -  y$) 

+  [sin  9g  cos  (<$j)g]  z$  (9.6a) 

yLG  =  '^C0S  96  C0S  s1n  +  Sin  ^6  C0S  ^G^G  '  XS^ 

-  [cos  9g  sin  <j)g  sin  <4>j )Q  -  cos  <t>6  cos  (<t>j)G3(yG  -  ys) 
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-  [sin  eg  sin  Uj)G]  Z$  »  0 


(9.6b) 


zLG  =  (sin  6g  cos  4>g)(xG  -  xs)  +  (sin  9g  sin  $g)(yG  -  y$) 


-  (cos  eg)  zs 


(9.6c) 


E.  From  the  results  of  A,  B,  C,  and  0,  cos  0gj^  and  6^  are  computed 
from  Eqs.  (9.1a,b)  (repeated  and  renumbered  below): 


cos  eg)  - 
bL  r$G 


sin  e(6)  ,  rtG 

S'n  ^  ’•sg 


(9.7a) 


(9.7b) 


F.  From  the  results  of  A,  B,  C,  D,  and  E,  compute  E^,  E^,  and  E^ 
from  Eqs.  (8.4a,b,c)  (repeated  and  renumbered  below): 


Eg^  =  ^cos  0g  cos  cos  ^J^G'sin  ^8  sin  cos  0SL^ 


-  (sin  eg  cos  <frs)  sin  0^3  Eg{^ 

-  [cos  0g  cos  0g  sin  ( j ) q  +  s i n  <f>g  cos  (<J»j)Q]  E^  (9.8a) 
Eg^  ■  [(cos  0g  sin  <frg  cos  )G  +  cos  <pg  sin  (d»j)s)  cos  0^ 

-  (sin  8g  sin  <pg)  sin  0^]  E^ 


i 


(9.8b) 


-  [cos  6g  sin  $g  sin  (<j>j)G-cos  $g  cos  (<bj )g3  E^ 

*  C-(sin  0g  cos  (4>j)g)  cos  6^-  (cos  6g)  sin  0^] 

+  [sin  9g  sin  (<t>j)g3  E^  (9.8c) 


G.  Return  the  results  of  f  to  the  main  program,  in  which  the  remaining 

part  of  the  computation.  Phase  2,  is  performed. 

***** 

To  continue  with  Phase  2,  we  must  calculate  the  fields  of  the  ground- 
reflected  plane-wave  due  to  the  scattering  process.  To  this  end,  we  invoke 
equations  analogous  to  Eqs.  (5.18a)  and  (5.19),  but  in  a  different  context. 

Equations  (5.18a)  and  (5.19)  contain  the  (x,  y,  z)  comoonents  of  the 
vectors  8+  and  3”.  The  analog  of  0”  in  the  present  development  is  a  vector 

-V  -+ 

85,  which  we  will  rename  8^,  and  which  is  the  unit  vector  parallel  to  a  vector 
originating  at  the  scatterer  center  and  terminating  at  the  ground- reflection 
point  G.  The  analog  of  £+,  which  will  be  called  j^,  is  the  unit  vector  parallel 
to  a  vector  originating  at  G  and  terminating  at  the  antenna  A.  It  is  easily 
deduced  from  these  definitions  that  (where  vectors  rSG  and  r^)  8^  and  8^  are 
illustrated  in  Figure  9.1. 


CsG  j(xG  -  *S>  *  -  V  '  l(zS> 

rsG  K  -  *s>2  *  <*a  -  *s>2  ♦  4 


(9.9a) 


and 


•**’ 


8Sh  and  ySz,  respectively. 

We  can  now  return  to  the  summary  of  steps  in  programming  of  the  calcu¬ 
lations.  The  steps  in  implementation  of  Phase  2,  continuing  after  Step  G,  are 
as  follows: 

H.  Compute  the  following  variables  which  will  be  needed  in  subsequent 
computations  [from  Eqs.  (9.9a,b)  and  (9.10a,b,c)] 


GA 


ZA  +  ZS 


XA  -  XS>2  +  <*A  •  >S>2  +  ^ZA  +  ZS)2 


(9.12a) 


(9.12b) 


2S 

ZS  +  ZA 


-  ys} 

rSG 


(9.12c) 


(9.1 2d) 


(9.12e) 


ZA  <XA  ‘  XS> 
ZS  +  ZA  rGA 


(9.1 2f ) 


ZA  ,  (*A  -  *S> 
ZS  *  ZA  rGA 


(9. >29) 


(@Sx  )2  +  (fijy  ) 


( 9 . 1 2h ) 


-1 


(9.121) 


(9.12J) 


(9.12k) 


I.  Using  the  computations  of  E^,  E^  ,  and  E^  performed  in  F  and 
the  variables  computed  in  H  as  inputs  to  this  step,  compute  [Ej^] 
through  the  matrix  equation  in  Eq.  (9.11)  (repeated  and  renumbered) 


[E^]  =  [R^GA)][E^G)] 

where 


Co^BSx’  6Sy) 


r(GA) 

kejk 


(9.13) 


Co^SSx’  BSy^ 


<BS«)2  -  <BSy)Z  -  141  y 


+ 

Sz 
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10.  TOTAL  SPECTRAL  FIELD  AT  THE  OBSERVATION  POINT 
ANO  ITS  INVERSE  FOURIER  TRANSFORM 


The  plane-wave  spectrum  of  the  x,  y,  and  2  components  of  the  total 
electric  field  at  the  observation  point  (receiving  antenna  position)  is  the 
end  result  of  the  calculations  described  In  the  previous  sections. 

Referring  back  to  Section  2,  the  plane-wave  spectrum  of  the  electric 
field  is  the  superposition  of  those  of  four  fields,  as  follows: 

(1)  The  contribution  from  the  cable  slots  as  if  the  cable  were  in  free 
space  (the  sum  of  the  fields  from  all  the  slots,  neglecting  the 
effect  of  the  ground).  The  individual  slot  fields  are  denoted  by 
E^(gh)  in  Section  2.  The  calculation  of  these  fields  is  described 
in  Section  4. 

(2)  The  contribution  resulting  from  ground  reflection  of  the  plane-wave 
spectrum  of  the  field  resulting  from  excitation  of  the  slots,  which 

is  denoted  by  E^*3^)  in  Section  2  and  whose  calculation  is  described 
in  Section  2. 

(3)  The  contribution  due  to  scattering  from  the  target  directly  toward 
the  antenna  of  the  superposition  of  direct  and  ground- reflected  plane- 
wave  fields  [i.e.,  E^(gh)  +  E^(g  )],  which  is  denoted  by  E^c^(g.) 

-►  -b-n  -*■  -*n  -*•  -h> 

in  Section  2  and  whose  calculation  is  described  in  Sections  6,  7, 
and  8. 

(4)  The  contribution  due  to  ground  reflection  of  the  scattered  field, 
which  Is  denoted  by  E^(gk)  in  Section  2  and  whose  calculation  is 
is  described  in  Section  9. 

The  final  step  In  the  analysis  Is  a  two-dimensional  Fourier  transformation 
which  transforms  the  spectral  field  components,  which  are  functions  of  8^,  into 
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actual  field  components  which  are  functions  of  position. 


In  what  follows  we  will  consider  this  problem  In  a  generic  sense.  Con¬ 
sider  a  spectral  field  component  denoted  by  l(|h).  We  would  like  to  perform 
a  double  Fourier  transform  to  convert  this  Into  a  field  component  E(r),  where 
r  *  (x,y,z)  3  arbitrary  point  In  space. 

The  transformation  equation  Is  (see  Appendix  III) 


E(r) 


e 


jk  8*r  .. 
-  * 


We  now  convert  both  r  and  |  to  polar  coordinates,  l.e.. 


(10.1) 


r  3  p( x  cos  <j>  +  y  sin  $)  +  zz 
■*  -*•  ■* 


(10.2) 


where 


/ 


x  +  y 


tan'1  (  £  ) 


Also: 


6  3  8h(x  cos  +  l  sin  ♦g)  +  z  8Z 


(10.3) 


where 


tan’ :  (  ^  );  8 


•±  le2l.  Iezl  • 


v 

f 

i 


We  express  £(8^)  as  a  function  of  $3  and  8^*  l.e.,  we  rename  the  function 
arguments  as  follows: 
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£(8h)  -  E(0h,  ♦g) 


00.4) 


From  Eqs.  (10.1),  (10.2),  (10.3),  and  (10.4),  we  can  write 


E(r)  -  (  -sr  ) 


j  k  ,  rl  ,2ir 

O 


d  Sh  6„  |  (I  *fl  E(Sh.  ♦„) 
0 


ejkop8hcos(0s'<(,)  .  ejkoz^'sh 


Jk«zJi-s! 


(10.5) 


which  reflects  the  fact  that  Bh  <  1. 

We  could  now  invoke  the  well-known  Bessel  function  relationship 


gjxsine  J  (x)  fijne 


n=-8 


00.6) 


which,  when  applied  to  Eq.  (10.5),  would  result  in  the  expression 


j  ko,2  *  *-j[n*-(ir/2)]  rl 


I.** 


n=-8 


d  6h  6h  Jn(ko  »  eh1 


2iT  jk  zll-6u  jnd>Q 

d  *8  E(8h>  ♦„)  e  °  h  e  B 


(10.7) 


Equation  (10.7)  is  perfectly  satisfactory  from  an  analytical  viewpoint, 
but  preliminary  investigation  shows  that  it  would  be  prohibitively  time- 
consuming  to  implement  on  the  computer.  As  an  alternative  technique,  we  write 
the  <j)g  integral  as: 
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v  !  J  i(|v 
u  0 


Kiy  ♦;  ♦  <')  * 


.1 ov;. 


(10.C) 


whore  =  <|)„  -  *i> 

0  p 

Recognizing  that  the  function 


jk  pB.  cos,>' 

P,  4->  flj  -=  E(Hh,  +  *)  e  B 


(10.9) 


can  be  expanded  in  an  exponential  Fourier  series  in  4>',  we  write 


F(4>A;  p.  <!>>  M  =  }’  c„(p.  Pu)  e 


(10.10) 


It  follows  from  Eqs.  (10.0),  (10.9),  arid  (10.10)  that 


00  r?.  v  jn*' 

I*  =  l  cn(fs  <!>;  f’h)  d  <p'  e  15  =  2v  c  (p,  3h)  (10.11) 

vi3  n=-ro  l 


where  we  note  that 


1  r?-i  i  ..  Jknp3hC0S<!,ft 

co(p’  -  21  d  ♦£  E<V  *3  +  e 


(10.12) 


It  follows  in  turn  from  Eq.  (10.11)  that  Eq.  (10.5)  can  be  expressed  in  the 


(j  kj2  fl 

E(r)=  -T-0-j  d«h* 


jk  z|l-Bu 
h  e  co(p’ 


(10.13) 


The  evaluation  of  cQ(p,  4>;  p>h)  through  Eq.  (10.12)  and  the  subsequent 
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evaluation  of  E(r)  through  Eq.  (10.13)  are  accomplished  numerically.  Un¬ 
fortunately,  for  high  accuracy  these  computations  require  a  very  large  amount 
of  computer  time,  particularly  that  of  Eq.  (10.12). 

In  the  inverse  Fourier  transformation  of  the  spectral  fields 
^3^(fh>  +  of  Section  2,  i.e.,  those  that  exist  in  the  absence  of 

the  scatterer,  the  point  in  space  at  which  the  fields  are  to  be  evaluated  is 

the  antenna  location.  In  the  coordinate  system  used  in  our  problem,  the 

/____ _ 

XA  +  >V  2  3 

located  near  the  point  x  3  0,  y  =  0.  For  that  reason,  the  exponent  in 
Jkoehpcos(<f>6"<*>) 

e  is  not  necessarily  large  in  those  calculations.  In  fact,  it 

is  zero  if  the  antenna  is  placed  exactly  at  the  origin.  Even  if  it  were 
assumed  that  the  factor  E( ,  <f>s'  +  q)  on  the  RHS  of  Eq.  (10.9)  is  "slowly 
varying"  in  the  angle  4>g  (which  would  seem  to  be  a  justifiable  assumption  be¬ 
cause  of  the  near-symmetry  of  the  circular  cable  ("near"  because  of  the  cable 
attenuation;  without  that  attenuation  the  same  amount  of  energy  would  be 
leaking  out  of  all  the  slots  in  the  cable.  The  attenuation  around  the  cable 
is  sufficiently  small  so  that  the  distribution  of  energy  from  different  slots 
is  quite  close  to  uniform)  the  stationary  phase  method  or  some  variation  there¬ 
of  would  not  be  easy  to  justify.  Hence,  the  full  numerical  integration  tech¬ 
nique  is  used  in  the  evaluation  of  the  field  components  in  the  absence  of  the 
scatterer. 

In  the  case  of  the  superposition  of  the  direct  scattered  field  and  the 

ground- reflected  scattered  field,  [(E^(Sh)  +  E^(Bu)]  in  Section  2,  .he 

-v  -v  ->n 

Fourier  transformation  is  taken  at  the  center  of  the  scatterer,  which  is 
located  near  the  cable.  Hence,  In  this  case 

p  *  p$  s  +  y\  -  R  +  &R  ,  where  |<5/?|  «  R  (10.14) 
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The  factor  (kQ  p$  8^)  In  the  exponent  in  Eq.  (10.12)  is  large  enough  to 

justify  a  stationary  phase  approximation  provided  8h  is  at  least  as  large  as 

0.01  and  provided  we  can  assume  that  E(8h,  +  <j>)  is  "slowly  varying"  in 

There  is  marginal  justification  for  this  assumption  in  this  case,  since  the 

scattered  field  may  be  significantly  variable  with  <j>g.  However,  it  is  still 

likely  that  it  does  not  vary  in  a  highly  oscillatory  fashion  with  $>g  and  that 

jk  p  8hcos(j)' 

the  factor  multiplying  e  p  is  still  "slowly  varying"  compared  with 

that  exponential  factor.  Hence,  the  stationary  phase  method,  based  on  the 
assumption  that  the  integrand  consists  of  a  slowly  varying  factor  and  a  highly 
oscillatory  exponential  factor  will  be  invoked  to  approximate  the  <j>g  integral. 
The  major  contribution,  in  this  case,  comes  from  the  stationary  phase  points. 

Because  of  time  and  financial  limitations  in  the  execution  of  this  project, 
it  is  not  possible  to  carry  out  an  extensive  study  to  determine  the  validity 
of  the  stationary  phase  approximation  for  this  problem.  Further  studies  might 
dictate  a  full-scale  numerical  integration  over  $g,  which  was  originally 
planned.  Such  an  integration  is  feasible,  but  preliminary  studies  indicate 
that  it  would  require  a  very  large  expenditure  of  computer  time  to  accomplish. 
Estimates  indicate  that  this  would  exceed  the  project’s  resources,  both 
temporal  and  financial.  As  a  compromise,  one  might  consider  a  numerical  inte¬ 
gration  over  a  small  range  of  angles  near  =  0  and  u  (which  are  the  stationary 

D 

phase  points)  which  seem  intuitively  to  be  the  directions  in  which  most  of  the 
incident  wave's  energy  would  propagate.  Even  this  approach  would  be  much 
more  computer-time  intensive  than  stationary  phase. 

Originally,  we  considered  a  FFT  algorithm.  After  a  careful  perusal  of  the 

problem.  It  was  decided  that  this  would  not  be  much  help  because  its  primary 

2 

utility  is  in  reducing  the  number  of  computation  points  from  N  to  N  log  N  in 
computations  of  the  discrete  Fourier  transform  at  N  points.  In  the  present 


computation,  the  Fourier  transform  is  evaluated  at  only  one  point  (the  scat- 
terer  center  position).  A  two-dimensional  Fourier  transformation  on  6X  and 
gy  is  performed  at  that  point,  but  because  Bh  =  +  By  is  between  0  and  1, 

the  polar  coordinate  form  (rather  than  the  rectangular  coordinate  form)  is  the 
natural  one,  and  when  cast  in  this  form,  the  problem  reduces  essentially  to  a 
single  Fourier  transform- type  integral,  that  over  the  angle  and  a  subse¬ 
quent  integration  over  Bh-  The  latter  is  not  necessarily  a  Fourier  transform- 
type  integration,  because  z%  is  so  small  and  might  even  vanish,  in  which  case 
this  last  step  reduces  to  an  ordinary  integration  problem  which  has  nothing 
whatsoever  to  do  with  Fourier  transformation. 

The  stationary  phase  technique  in  this  case  involves  the  differentiation 
of  the  phase  p  cos  <#>g  in  Eq.  (10.12)  with  respect  to  and  the  setting 
of  this  derivative  to  zero  in  order  to  find  the  stationary  phase  points. 


(cos  *£) 


sin  3  0 


*  0  or  tt 


For  4>g  near  zero  and  8h  t  0 

jk0pBhcos^  jk0pBh{l-[($^)/2]2} 
e  =  e 


(10.15) 


(10.16a) 


For  <f>g  near  it  and  8h  f  0 

(^  =  TT  +  A<pg ,  |A$g|  «  tt) 

Vc,peha-[Ca^)/232>  (1(Ua>) 


From  Eqs.  (10.16a,b)  and  (10.12), 


B fcnwfita*  al***-»  -  - 
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Cq(p»  0,  8^)  s  2^  ® 


LZi.  .1  Apsh  r  ,  /jkop6h(*e/2) 


d  e 


•  /o\2 


1  _  -Jkftp8h  r  jk0PSh(A<{»a/2) 

^  E(8h,  7T  +  0)  e  0  h  j  d(A»*)  e  0  hV  ^ 


/ — o !  ~  jknp6h  f00 

^r\s(E(8h-  ♦)•  T 


-jknpSh  r»  iy2 
+  E( 0^ ,  it  +  0)  e  I  dx  e^  } 


(10.17) 


f  dx  e+J'x2  =  2{|  dx  cos(x2)  +  j°°  dx  sin(x2)}  =  j\  { 1  +  j) 


(10.18) 


Substitution  of  Eq.  (10.18)  into  Eq.  (10.17)  results  in 


c0(p.  0»  8^) 


1  -  j(k0P6u“|-) 

] - (E(6h,  0)  e  0  h  4 


/2ir  k0  P  8h 


+  E($h,  it  +  0)  e 


•j(kop6h‘l)1  .. 


)  if  6h  t  0 


(10.19a 


If  8h  *  0,  then  0g  is  arbitrary.  Setting  0g  equal  to  zero  in  Eq.  (10.12),  we 


cjp,  0;  0)  *  E(0,  0)  if  &u  *  0 


(10.19b 


ti-Timr: 


The  numerical  integration  process  is  now  reduced  to  a  single  rather  than 
a  double  integration.  In  practical  computation,  we  must  choose  a  small  value 
of  8^»  which  we  may  call  8^.  below  which  8^  may  be  approximated  as  zero, 
(this  may  be,  for  example,  at  8h  *  0.01),  i.e.,  Eos.  (10. 19a, b)  may  be  re¬ 


written  as 


co(p»  ♦*»  8h>  “  — ,  =Z  I  (E(sh,  <0)  ej(k°P6h'?) 


^  Ko  p  8h 


”j(kftp6h-T) 

+  E(8h,  d>  +  ?r)  e  0  )  if  &.  > 


h  -  8ho 


*  E{0‘  ♦>  if  8h  <  8ho 


(10.20) 


Using  Eq.  (10.20)  in  Eq.  (10.13),  we  obtain 


E(r)  = 


(j  k J 


^  (E(0,  «)  fx*  d  8h  sh  ejk“Z<R 


‘”K0P  e 


f  1  sh  ^  S(Bh,  *)  eJ"'°PB'r4) 


?/  .  -a(k  p6h-J) 

+  E(eh ,  0  +  ir)  e  0  n  4  ]  e  0  '  h} 


(10.21) 


But  (since  8ho  «  1) 
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ho 


d  (ih  r-}i  o 


i  y 


» ..  .  -i  (.VV’-'W  /,  . 

J  ko  2  ^ 


,2  -  .  „JV 

no 


-  c  ”  ) 


i  jk  y.h-(\  jk  2  Bf'  jk  /. 

+  —  ■*  [e  0  «  ho  -  e  0  ]  -  e  0 

(j  K  *)'  l 

'  J  0 


,2 


(10. 22) 


Substituting  Fq.  (10.22)  into  Eq.  (10.21),  we  obtain  the  final  approxi- 

rule  form  of  E(r),  i.e., 

•>  -> 


E(r) 


(j  kj 

"  2  it 


2 

o'  _  {(‘‘ho,2  r/A  ^  _JV 


2 


■)-  E(0,  t)  o 


1 


+  . 1  - 

/2”  V  •>  L 


d  ph  /fih  W*V  ^  e 


j(koPiV^ 


ho 


+  E(f  ,  ,  <|>  +  ii)  e 


(10.23) 


The  integration  on  ;’j  indicated  in  Eq.  (10.23)  is  carried  out  on  the 
computer  using  a  Simpson’s  rule  algorithm.  In  the  case  of  the  scattered 
fields,  the  parameters  p  and  ij,  are  those  of  the  scatterer  center 


P  "  f>s  =/x \  +  yf  ;  <J>  "  <!'s  =  tan_1(  )  (10.24) 

*  In  the.  numerical  computations,  only  one  of  the  two  terms  in  (10.23) 

is  used  iri  any  single  computation.  The  first  term  wherein  <f,  =  <f$,  is  used 

for  scatterer  positions  outside  the  cable  and  the  second  term 
wherein  t  +  it,  is  used  for  positions  inside  the  cable. 

This  choice  issbased  on  intuition  and  the  knowledge  that  the  predominant 
contributions  come  from  the  portion  of  the  cable  near  the  scatterer,  and 
from  values  of  the  propagation  vector  R  nearly  parallel  with  the  line- 
of-sight  from  the  source.  'y 
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11.  NUMERICAL  RESULTS  AND  CONCLUSIONS 


11.1  Computer  Problems 

The  numerical  results  for  a  number  of  cases  that  were  run  on  the  North¬ 
eastern  University  Faculty  VAX  are  shown  in  Subsection  11.4.  There  have  been 
extensive  computation  problems  due  to  the  transition  from  the  CDC  Cyber  70 
to  the  VAX.  Among  these  were:  overflow  problems  due  to  the  reduced  range  of  the 
VAX;  reprogramming  required  to  adapt  programs  already  running  on  the  Cyber  70 
to  the  new  machine,  thus  requiring  extensive  debugging  time;  hardware  errors 
on  the  VAX  which  have  resulted  in  considerable  down-time;  and  the  necessity  to 
use  a  batch  processing  system  which  severely  limits  the  number  of  cases  that 
can  be  run  per  day. 

Because  of  these  computer  transition  problems,  which  were  beyond  the 
control  of  the  technical  staff  working  on  this  project,  there  has  been  a  delay 
of  several  weeks  in  the  execution  of  production  runs  on  the  computer. 

The  computations  performed  on  this  problem  are  quite  extensive  and 
require  considerable  running  time.  For  this  reason,  only  a  limited  number  of 
cases  were  computed  for  presentation  in  this  report.  The  computer  program  is 
available  and  could  be  used  to  study  many  more  parameter  regimes  that  might  be 
of  practical  interest. 

11.2  Program  Variables 

For  the  results  presented  in  this  report,  the  following  variables  had 
fixed  value  ranges  based  on  use  of  a  specific  cable  configuration  and  some 
empirically  determined  values  of  certain  parameters  given  to  us  by  the  Contract 
Monitor. 
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Parameters  associated  with  the  cable: 


Fortran 

Name 

Algebraic 

Symbol 

Definition 

Value 

i 

R 

Radius  of  cable 
configuration  in 
meters 

24 

(approximation 
based  on  circum¬ 
ference  of  151 
meters ) 

A 

a 

Inner  radius  of 
coaxial  cable- 
meters 

.00476 

B 

b 

Outer  radius  of 
coaxial  cable- 
meters 

.0127 

Length 

L* 

Slot  length  (along 
cable);  same  for 
all  slots-meters 

.003 

Width 

\ 

Slot  width  (around 
cable);  same  for  all 
slots-meters 

.0155 

EPSILONCA 

eca 

Permittivity  of 
cable  material  - 
farads/meter 

15.05(10"12)  * 

SIGMACA 

°ca 

conductivity  of 
cable  material  - 
mhos/meter 

0 

NSLOTS 

Ns 

Number  of  slots 
on  cable 

31,723 

PHI BAR 

1 

*10 

Azimuthal  angle 
along  the  cable 
of  Slot  #1.  (The 
slot  nearest  the 
power  source)  in 
ground  frame 

0 

PHI 

<p .  "  =  <t>  " 

Angle  of 
slot  center 
around  periphery 
of  cable 

— J  *  45® 
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Fortran 
_ Name 


AT gebraic 
Symbol 


Definition 


Value 


NMODES 

Number  of  modes 
of  each  category 

1  TEM  mode 
(quasi -TEM  or 
coaxial  mode 
observed  to  propa¬ 
gate  in  this  cable 
as  the  principal 
mode) 

ALTEM 

%  (TEK) 

Attenuation  of 

TEM  mode-nepers/ 
meter 

.002 

ALTEN 

%  'TE! 

Attenuation  of 

TE  mode-  nepers/ 
meters 

— 

ALTMN 

on(TM) 

Attenuation  of 

TMn  mode -nepers/ 

meter 

— 

F 

f 

Radio  frequency 
in  Hertz 

57(1 06)= 

57  MHz 

ATEN 

A<TEn> 

Amplitude  of  TE 
mode-units  not  n 
indicated  here,  be¬ 
cause  this  is  treat¬ 
ed  as  a  scale  fact¬ 
or  and  its  units 
are  not  important 

— 

ATMN 

A<™n> 

Amplitude  of  TM 
mode-same  re-  n 
marks  as  above 
concerning 
units 

Fortran 

Name 

Algebraic 

Symbol 

Definition 

Value 

VO 

Vo 

Voltage  between 
inner  and  outer 
convectors  of 
coaxial  cable- 
used  in  amplitude 
of  TEM  mode- volts 

1 

E 

e 

Permittivity  of 

ground-farads/ 

meter 

35 .416(1 0-1 ^ )  * 

4en 

0  » 

SIGMA 

a 

Ground  conductivity 
-mhos/meter 

0.002 

XAP 

XA 

x-coordinate  of 
antenna-meters 

0 

(antenna  at  center 
of  circular  cable 
configuration) 

YAP 

y-coordinate  of 
antenna-meters 

0 

(same  remarks 
as  above) 

ZAP 

ZA 

z-coordinate  of 
antenna-meters 

0.5 

BH 

6h 

Magnitude  of  (x-y) 
plane  projection 
of  B  , 

Varied  from 

0  to  1 

i.e.  VBx  +  By 

PHIB 

*6 

Azimuthal  angle 
of  (x-y)  plane  pro- 

Varied  from  0° 
to  360° 

jection  of  B  .  i.e. 


Parameters  Associated  with  the 

Scatterer 

Fortran 

Name 

Algebraic 

Symbol 

Definition 

Values 

XSP 

XS 

x -coordinate  of 
scatterer  center- 
meters 

Varied  from 
case  to  case 

YSP 

y-coordinate  of 
scatterer  center- 
meters 

Varied 

ZSP 

zs 

z- coordinate  of 
scatterer  center- 
meters 

Varied 

THB 

0b 

Spherical  polar 
angle  of  long 
axis  of  prolate 
spheroidal  scatterer 
(in  ground-frame) 

Varied: 

0  for  "upright 
man"  cases; 

-y-  for  "crawling 

man"  cases 

PHB 

*b 

Azimuthal  angle 
of  long  axis  of 
scatterer  (in 
ground-frame) 

Varied 

ES 

£s 

Permittivity  of 
scatterer  material - 
farads/meter 

408.992(10”^ )  = 
46.2c0= 

%<60£o>. 

based  on 
Reference  D. 19 

SIGMAS 

as 

Conductivity  of 
scatterer  material- 
mhos/meter 

0.592 

based  on 

Reference  D.19 

RS 

(computed 

variable) 

rs 

y/xsV  =  radia1 

coordinate  of 
scatterer  center 
in  ground  frame 

Varied 

PHIS 

(computed 

variable) 

*s 

tan_1/ys  \ 

1— — J-  azimuthal  Varied 

angle  of  scatterer 

center  in  ground  frame 


11-5 


Fortran 

name 


Algebraic 

Symbol 


Definition 


Values 


RS 

Rs 

Radius  of  prolate 

spheroidal 

scatterer-meters 

0.25 

LS 

Ls 

Length  of  prolate 

Varied 

spheroidal 

scatterer-meters 
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11.3  Delineation  of  Cases  Studied 


The  outputs  of  the  computations  are: 


amplitude  of  x-component  of  the  electric  field  of  the  signal  at 


the  antenna  where  q  =  o  for  the  field  in  the  absence  of  the  scatterer,  q  =  s  for 
the  scattered  wave  field  and  a  =  (0+s)  for  the  superposition  of  the  field  without 
the  scatterer  and  that  due  to  the  scatterer. 

lEy^i  =  amplitude  of  y  component  of  field  at  antenna 
I |  ■  amplitude  of  z  component  of  field  at  antenna 

iE<q)i  -  V|E*(qii2  *  |Ey<,)|2  *  |ez(9)i2 


The  parameters  to  be  varied  in  these  computations  are: 
x$  =  x  component  of  scatterer  center-meters 
y$  =  y  component  of  scatterer  center-meters 
zs  =  z  component  of  scatterer  center-meters 
L$  =  Long  dimension  of  scatterer  =  Height  of  "man"-meters 
P-s  =  Height  of  "stilts"  in  cases  where  the  intruder  is  assumed  to  be 
elevated-meters 

*s  "  tan  As\a  Azimuthal  angle  of  scatterer  center 


Ax$  *  xs  "  xc  *  where  xc  *  x  coordinate  of  cable-meters 

Ays  *  y$  -  yc  ,  where  yc  =  y  coordinate  of  cable-meters 

Ar$  ■  r$  where  r$  =  r  coordinate  of  scatterer  center-meters 
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The  first  set  of  results  which  are  of  course  the  same  for  all  scatterer  positions 


and  orientations,  are  the  field  amplitudes  of  the  antenna  in  the  absence  of  the 
scatterer.  These  are  as  follows:* 

|Ei0>l  *  4.8166{10'4) 

|E^o)|  *  2.38S8O0'2) 

]E^o)|  =  5.0226(10'“) 

The  computed  field  in  the  absence  of  a  scatterer  is  predominantly  in  the 
y-direction.  The  interpretation  of  the  result  will  be  discussed  in  Subsection  11.4.6. 
The  cases  involving  an  intruder  are  delineated  below. 

Cases  1  :  "Upright  Han"  Model 

These  cases  involve  a  model  of  an  intruder  walking  in  upright  position  across  the 
cable.  The  intruder  starts  at  a  point  outside  the  cable  and  walks  radially  Inward, 
crosses  the  cable  and  continues  to  walk  toward  the  center  of  the  configuration. 

P<*cause  of  time  limitations  computations  were  made  for  only  a  few  positions 
on  either  side  of  the  cable.  Since  the  effects  should  not  be  entirely  circularly 
symmetric,  due  to  cable  attenuation  and  other  effects,  these  computations  were  made 
at  four  different  azimuthal  angles. 

Case  I-fl 

"Upright  man",  walking  radially  inward  across  cable;  0^  *  0°,  =  0°,  r$  varied 

from  28  to  17.5  meters  in  stens  of  1.5  meters,  zs?Ls/2(see  diaoram  belnw);  *s=  n« 
path  along  x  axis  in  -x  direction). 

*  These  "field  amplitudes"  are  actually  relative  quantities  and  are  not  intended  to 
represent  the  actual  field  amplitudes.  If  the  true  values  of  the  power  generated 
by  the  source  in  the  cable  were  used  here,  these  numbers  would  be  true  field  ampli¬ 
tudes  in  volts  per  meter. 
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Case  1-A  (continued) 


Side  View 


Downward  View  * 


Path  of  intruder 


To  antenna 


Ground ' 


Cable  cross-section 


Path  of 
r  intruder 


Case  I-B;  same  as  I-A,  but  $  *  it  ;  path  along  y  axis  in  -y  direction 

*> 


Downward  View 


28  m 


Case  I-C,  same  as  I-A  and  I-B  but  =  7r;  path  along  x  axis  in  +x  direction 


Downward  View  y 

*-.28  m  — 

Path  of  ' 

intruder  /  U 


kSMEI 


Case  I-D:  same  as  I-A,  I-B  and  I-C,  but  a  =  3ti  ;  path  along  y  axis  in  +y 

s  4 

direction. 

Cases  II;  "Elevated  Upright  Man”  Model 

In  these  cases  the  intruder  is  again  upright,  but  elevated,  i.e.,  on  stilts. 

The  stilts  are  assumed  to  scatter  no  energy;  hence  the  only  change  in  the  mathematical 
model  is  that  the  scatterer  center  is  raised  (relative  to  Cases  I)  by  an  amount  zs, 
the  height  of  the  stilts. 

These  cases  are  evaluated  with  the  intruder  at  a  few  fixed  positions,  to 
determine  at  those  positions  the  effect  of  various  elevations  on  the  field 
components  seen  at  the  antenna. 

Case  I  I-A  Upright  man  on  stilts;  ©b  =  0°  ,  4>b  =  0°;  *  Height  of  stilts; 

z$  *  +  l_s  ;  Zs  varied  from  0.4  to  2.0  meters;  xs  *  22,  ys  *  0. 


Case  I I-B.  Same  as  Case  I I-A  but  x$  =  23.5 
Case  I I-C .  Same  as  Cases  I I-A  and  I I-B  but  x$  =  25 
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Case  III:  "Upright  Man",  standing  at  ground  level;  ©b  =  0°  ,  $b  =  0°  ,  zs  =  L$  , 

T 

y$  =  0  ;  Ls  varied  from  1.6  to  2.08.  This  case  is  for  the  purpose  of 
examining  the  variation  of  the  field  components  seen  at  the  antenna  with  the  height 
of  the  intruder.  Case  III-A,  x_  =  22;  Case  III-B,  xc  =  23.5;  Case  III-C,  xe  =  25. 

^  i  5 

Case  IV:  "Crawling  man"  Models 

In  these  cases  the  intruder  is  "crawling"  radially  inward  from  a  position 
outside  the  cable  to  a  position  inside  the  cable.  The  spheroid  in  the  model  is  lying 
on  the  ground  surface  with  the  long  axis  pointed  in  the  radial  direction. 


Case  IV-A:  "Crawling  man",  ©b  =  I  .  4>b  =  0°,  z^R^,  <J>s  =  0,  rs  varied  from  28  to 
17.5  meters  in  steps  of  1  meter;  path  same  as  in  Case  IA. 


Downward  View 


Scatterer 
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CASE  IV-B  Same  as  Case  IV-A  except  that  4>h  =  tt.  ,  path  same  as  in  Case  I-B. 

D  2 


Downward  View 


Case  IV-C:  same  as  Cases  IV-A  and  IV-B  except  that  ^  =  *  •>  path  same  as 


Case  I-C. 


Case  IV-Q:  same  as  Cases  IV-A,B,C  except  that  =  -rp,  path  same  as  Case  I-D. 


Case  V:  Prone  Position  Transverse  to  Radial  Direction 
In  these  cases,  the  intruder  is  in  a  prone  position,  but  the  long  axis 
is  transverse  to  the  radial  direction.  Otherwise,  same  as  "crawling  man"  cases 
In  cases  V- A,  V-B,  V-C  and  V-D,  =  0,  VZ  ,  71  and  3*/2,  respectively. 


•  /* *y r **;*+"+ it ■  ■  **-»^**  i 


In  addition  to  Cases  I  through  V  above,  two  other  sets  of  auxiliary  computations 
were  made.  Their  purpose  was  to  determine  the  magnitude  and  polarization  of  the 
fields  illuminating  the  scatterer  at  the  positions  assigned  to  its  center  in  Cases  I 
through  V. 

The  first  set  of  computations  were  of  the  field  components  due  to  the  cable 
without  ground  reflections  at  the  values  of  p  and  $  corresponding  to  the  scatterer 
center  values  p$  =  17.5,  19,  20.5,  22,  23.5,  25,  26.5  and  28  meters.  Four  sets  of 
these  runs  were  made,  at  $  =  0°,  90°,  180°  and  270°.  These  values  of  p  and  <t> 
correspond  to  the  scatterer  positions  in  Cases  I-A,  B,  C,  D,  IV-A,  B,  C,  D  and 
V-A,  B,  C,  D. 

The  second  set  of  these  illuminating  field  computations  were  identical  with  the 
first  but  with  the  inclusion  of  ground  reflections,  i.e.,  in  the  first  set  of 
computations  the  ground  permittivity  is  that  of  free  space  and  the  ground  conductivity 
is  zero,  those  assignments  being  equivalent  to  the  removal  of  the  effect  of  the 
ground.  In  the  second  set  the  values  of  qround  permittivity  and  conductivity  were 
those  indicated  in  the  table  in  Section  11.2  above.  Comparison  of  the  first  and 
second  sets  of  results  should  allow  a  determination  of  the  effect  of  ground  reflections 
on  the  fields  illuminating  the  scatterer. 


11.4  Presentation  and  Discussion  of  Numerical  Results 


The  available  numerical  results  will  be  presented  In  this  section,  beginning 
with  the  results  of  the  study  of  the  fields  illuminating  the  scatterer  (Subsection 
11.4-1),  followed  by:  (Sec.  11.4-2)  the  "upright  man"  (or  equivalently  "radial 
walk")  cases  (Cases  1-A,  b,  C,  U);  the  “elevated  upright  man"  or  “man  on  stilts" 
cases  (Cases  II-A,  8,  C;  Subsection  11.4-3);  the  "height  variation"  cases  (Cases 
III-A,  B,  C;  Subsection  11.4-4)  and  the  "crawling  man"  cases  (Cases  IV-A,  B,  C,  0  and 
V-A,  B,  C,  D;  Subsection  11.4-5). 

In  discussing  the  results  we  have  tried  to  take  cognizance  of  available 
experimental  results  whenever  they  exist.  Those  experimental  results  mentioned  in 
what  follows  are: 

Reference  B.2,  Figure  5,  Page  10  (circumferential  walk) 

Figure  8,  Page  13  (  “  ) 

Reference  B.3,  Figure  9,  Page  15  (circumferential  walk) 

Figure  10,  Page  16  (radial  walk) 

Figures  13,  15,  16,  17,  Pages  18,  20,  21,  22  (circumferential  walk) 

Reference  B.6,  Figure  7,  Page  12  (circumferential  walk) 

Figures  4,  5,  6,  Pages  1,  12  (circumferential  walk;  analytical 
results  based  on  L.  Poirier's  simplified  analysis  neglecting 
field  polarization) 

In  all  of  the  circumferential  walk  cases  studied  experimentally,  a  very  sub¬ 
stantial  oscillation  is  observed  in  the  results  as  the  intruder  walks  around  the 
periphery  of  the  cable.  This  is  also  exhibited  by  L.  Poirier's  analysis  reported 
in  Reference  B.6.  Roughly  speaking  the  oscillatory  behavior  is  due  to  alternations 
between  constructive  and  destructive  interference  as  the  scattering  object  changes 
azimuthal  angle. 

The  experimental  curves  plotted  in  the  above  references  are  of  the  total 
signal  voltage  received  at  the  antenna  vs.  $s,  the  scatterer’s  azimuthal  angle. 

These  are  compared  with  the  signal  voltage  in  the  absence  of  the  scatterer  (quiescent 
level).  The  signal  voltage  (assuming,  of  course,  that  the  antenna  is  linear)  is 


proportional  to  the  field  strength  of  the  component  polarized  in  the  direction  to 
which  the  antenna  is  receptive.  In  the  case  of  the  monopole  antenna  used  in  making 
_  the  measurements,  that  direction  is  the  vertical.  Hence,  only  our  results  on  |EZ|, 
the  amplitude  of  the  vertical  component  of  the  field  at  the  antenna,  could  be  com¬ 
pared  with  those  results.  The  horizontal  components  j E  |  and  )E  |  cannot  be  so 

x  y 

compared,  since  the  antenna  is  not  responsive  to  horizontal  fields. 

Two  points  must  be  made  before  proceeding  with  the  numerical  results.  First, 
it  was  impractical  under  our  time  limitations  to  plot  circumferential  walks;  con¬ 
sequently  all  results  where  the  scatterer  moves  and  also  those  results  for  the  fields 
illuminating  the  scatterer  at  various  positions  are  presented  as  "radial  walks",  i.e., 
the  field  components  are  presented  for  various  radial  positions  p$  at  a  fixed  azimuthal 
position  Moreover,  the  angles  $$  are  confined  to  0°,  90°,  180°,  and  270°.  It 
would  be  highly  desirable  to  have  results  for  many  intermediate  angles.  There  was 
insufficient  time  to  compute  enough  points  around  the  cable  to  obtain  meaningful 
"circumferential  walk"  results.  To  do  so  would  require  computations  for  every  few 
degrees(e.g.  2°  spacing  would  require  180  computations  for  each  fixed  value  of  p$). 

For  this  reason  it  was  decided  to  confine  attention  to  the  radial  varation  at  fixed 
azimuth  angle. 

A  second  point  is  that  there  is  some  question  in  our  minds  (and  no  remaining 
tine  to  study  the  question  and  try  to  improve  the  model  to  rectify  the  problem  if  it 
exists)  on  the  accuracy  of  the  ratio  of  the  computed  scattered  fields  at  the  antenna 
(i.e.,  direct-plus-ground-reflected  scattered  fields)  to  the  computed  fields  at  the 
antenna  in  the  absence  of  the  scatterer  (direct-plus-ground-reflected).  According 
to  our  results  the  fields  at  the  antenna  (at  the  center  of  the  configuration; 
xA  =  yA  =  0,  zA  =  0.5)  in  the  absence  of  the  scatterer  are  (in  dB  relative  to  an 
arbitrary  reference  level;  the  same  reference  level  used  in  all  of  the  results  with 
end  without  scatterers  presented  in  this  report): 
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|Ex|  =  -66.3  dB 
)Ey|  =  -32.4  dB 
|EZ|  =  -65.9  dB 

These  results  will  be  discussed  in  Subsection  11.4-6. 

Certain  approximations  used  in  our  analysis  if  not  as  accurate  as  we  originally 
believed,  might  have  a  different  effect  on  the  antenna  fields  in  the  absence  of  the 
scatterer  than  they  do  on  the  scattered  fields.  If  indeed  that  is  the  case,  then 
at  the  very  least  there  might  be  a  scale  factor  error  in  the  results  without  the 
scatterer  present  or  alternatively  a  scale  factor  error  in  the  scattered  fields,  or 
possibly  in  both  sets  of  fields.  If  that  were  the  case,  then  coherent  addition  of 
the  fields  with  the  scatterer  present  to  the  fields  without  the  scatterer  present 
could  produce  totally  meaningless  results.  The  two  important  constituents  of  the 
analysis  are  the  fields  generated  by  the  cable  and  the  perturbation  of  the  field  by 
the  scatterer.  Rather  than  present  large  numbers  of  results  on  the  composite  fields 
(all  of  which  have  been  computed  for  each  case  studied)  which  may  be  greatly  in 
error  because  of  an  erroneous  scale  factor,  it  was  decided  to  concentrate  on  the 
direct-plus-ground-reflected  scattered  field  components  and  show  their  behavior  as 
a  function  of  various  parameters. 


11.4.1  Numerical  Results  on  Illuminating  Fields 

The  first  results  presented  are  those  of  the  field  components  illuminating 
the  scatterer.  These  appear  in  tabular  form  in  Tables  11.1,  11.2  and  11.3,  Columns 
1  through  8.  Some  of  the  results  are  shown  graphically  in  the  curves  of  Figures  11.1 
through  11.16. 

Results  are  shown  both  with  and  without  the  effects  of  ground  reflections. 

The  format  for  each  of  the  curves  shown  is  a  plot  of  a  field  component  in  dB 
relative  to  an  arbitrary  reference  level  against  p,  the  radial  coordinate,  from 
p  =  17.5  to  p  =  28  meters.  The  reference  level  has  no  significance  because  there 
are  arbitrary  factors  which  cause  the  computed  field  component  magnitude  to  differ 
from  the  actual  field  magnitude  in  volts  per  meter.  All  plots  in  these  figures 
and  also  those  in  subsequent  figures  are  in  dB  relative  to  the  same  reference  level, 
hence,  it  is  the  relative  magnitudes  of  field  strength  rather  than  their  absolute 
magnitudes  on  which  attention  should  be  focussed.  Since  the  cable  is  at  a  radial 
distance  of  24  meters,  the  plots  may  be  considered  to  be  from  a  distance  of  6.5 
meters  inside  the  cable  to  a  distance  of  4  meters  outside  the  cable. 

F-om  very  simple  reasoning  based  on  a  single  plane  wave  propagating  from  the 
nearest  cable  slot  to  the  observation  point,  we  would  expect  something  close  to  a 
decay  in  amplitude  proportional  to  the  inverse  distance  from  the  nearest  cable  slot. 
The  coherent  addition  of  contributions  from  different  parts  of  the  cable,  with  the 
attendant  variability  in  the  relative  phase  etween  contributions  from  different 
cable  slots,  would  preclude  exact  inverse  distance  behavior.  However,  we  look  for 
some  suggestion  of  such  a  variation  with  distance.  Unfortunately  we  do  not  always 
attain  behavior  close  to  the  above. 


Focussing  for  the  monent  on  the  "best"  results  in  terms  of  correlation  with 
expectations  based  on  simple  reasoning,  we  examine  the  results  in  Table  11.3  and  in 


Figures  11.2,  11.4,  11.7  showing  the  z-components  of  *  '  field  as  a  function 
of  radial  distance  without  ground  reflection. 

At  $  3  0°,  the  field  amplitude  decays  from  a  value  of  -5.31dB  at  a  radial 
distance  of  0.5  meters  from  the  cable  to  a  value  -39dB  lower  at  6.5  meters  inside 
the  cable  and  to  a  value  -12dB  lower  at  4  meters  outside  the  cable.  The  inverse 
distance  behavior  would  give  us  a  decay  of  roughly  -22dB  at  6.5  meters  inside  and 
-18dB  at  4  meters  outside.  Thus  the  rate  of  decay  is  considerably  larger  than  would 
be  dictated  by  an  inverse  distance  law.  This  could  be  at  least  partially  due  to 
the  importance  of  the  approximately  inverse  squared  contribution  at  a  distance  very 
close  to  the  scatterer. 

The  z-field  amplitudes  are  not  changed  appreciably  when  ground  reflections 
are  added  at  4  =  0°,  The  z-field  amplitudes  are  slightly  smaller  but  within  about 
2  or  3dB  of  their  corresponding  values  without  ground  reflections  and  the  trends 
are  qualitatively  and  quantitatively  the  same.  Since  interference  between  direct 
and  reflected  waves  can  be  either  constructive  or  destructive,  there  is  no 
certainty  that  the  ground  reflection  contribution  will  increase  or  decrease  the 
field  magnitude.  In  this  case,  at  the  reflection  angles  at  which  these  effects 
occur,  the  interference  appears  to  be  mostly  destructive. 

At  <|>  =  90°  (Table  11.3  and  Figures  11.9,  11.11,  11.14)  the  trends  in  the  z- 
field  are  qualitatively  similar  but  the  rates  of  decay  with  radial  distance  from 
the  cable  are  not  as  pronounced.  The  decay  between  0.5  and  6.5  meters  inside  the 
cable  and  that  between  0.5  and  4  meters  outside  the  cable  are  respectively  19 . 6dB 
and  16.9dB  without  ground  reflections  and  24 . ldB  and  22 . ldB  respectively  with 
ground  reflections.  These  values  are  much  closer  to  the  respective  22dB  and  18dB 
that  would  prevail  for  an  inverse  distance  decay  law. 
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the  decay  levels 


At  =  180°  (Table  11.3  and  Figure  11.14) 
corresponding  to  the  above  comparisions  are  17.3dB  and  17.5dB  respectively  without 
ground  reflections  and  2l.9dB  and  19. 6dB  respectively  with  ground  reflections. 

Again  these  values  are  within  a  few  dB  of  what  would  be  expected  for  an  inverse 
distance  dependence.  Again,  as  at  41  =  0°  and  <f>  *  90°,  the  z-field  magnitudes  with 
ground  reflection  tend  to  be  slightly  lower,  usually  by  less  than  2dB,  than  their 
counterpart  values  without  ground  reflection,  indicating  mostly  destructive  inter¬ 
ference  between  direct  and  ground  reflected  waves. 

The  z-field  results  for  $  =  270°  (Table  11.3)  appear  anomalous  with  respect 
to  the  above  considerations.  The  same  comparison  exists  between  results  with  and 
without  ground  reflections  (nearly  the  same  in  both  cases),  but  there  is  no  obvious 
explanation  for  the  strange  behavior  of  the  z-field  magnitude  outside  the  cable 
and  its  erratic  behavior  inside  the  cable. 

The  variation  in  the  average  z-field  amplitudes  as  the  angle  $  changes 
through  90°  increments  is  substantial,  being  nearly  30dB  in  some  cases.  This  is 
not  anomalous,  since  experimentally  observed  scattered  field  magnitudes  (See  e.g. 
circumferential  walk  results  in  Ref.  B.2;  Fig.  5,  Pg.  10;  Fig.  8,  Pg.  13,  Ref.  B.3; 
Figs.  9,  13,  15,  16,  17,  Pgs.  15,  18,  20,  21,  22;  Ref.  B.6;  Fig.  7,  Pg.  12)  exhibit 
a  highly  oscillatory  azimuthal  variation.  Since  nulls  exist  at  some  angles  no 
variation  in  dB  would  be  anomalous.  This  can  be  explained  at  least  partially  by  the 
highly  complicated  nature  of  the  superposition  of  contributions  with  widely  differing 
phases  from  various  portions  of  the  cable.  More  will  be  said  of  this  in  later 
discussions  of  the  azimuthal  variation  of  our  numerical  results  on  scattered  field 
ampl i tudes. 

The  horizontal  components  of  the  computed  illuminating  fields  are  not 
generally  as  well-behaved  as  the  vertical  components  with  respect  to  any  resemblance 


to  inverse  distance  behavior.  This  is  possibly  due  to  ....a  fact  that  the  waves  from 
different  slots  along  the  cable  not  only  add  coherently  with  complicated  relative 
phases,  but  also  add  vectorially.  In  all  cases,  the  y-compone’it  seems  to  predom¬ 
inate. 


11.4.2  Results  for  the  “Upright  Man11  Cases 

In  Tables  11.4,  11.5  and  11.6,  Column  1  through  4  and  in  Fiqures  11.17  -  11.23 
some  scattering  results  for  a  man  1.78  meters  tall  (about  5' 10")  and  about  0.5 
meters  wide  in  upright  position  proceeding  from  (or  to)  a  point  17.5  meters  inside 
the  cable  to  (or  from)  a  point  28  meters  outside  the  cable.  These  are  the  cases 
I-A,  B,  C,  0  delineated  in  Section  11.3,  where  A,  B,  C  and  D  refer  to  scatterer 
center  azimuthal  angles  =  0°,  90°,  180°  and  270°  respectively. 

The  first  point  to  be  made  is  that  all  of  the  field  component  amplitudes 
exhibit  a  very  large  decrease  (from  40  to  70dB)  at  $s  =  90°  and  270°  relative  to 
their  values  at  <j>$  =  0°  and  180°.  Magnitudes  at  180°  are  generally  between  20 
and  30dB  below  those  at  0°,  while  those  at  270°  and  those  at  90°  are  comparable 
at  some  positions  and  separated  by  10  to  20dB  at  other  positions.  None  of  these 
variations  appear  to  come  directly  from  attenuation  along  the  cable,  which  is  a 
maximum  of  2.62dB  around  the  entire  circumference  of  the  cable,  Ap  was  discussed 
in  Section  11.4.1  in  connection  with  illuminating  fields,  the  circumferential 
variations  in  fields  illuminating  the  cable  are  probably  at  least  partially  due 
to  the  wide  variations  in  phase  of  the  contributions  from  different  portions  of 
the  cable.  If  one  adds  to  that  the  phase  variations  in  the  scattering  process 
itself,  one  cannot  expect  the  scattered  fields  from  various  azimuthal  angles  to 
follow  a  smooth  variation  or  to  be  comparable  in  magnitude,  in  spite  of  superficial 
indications  of  symmetry.  If  all  cable  field  contributions  were  added  noncoherently , 
i.e.,  if  their  amplitudes  or  squared  amplitudes  were  added  rather  than  their 
complex  fields,  the  results  would  undoubtedly  show  much  greater  azimuthal  symmetry. 
Experimental  results  for  the  vertical  field  components  (Ref.  B.3,  Pg.  10;  Ref.  B.3, 
Pp.  15,  18,  20,  21,  22;  Ref.  B.6,  Pg.  12)  show  substantial  azimuthal  fluctuations. 
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Turning  our  attention  now  to  the  radial  variation  in  field  amplitude  at  a 


fixed  azimuthal  angle  most  of  the  results  in  Tables  11. a,  11.5  and  11.6  and  Figures 
11.17  through  11.23  show  a  substantial  decay  as  the  scatterer  recedes  from 
the  cable  outside  the  cable  but  a  flat  or  sometimes  slightly  increasing  response 
as  it  recedes  from  the  cable  inside  the  cable.  Also  the  peak  doesn't  always  occur 
at  23.5,  the  closest  radial  distance  to  the  cable  for  which  a  computation  was  made, 
but  sometimes  at  25  or  at  22  or  even  further  away  in  a  few  cases.  Since  computations 
were  made  at  intervals  of  1.5  meters,  the  actual  peaks  might  occur  between  the 
computed  points.  But  in  general  the  occurrence  of  the  peak  at  a  point  other  than 
the  point  directly  above  the  cable  is  consistent  with  the  same  observation  made 
experimentally  and  discussed  on  Page  16  of  Reference  B.3.  Figure  10  on  Page  16 
of  that  reference  (radial  walk)  illustrates  the  point,  wherein  peak  responses 
sometimes  occur  a  meter  or  more  away  from  the  cable  radius. 

Based  on  the  simple  reasoning  of  Section  11.4.1  we  should  look  for  two 
contributions  to  variation  of  field  amplitudes  with  radial  distance  in  the  present 
case.  First,  this  simple  reasoning  would  predict  a  variation  as  the  reciprocal  of 
the  separation  distance  between  the  nearest  cable  slot  and  the  scatterer  position 
due  to  the  expected  variation  in  the  illuminating  field.  Also  there  should  be  an 
inverse  radial  distance  dependence  due  to  the  factor  appearing  in  the  scattered 
field  component.  For  example, (Subsection  11.4.1)  we  should  expect  from  the  first 
of  these  mechanisms  about  22dB  of  decrease  in  amplitude  between  23.5  and  17.5  meters 
and  about  18dB  between  23.5  and  28  meters.  From  the  second  effect  we  should  expect 
a  2.56dB  increase  as  we  go  from  23.5  to  17.5  and  a  1.52dB  further  decrease  as  we 
proceed  from  23.5  to  28  meters.  This  means  that  the  net  decay  should  be  about  19 . 5dB 
from  23.5  to  28  meters  and  about  20 . 5dB  from  23.5  to  17.5  meters.  If  we  were  to  base 
our  interpretation  of  the  results  on  these  mechanisms  alone,  we  would  conclude  that 
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the  decay  of  between  about  20  and  25dB  that  we  observe  on  most  of  the  field  com¬ 
ponents  between  scatterer  positions  at  23.5  and  28  meters  is  quite  reasonable,  but 
that  the  flat  or  even  increasing  response  we  see  as  we  proceed  from  23.5  to  17.5 
meters  inside  the  cable  cannot  be  explained  by  these  mechanisms. 

There  is  another  feature  of  the  geometry  which  could  help  to  explain  the 
great  asymmetry  between  departures  from  the  cable  radius  inside  and  outside  the 
cable.  Outside  the  cable,  as  we  recede  further  from  the  cable  radius,  the  cable 
appears  convex  from  the  observation  point  and  the  distances  to  the  cable  slots 
other  than  the  nearest  one  become  progressively  greater  as  we  move  out.  As  we 
move  inward  inside  the  cable,  however,  the  cable  appears  concave  and  the  distances 
from  the  various  parts  of  the  cable  nearest  the  observer  become  more  nearly  com¬ 
parable,  leading  to  a  greater  tendency  toward  constructive  interference  between 
contributions  from  various  slots. 

Some  of  the  results  shown  in  Tables  11.4,  11.5  and  11.6  appear  anomalous, 
and  there  is  no  obvious  explanation  for  them  based  on  any  of  the  mechanisms  dis¬ 
cussed  above.  More  will  be  said  of  this  in  Subsection  11.5. 
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These  results  correspond  to  cases  II-A,  B,  C  in  Section  11.3.  The  "man"  is 
placed  on  stilts  in  three  different  fixed  positions  (<!>s  =  0°,  d$  =  22,  23.5  and  25, 
those  radial  positions  nearest  the  cable  at  which  fields  were  computed).  The  results 
showing  field  component  amplitude  vs  height  of  stilts  at  each  position  are  presented 
in  Tables  11,7,  11.8  and  11.9.  Selected  results  are  also  shown  in  Figures 
(p  =  22),  Figure  11.26  (pe  =  23.5)  and  Figures  11.27,  11.28(pc  =  25). 

In  all  three  cases,  the  variation  with  stilt  height  is  from  0  to  2.0  meters. 

The  field  component  amplitudes  in  all  three  coordinate  directions  (except  for 
the  vertical  field  component  at  p$  =  25m,  which  exhibits  anomalous  behavior)  there 
is  a  decrease  of  between  2dB  and  6dB  between  the  case  of  zero  stilt  height  and  the 
case  of  2  meter  stilt  height.  This  variation  is  in  the  proper  direction,  obviously, 
since  the  expectation  is  for  a  reduction  in  field  strength  as  the  scatterer  is 
further  elevated  and  therefore  moves  further  away  from  the  cable  slots  generating 
the  fields.  To  examine  the  issue  of  how  close  the  variation  is  to  the  inverse 
distance  variation,  we  note  that  the  difference  between  the  "no  stilts"  case  and 
2  meter  stilts  (since  the  "man"  is  1.78  meters  tall)  is  a  change  in  the  z-coordinate 
of  the  scatterer  center  from  0.890  to  1.890.  With  a  pure  inverse  distance  law, 
this  would  result  in  roughly  a  6dB  decay  between  these  two  cases.  This  is  not 
very  far  from  the  amount  of  decrease  observed  in  most  of  the  computed  results. 


11.4.4  Height  Variation 


A  set  of  results  showing  field  component  amplitude  variation  versus  intruder 
height  at  ♦  *  C°  and  a.  -  22,  23.5  and  25  meters  (Cases  1 1 1 - A ,  1II-B  and  II1-C 
respectively  in  Section  11.3)  are  shown  in  Tables  11.10,  11.11  and  11.12.  Selected  results 
also  appear  in  Figs.  11.29,  11.30,  Figs.  11.31,  11.32  and  Fiqs.  11. 33-11.35  respectively. 
In  each  case  the  height  is  varied  from  1.60  to  2.08  meters  (about  5 ' 3"  to  nearly 
6'10",  a  more-than-reasonable  range  of  heights  for  an  adult  male).  In  all  of  these 
cases,  an  increase  in  field  amplitude  of  about  6dB  is  observed  as  the  height 
increases  over  the  indicated  range. 

This  variation  is  certainly  in  the  proper  direction.  The  height  range 
represents  an  increase  of  about  30%,  if  we  want  to  consider  the  field  amplitude 
as  proportional  to  the  length  of  the  scatterer,  this  would  imply  an  increase  of 
about  2.27dB  between  the  response  of  the  short  fellow  of  5' 3"  and  the  basketball 
player  of  6'10“.  The  6dB  difference  we  actually  observe  would  be  more  nearly 
consistent  with  a  proport.onality  to  the  square  of  the  scatterer  length,  implying 
a  4.56dB  increase. 

In  fact,  the  scatterer  long  dimension  is  about  a  third  of  the  wavelength  at 

O  a 

the  frequency  used  in  these  computations  and  tr parameter  is  about  2  in  free 

X 

space.  Hence  no  simple  proportionality  between  the  field  strength  and  the  scatterer 

length  can  be  predicted,  such  relationships  being  usually  applicable  only  to 

Rayleigh  scatterers,  very  small  compared  to  wavelength.  However,  for  the  cases 

computed,  we  do  see  a  positive  variation  with  intruder  height  L$  corresponding 

2  3 

roughly  to  somewhere  between  and  L$  . 
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11.4.5  Results  for  the  "Crawling  Man"  Cases 

Unfortunately,  not  all  of  the  crawling  man  cases  (Cases  IV-A,  B,  C,  D  and 
V-A,  B,  C,  D  of  Section  11.3)  were  run  successfully.  Attempts  were  made  to  run 
these  cases  immediately  before  the  contract  deadline.  Some  of  the  cases  were 
aborted  due  to  computer  difficulties  and  there  was  no  time  left  to  correct  these 
difficulties. 

Most  of  the  intended  results  were  obtained  and  these  are  presented  in 
Tables  11.4,  11.5,  and  11.6,  Columns  5  through  12.  A  subset  of  these  results  are 
shown  graphically  in  Figures  11.36  through  11.39. 

Comparisons  of  the  "radially  crawling  man"  (RCM),  the  "transverse  crawling 
man"  (TCM)  and  the  "upright  man"  (UM)  cases  are  best  made  from  Tables  11.4,  11.5 
and  11.6,  where  for  each  field  component,  the  amplitudes  for  $s  =  0°  are  shown 
in  Columns  1,  5  and  9  for  UM,  RCM  and  TCM  respectively.  At  $>s  =  90°,  the  cor¬ 
responding  columns  are  2  (UM),  6  (RCM)  and  10  (TCM).  At  =  180°,  we  have  3  (UM), 

7  (RCM)  and  11  (TCM)  and  at  $$  =  270°,  4  (UM),  8  (RCM)  and  12  (TCM).  We  now  focus 
attention  on  the  variation  of  the  scattered  field  components  at  the  antenna  with 
distance  of  the  scatterer  from  the  cable.  This  was  discussed  for  the  UM  cases  in 
Section  11.4.2.  Certain  RCM  and  TCM  cases  have  the  expected  behavior,  beginning 
with  low  values  at  p$=  17.5,  peaking  at  or  near  23.5  and  decaying  to  low  values  at 
28.  At  *s  =  0°  this  behavior  is  exhibited  by  |Ex|  and  |Ey|  for  RCM,  but  not  by 
1 E  [ ,  which  shows  a  more  nearly  flat  response  as  distance  changes  inside  the  cable, 
but  a  substantial  decay  as  the  scatterer  recedes  from  the  cable  radius  outside  the 
cable.  The  available  results  for  TCM  at  $s  =  0°  (not  all  of  these  are  available,  as 
indicated  in  Column  9  of  Table  11.4)  do  not  seem  to  indicate  that  kind  of  behavior, 
in  fact  they  depart  from  it  very  substantially.  At  =  90°,  there  is  a  small  decay 
in  all  the  components  of  both  RCM  and  TCM  between  17.5  and  23.5,  opposite  to  the 


11-27 


desired  trend,  but  again  a  substantial  decay  from  23.5  to  28.0.  The  rates  of  decay 
outside  the  cable  for  both  C's  =  0°  and  $s  =  90°  are  comparable  to  those  observed 
with  UM,  which  have  already  been  discussed  in  Subsection  11.4.2. 

The  RCM  and  TCM  amplitudes  for  <>  =  180°  and  <j>$  =  270°  are  in  general  more 
erratic  than  their  counterparts  for  <t> s  =  0°  and  $  =  90°.  Some  of  these  runs  were 
aborted  due  to  a  difficulty  in  one  of  the  subroutines  in  the  Barber  scattering 
program.  This  problem  had  not  been  encountered  before.  Since  it  occurred  very 
near  the  end  of  the  contract  period,  there  was  no  time  to  correct  the  difficulty. 
Some  very  small  changes  in  angle  assignments  were  made  to  overcome  the  difficulty 
in  doing  re-runs  of  these  cases.  Some  of  the  cases  then  ran  successfully  but 
produced  very  strange  results  which  we  believed  to  be  unreliable.  Consequently, 
some  of  the  points  are  not  presented  in  the  tables.  In  one  case  (TCM,  <t>s  =  270°) 
the  entire  run  was  aborted  and  we  have  no  results  on  this  case. 

Comparison  of  the  overall  order  of  magnitude  differences  between  UM,  RCM 

and  TCM  cases  at  a  given  value  of  p$  and  <t>s  can  be  made  from  the  tables.  At 

$s  =  180°  and  <t>s  =  270°,  too  few  computed  points  are  available  for  meaningful 

comparisons.  It  is  observed  that  for  a  =  0°  and  a  =  90°  (comparing  column  1, 

5  and  9  results)  the  values  of  jE  j  and  |E  |  for  RCM  cases,  through  the  range  of 

x  y 

values  of  p$  covered  tend  to  be  roughly  comparable  (on  the  average)  to  the  cor¬ 
responding  UM  cases,  but  the  corresponding  |E  |  values  in  the  TCM  case  (those 
available)  are  higher  by  between  10  and  20dB.  The  |E  |  values  exhibit  somewhat 
different  trends  at  $  =  0°.  The  TCM  values  are  nearly  comparable  to  the  UM  values 
but  the  RCM  values  tend  to  be  higher,  by  as  much  as  20dB  in  some  cases.  At 
*s  =  90°,  the  jE  |  values  for  TCM  are  as  much  as  70  or  80dB.  This  would  imply  that 
at  $  =  90°,  the  scattering  is  negligibly  small  for  the  RCM  case  relative  to  the 


The  difference  in  orientation  between  the  UM,  RCM  and  TCM  cases  would  not  alter 
expectations  for  inverse  distance  behavior,  which  are  based  on  the  naive  thought  that 
the  intruder  has  the  properties  of  a  point  scatterer,  i.e.,  subtends  a  small  angle  at 
both  the  source  and  the  receiving  antenna.  That  is  certainly  the  case  at  the  antenna, 
reasonably  well  into  the  scatterer's  far  zone,  but  not  so  at  the  nearest  illuminating 
slot,  which  sees  the  scatterer  as  an  extended  object  subtending  a  large  angle  and 
lying  in  its  near-zone,  particularly  at  radial  positions  near  24  meters. 

The  differences  in  orientation  do  produce  substantial  differences  in  the 
way  the  scatterer  is  illuminated  and  the  angular  position  of  the  antenna  in  the 
object's  scattering  pattern.  For  the  UM  cases,  the  illumination  and  scattering  are 
both  predominantly  broadside.  For  the  RCM  cases,  both  illumination  and  scattering 
are  primarily  "nose-on".  For  the  TCM  cases,  illumination  and  scattering  are  again 
largely  broadside,  but  with  a  much  smaller  illunimation  coverage  in  the  vertical 
direction  than  in  the  UM  cases.  All  that  can  be  said  is  that  these  differences  in 
geometry  should  result  in  some  differences  in  order-of-magnitude,  as  is  observed 
in  our  results.  But  these  observed  differences  are  not  always  in  the  same 
direction.  That  might  be  considered  anomalous  if  we  were  dealing  with  a  Rayleigh 
scatterer  illuminated  by  a  single  plane  wave.  In  our  model,  however,  we  have  a 

scatterer  (free  space  wavelength  -  5.26  meters)  whose  parameters  ^s  and  ^2Rs^ 

X  X 

in  free  space  are  respectively  about  2.125  and  0.596.  With  the  assumed  scatterer 
permittivity  (e$  =  46.2eQ),  neglecting  its  conductivity,  the  wavelength  in  the 
scatterer  medium  is  reduced  to  about  .774  meters.  Hence  the  parameters  Ls  and 

o  / op  X  X 

n4*s;  referred  to  the  scatterer  medium  are  now  about  14.4  and  4.06  respectively, 

x 

These  numbers  indicate  that  the  scatterer  is  large  enough  compared  with  wavelength 
to  produce  scattered  fields  that  may  be  highly  oscillatory  and  highly  directive. 


Changes  in  orientation  such  as  those  between  UM,  RCM  and  TCM  coes  may  well  produce 
very  pronounced  differences  in  both  the  orders  of  magnitude  of  the  scattered  fields 
and  their  sensitivity  to  scatterer  positions. 

■  Another  effect  which  complicates  attempts  at  interpretation  of  these  differences 
is  that  of  ground  reflection.  The  assumed  ground  parameters  (eg  =  4eo,  =  .002, 

°g  =  0.158)  are  sufficiently  large  to  produce  very  significant  ground  reflection 

“eg 

terms.  Coupled  with  the  highly  angle-dependent  scattered  fields,  the  ground 
reflections  of  both  illuminating  and  scattered  fields  should  further  increase  the 
sensitivity  of  the  results  to  small  changes  in  the  geometry. 
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11.4.6  Relative  Magnitudes  of  Components 

The  first  result  that  warrants  discussion  in  correction  with  this  topic  is 

that  near  the  beginning  of  Subsection  11.4  giving  the  amplitudes  of  the  x,  y  and  z 

components  of  the  electric  field  at  the  antenna  in  the  absence  of  a  scatterer.  In 

summary,  this  result  says  that  the  y-component  predominates,  the  x-component  is 

between  one  and  two  orders  of  magnitude  below  the  y-component  and  the  z-component 

is  comparable  in  magnitude  to  the  x-component.  To  attempt  to  explain  this,  we  first 

examine  Eq.  (I-3a)  in  Appendix  I.  We  also  observe  that  we  have  assumed  that  the 

TEM  mode  is  the  only  propagating  mode  in  the  cable.  For  the  TEM  mode,  the  electric 

field  is  entirely  radial  and  the  magnetic  field  is  entirely  azimuthal.  That  implies 

that  the  (n  x  E)  component  in  (I-3a)  vanishes,  the  (n  x  H)  component  is  entirely  in 

*  "  "  "  /IjkRN 

the  direction  of  the  cable  axis  and  the  (n  .  E)  term  is  in  the  direction  of  v'  Ur  j 

For  any  given  value  of  the  latter  is  in  the  direction  of  £. 

If  all  slot  contributions  were  in  phase  and  if  there  were  no  attenuation  along 
the  cable,  then  circular  symmetry  would  prevail.  Since  the  antenna  is  at  the  center 
of  the  configuration,  one  would  expect  a  cancellation  of  the  horizontal  electric 
field  contributions  from  the  cable  slots  and  a  summation  of  the  vertical  contributions 


based  on  the  observations  made  above.  The  assumed  attenuation  of  the  TEM  mode  is 
about  .000728dB  per  degree  of  angle  along  the  cable.  This  implies  the  following: 

(a)  if  the  field  strength  of  the  TEM  mode  were  unity  at  <j>  =  0°,  it  would  decay  to  a 
value  of  about  .927  at  90°,  .860  at  180°,  .798  at  270°  and  .739  at  360°.  These 
values  would  not  explain  entirely  the  x-y  asymmetry  that  is  observed  in  the  results. 

In  order  that  the  y-directed  field  at  the  antenna  will  greatly  exceed  the 
x-directed  field,  the  (n  x  H)  contribution,  which  is  the  predominant  one  for  the  TEM 
mode,  must  be  primarily  from  the  portion  of  the  cable  near  <f>  s  0°,  i.e.,  the  contri¬ 
butions  from  that  region  must  be  more  heavily  weighted  than  those  from  other  angular 
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regions.  Although  attenuation  could  play  some  role  in  that  predominance,  it  is 
unlikely,  as  indicated  above,  to  play  a  major  role.  The  relative  phases  of  the 
contributions  from  the  other  portions  of  the  cable  might  cancel  out  the  contributions 
from  these  regions  in  such  a  way  as  to  leave  the  predominant  contribution  from  the 
region  near  $  =  0°.  If  that  occurs,  then  that  would  explain  the  observed  predominance 
of  the  y-field. 

In  preliminary  results  presented  a  few  months  ago,  this  same  asymmetry  was 
observed,  except  that  the  2-component  was  even  smaller  relative  to  other  components 
than  in  the  present  results  and  the  y-component  was  larger  compared  with  the  x- 
component  than  in  the  present  results.  It  was  suggested  at  the  time  that  the  pre¬ 
sence  of  an  additional  mode,  subsequently  removed,  could  be  partially  responsible 
for  the  asymmetry.  If  that  was  true,  then  the  removal  of  that  mode  did  not  com¬ 
pletely  remove  the  asymmetry.  Moreover,  the  method  of  summation  of  the  contri¬ 
butions  around  the  cable  slots  was  changed  since  those  results  were  obtained  and  the 
asymmetry  still  persists.  Hence,  barring  an  error  in  the  computations,  it  must  be 
due  to  some  features  of  the  source  model  itself. 

In  the  numerical  results  for  the  illuminating  fields  with  ground  reflections 
(Subsection  11.4.1;  Tables  11.1,  11.2,  11.3  and  Figures  11.9-11. 16^ the  same  relative 
orders  of  magnitude  between  components  is  usually  observed  qualitatively  i.e.,  the 
x  and  z  components  are  of  comparable  magnitudes  and  the  y-component  is  about  one 
order  of  magnitude  (about  20dB)  above  the  x  or  z  component.  Without  ground  re¬ 
flections  (Figuresll.2,11.4,11.7),  the  z  component  tends  to  be  about  10  to  20dB 
below  the  x-component  and  the  y-component  is  again  larger  than  the  x-component,  but 
usually  not  as  much  larger.  This  indicates  that  ground  reflections  seem  to  enhance 
the  x-y  asymmetry. 
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Evidently,  the  same  mechanism  is  producing  these  asymmetries  in  the  results 
for  the  illuminating  fields  near  the  cable  as  is  producing  them  at  the  antenna 
location  far  from  the  cable.  Hence,  it  cannot  be  explained  by  the  differences 
between  the  fields  at  the  antenna,  which  are  entirely  far  zone  fields,  and  those 
at  the  potential  scatterer  locations,  which  are  near  zone  fields. 

The"Upright  Man"  cases  (Subsection  11.4.2;  Tables  11.4,  11.5  and  11.6  and 
Figures  11.17  -  11.23  )  and  most  of  the  "Crawling  Man"  cases  (Subsection  11.4.5; 
Tables  11.13,  11.14  and  11.15  and  Figures  11.36  -  11.39  )  do  not  exhibit  the  x-y 
asymmetry  observed  in  the  illuminating  fields  except  in  the  anomalous  case  <f>s  =  270° 
where  some  of  the  y-components  are  lower  than  their  corresponding  x-components  by 
as  much  as  20dB.  Otherwise,  x  and  y  components  tend  to  be  of  comparable  magnitude 
and  z-components  tend  to  be  below  the  x-components  by  between  10  and  20dB. 

It  would  seem  from  the  geometry  that  the  horizontal  field  contributions 
from  the  scatterer  at  $s  =  0°  and  180°  should  be  predominantly  y-directed  and  that 
those  from  <j>s  *  90°  and  270°  should  be  predominantly  x-directed  (since  the  scattered 
wave  appears  as  a  spherical  wave  emanating  from  the  scatterer  center).  Considering 
the  case  of  a  scatterer  at  <j>s  =  0°  or  180°,  the  x-contribution  comes  from  the  fact 
that  the  scatterer  center  and  antenna  are  at  different  elevations  and  hence  the 
wave  propagation  vector  for  the  scattered  wave  has  a  vertical  component.  This 
produces  an  x-component  of  the  free-space  scattered  wave  field. 

Without  ground  reflections,  this  x-component  might  not  be  as  large  as  the 
y-component.  There  would  be  only  a  y-component  and  a  z-component  if  the  antenna 
and  scatterer  center  were  at  the  same  height  and  only  a  small  x-component  if  the 
height  discrepancy  were  small.  The  ground  reflections  should  favor  vertical  over 
horizontal  polarization,  since  the  effective  ground  conductivity  is  reasonably  high 
at  57MHz.  The  antenna  and  scatterer  center  heights  are  large  enough  to  produce  a 


significant  vertical  component  of  the  wave  propagation  vector  for  the  ground- 
w  reflected  scattered  wave  which  is  vertically  polarized  on  ground  reflection.  The 
ground  reflection  of  the  y-component  of  the  scattered  wave  is  smaller  because  it  is 
horizontally  polarized.  The  vertically  polarized  ground  reflected  wave  contains  an 
x-component  which  is  added  to  the  x-component  already  present  in  the  free-space 
scattered  wave,  whereas  the  ground  reflected  contribution  to  the  y-component  is 
smaller  because  it  was  horizontally  polarized  on  ground  reflection.  All  of  these 
effects  are  in  the  direction  of  reducing  the  predominance  of  the  y-component  in  the 
wave  from  a  scatterer  at  <f>s  =  0°  and  180°.  The  same  arguments  can  be  applied  to  a 
scatterer  at  $s  =  90°  and  270°  except  that  x  and  y  are  interchanged. 


11.5  Conclusions  and  Suggestions  for  Future  Research 


The  range  of  parameter  values  for  which  results  have  been  obtained  is 
extremely  limited.  It  is  not  possible  to  draw  sweeping  conclusions  about  the 
generic  class  of  RF  intrusion-sensing  systems  exemplified  by  the  system  we  have 
attempted  to  analyze  on  this  project  (or  even  about  this  specific  configuration). 
Because  of  time  and  cost  limitations,  it  was  necessary  to  confine  ourselves  to  a 
very  specific  configuration  with  tightly  specified  parameter  values. 

The  electromagnetic  effects  we  have  investigated  here  are  extremely  com¬ 
plicated,  involving  near-zone  effects  and  a  scatterer  with  dimensions  comparable  to 
wavelength.  There  are  a  number  of  parameters  to  which  the  results  might  be  quite 
sensitive.  Among  these  are:  constitutive  parameters  of  the  ground,  which  are 
highly  sensitive  to  the  chemical  and  physical  composition  and  moisture  content 
of  soil;  the  constitutive  parameters  and  dimensions  of  the  “intruder"  (  the 
spheroidal  model  of  which  is  an  extremely  crude  approximation  to  a  human  body)  and 
the  geometric  parameters  of  the  cable,  small  changes  in  which  might  have  a  signif¬ 
icant  effect  on  the  fields  illuminating  the  scatterer. 

The  decision  was  made  early  in  the  project  to  construct  an  analytical  model 
based  on  the  concept  of  the  plane-wave  spectrum  of  an  electromagnetic  wave.  This 
was  based  partially  on  the  fact  that  the  scattering  program  we  planned  to  use, 
and  did  use,  presupposed  plane-wave  illumination.  The  decision  to  use  that  partic¬ 
ular  scattering  program,  in  turn  was  based  on  the  following  reasons.  First,  the 
Barber  scattering  program  appeared  to  be  the  best  one  available  for  the  purpose  at 
hand,  i.e.,  the  best  compromise  between  extremely  accurate  modelling  of  a  human 
frame,  requiring  absolutely  prohibitive  computer  time,  and  an  extremely  crude  but 
highly  inaccurate  model  requiring  very  little  computer  time.  Secondly,  the  program 
had  been  used  successfully  many  times  and  appeared  to  have  been  thoroughly  checked 


out.  Thirdly,  we  had  insufficient  time  and  resources  to  devise  a  different  program 
which  would  have  allowed  an  illuminating  field  other  than  a  plane  wave. 

Another  reason  for  using  the  plane-wave  spectrum  concept  was  the  desire  to 
include  ground  reflection  effects.  For  each  plane  wave  in  the  spectrum,  the  simple 
theory  of  reflection  of  a  plane  wave  by  an  infinite  boundary  can  be  used. 

For  the  above  reasons,  no  apology  is  made  for  using  this  form  of  analytical 
modelling.  However,  the  negative  aspect  of  this  choice  is  the  necessity  for 
performing  two-dimensional  Fourier  transformation  as  the  final  stage  of  the 
process.  At  the  initial  stage  of  the  analysis,  there  was  no  way  of  knowing  how 
computer-time  intensive  that  process  would  be.  After  it  was  too  late  to  turn  back, 
it  was  found  that  (1)  performance  of  the  two  dimensional  inverse  F.T.  on  the  scattered 
field  spectra  with  enough  points  for  high  accuracy  and  with  no  approximations 
reducing  the  number  of  points  needed  required  a  prohibitive  amount  of  computer  time, 
(2)  conventional  F.F.T.  algorithims  available  as  packaged  programs  would  not  be  of 
much  help  and  there  was  not  time  to  devise  an  F.F.T.  program  that  would  fit  our 
needs.  For  these  reasons  the  stationary  phase  method  was  used  as  a  crude  approxi¬ 
mation.  It  probably  reduces  accuracy  somewhat  but  also  reduces  the  required  computer 
time  very  substantially  and  brings  it  within  reasonable  ranqe.  However,  it  is  still 
too  high  to  allow  enough  computations  for  the  very  extensive  parametric  study  that 
we  believe  is  needed  to  develop  a  really  solid  analytical  understanding  of  the 
behavior  of  RF  intrusion  sensor  systems  of  this  class. 

A  set  of  preliminary  results  were  presented  several  months  ago,  together 
with  suggestions  for  ways  to  improve  those  results.  Nearly  all  of  the  suggestions 
were  carried  out  in  the  extension  of  the  contract  that  followed  (and  some  other 


improvements  not  suggested  were  also  made)*  except  for  the  suggestion  that  the 
use  of  stationary  phase  might  not  be  a  sufficiently  good  approximation  and  hence 
should  be  replaced  by  a  more  rigorous  method  of  performing  the  inverse  Fourier 
transformation.  An  attempt  was  made  to  follow  that  suggestion,  but  again  there 
was  insufficient  time  to  develop  a  new  and  possibly  more  accurate  approximation 
and  attempts  to  use  the  rigorous  F.T.  method  used  to  get  fields  in  the  absence 
of  the  scatterer  required  prohibitive  computer  time.  Hence  we  reverted  back  to  the 
use  of  the  stationary  phase  method  with  only  a  minor  modification  just  prior  to  the 
contract  deadline,  in  order  to  allow  sufficient  time  to  produce  the  planned 
numerical  results.  These  results  were  improved  relative  to  the  preliminary  results, 
but  we  are  still  not  satisfied  with  their  accuracy. 

A  recommendation  for  future  research,  which  we  plan  to  carry. out  in  continuing 
work  on  this  problem,  is  to  devise  an  approximate  method  of  inverse  Fourier  trans¬ 
formation  which  will  produce  hopefully  more  accurate  results  in  a  much  shot  ter 
running  time  per  point.  Once  that  is  accomplished,  we  will  re-run  all  of  the 
results  in  this  report  and  also  study  a  great  many  more  parameter  variations. 

Separation  of  the  various  constituents  of  the  field((l)  source  field  in  free 
space;  (2)  ground- reflected  source  field;  (3)  scattered  field  in  free  space;  (4) 
ground-reflected  scattered  field)  in  order  to  assess  the  effects  of  each  part  of 
the  process  will  be  one  aspect  of  a  future  study.  Another  aspect  will  be  an  improve¬ 
ment  in  the  accuracy  of  the  method  of  summing  over  the  cable  slots,  which  should  in 

turn  improve  the  accuracy  of  the  ratio  of  the  fields  at  the  antenna  with  and  without 
★ 

A  "field  mapping"  to  determine  the  spatial  variation  of  illuminating  fields;  removal 
of  a  second  propagation  mode  in  the  cable  which  was  used  in  the  earlier  computations 
reassessment  of  the  use  of  the  stationary  phase  method;  correction  of  an  erroneous 
value  of  scatterer  radius;  removal  of  all  phase  ambiguities  in  the  program;  improve¬ 
ment  of  the  accuracy  of  the  analytical  model  of  the  field  from  the  source. 
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the  scatterer.  This  problem  was  discussed  in  Section  11.4.6.  It  was  the  reason 
we  did  not  directly  address  the  issue  of  target  detectability  in  this  report. 
That  is  another  task  that  is  planned  for  the  continuing  work.  As  a  final  step, 
we  will  present  those  results  for  which  experimental  comparisons  are  possible  in 
a  form  amenable  to  direct  graphical  comparison  with  these  experimental  results. 
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Appendix  I:  KIRCHHOFF-HUYGHENS  (OR  STRATTON-CHU)  INTEGRAL  FORMULAS 

FOR  FIELDS  FROM  APERTURES 


A-1 

The  Ki rchhoff-Huyghens  Formula  for  the  fields  from  an  aperture  with  area 
AS'  in  an  unbounded  charge- free,  current-free  linear,  homogeneous,  isotropic 
medium  (neglecting  contributions  from  the  aperture  edges)  is  (see  Figure  1-1): 


V(r) 
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dS * [V(r ' )  n'  •  7*  ( 
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jkR  3V 
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Formulas 


where 


AS'  is  the  aperture  area 

dS'  ■  d£.|  d?2  *  aperture  area  element,  where  and  ^  are  orthogonal 
coordinates  along  the  aperture 

V(r)  ■  generic  vector  Indicating  electric  field  vector  E(r)  or 
•*  ■* 

magnetic  field  vector  H(r) 

0  *  origin  or  coordinates 
r  »  field-point  vector  *  xx  +  yy  +  zz 


r'  *  source-point  vector  (on  aperture)  s  xx'  +  yy'  +  zz' 
n'  3  outward  normal  unit  vector 


7' 


■  n*  •  7'  =  normal  derivative  operator 

R  3  r  -  r* 

R  3  |R | 


Note  that 


r'  3  r'  +  6r‘ 
-*■  -+o 


(1-2) 


where  r^,  position  vector  of  aperture  center  3 xx^  +  yy^  +  zz^ 

6r‘  =  vector  from  the  center  r'  to  a  point  on  the  overture 

-*-0 


'i*:'  n'l , 


n'*0 


where  (^,  c2)  are  the  unit  basis  vectors  in  the  and  c2  directions  respec- 
tively,  and  n‘  is  the  coordinate  in  the  n'  direction  (which  is  set  to  zero 

-b 

along  the  aperture). 

Note  that  Eq.  (1-1)  applies  to  either  the  electric  field  vector  E(r)  or 


the  magnetic  field  vector  H(r). 


Another  Integral  formula  similar  to  (and  equivalent  to)  Eq.  (1-1)  is 
sometimes  referred  to  as  the  Stratton-Chu  Integral  formula  [Stratton,  p.  466]. 
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Where  in  Eqs.  (I-3a,b)  we  neglect  contributions  from  aperture  edges 
[Stratton,  p.  469]  and  assuming  volume  current  and  charge  densities  to  be  zero. 

We  can  easily  derive  an  alternative  form  for  Eqs.  (I-3a,b)  in  which  the 
gradients  are  taken  with  respect  to  the  field  point  (unprimed)  coordinates 
rather  than  the  source  point  (primed)  coordinates.  We  recognize  that 
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(1-4) 


and 


E(r'),  H(r'),  n'  are  all  independent  of  the  field  point 

-►  -f* 

(unprimed)  coordinates  (1-5) 


1-3 


Using  Eqs.  (1-4)  and  (1-5)  In  Eq.  (I-3a,b).  we  have 
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Appendix  II:  PLANE-WAVE  SPECTRAL  REPRESENTATION  OF  FIELDS 


The  representation  of  electric  and  magnetic  fields  as  a  superposition  of 
plane  waves  is  amply  covered  in  the  literature*"1 ,A"2  It  is  a  well-established 
technique  for  dealing  with  arbitrary  electric  and  magnetic  fields  in  a  homo¬ 
geneous  medium.  In  principle,  it  Involves  no  more  than  the  three-dimensional 
spatial  Fourier  transformation  of  the  field  vector  or  one  of  its  components. 

As  in  Figure  II-l,  consider  a  right-handed  rectangular  coordinate  system 
(x,y,z)  with  unit  vectors  (x,y,z).  Assume  harmonic  time  dependence  e”ju)t. 

The  vector  r  *  xx  +  y£  +  zz  is  the  position  vector  referred  to  the  origin  of 
coordinates  0.  Electric  or  magnetic  field  vectors  are  denoted  generically 

by  V(r)  and  their  three-dimensional  spatial  Fourier  transforms  by  V(e).  The 
-*•  -*■  -*■ 

vector  k  is  the  propagation  vector  associated  with  a  plane  wave  propagating 
in  a  particular  direction.  The  angle  between  r  and  k  is  denoted  by  y. 


Geometry  for  Plane-Wave  Spectral 
Representative  of  Fields 


The  propagation  vector  is  denoted  by 


( J I  - 1 ) 


WlUV'  v  '  prop-iga  i' ion  cn!'r»tarit  of  Lht*  ir.'.i! i uni ,  where  i.i  ~  angular  f rcquer.ey , 
v  -  veior.i ty  of  .•l»ctro:.o<j:ietic  waves  in  the  medium,  which  in  the  cci:  e  of 

p 

f>  oo-'.p-roo  is  c  .5(10  ’)  n.e  ter/second  and  which  in  gr-noral  could  ho  real  (in  a 
n.'.n  dissipative  radium)  or  complex  (in  a  dissipative  ir.ydiu.a,  c.y.,  the  earth), 
md  3  is  the  direction  vector,  given  by 


B  -  Bu  +  zp. 


(II-2) 
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and  where  the  quantities  in  flq.  (!I-2)  are  given  in  terms  of  angles  by  the 
relationships 
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-  if  £<  V<  * 


( II-2-4) 1 


0^  -  spheric.il  polar  angle  of  propagation  direction 


( II-2-5) ' 


1 1-2 


*  azimuthal  angle  of  propagation  direction 


( I 1-2-6) ' 


and  the  plus  and  minus  signs  correspond  to  upward  and  downward  propagation, 
respectively. 

Integration  over  all  direction  angles  results  in  the  representation 


V(r)  •  V(£,z)  *  fj  d2  8„  e 
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( 1 1- 3) 


where  d28h  =  dBx  d0y,  o  =  xx  +  yy,  V_(|h)  and  V+(|h)  denote  the  two-dimensional 

spatial  Fourier  transforms  for  downward  and  upward  propagation  respectively 

and  where  the  z-component  6Z  has  been  eliminated  in  the  arguments  of  V_  and  V+ 

by  invoking  the  relationship  ( I 1-2-4) ' ,  by  virtue  of  which  0Z  can  be  expressed 

as  a  function  of  6  and  0  . 

a  y 

By  differentiating  Eq.  (II-3)  with  respect  to  z  and  setting  z  equal  to 
zero,  we  obtain 

jk0. -p  „ 

V(p.O)  ■  jj  dz  ^  e  ^  -  {V+(0 h)  +  V.(6h)}  ( I I -4a ) 
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It  is  easily  shown  from  Eq.  (II-4a)  (by  multiplying  by  e  *  *  and  inte¬ 
grating  on  x  and  y)  that 
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where  dS  *  dx  dy. 


An  alternative  representation  equivalent  to  Eq.  ( I 1-3)  in  terms  of  the 
spherical  angles  of  the  wave  vector  and  the  spherical  coordinates  of  r  (given 
by  r,  0,  <j>)  is: 


f2ir  rfr/2  „ 

V(r)  *  |  dq>0  {  J  deg  sin  6g  V_(0fil  q>6 )  ejkr  cos  y 


+  j  de6  sin  e6  V+(0e,  ♦g)  ejkr  cos  Y 


( II— 6) 


where 


cos  y  *  cos  e  cos  eg  +  sin  0  sin  eg  cos(q>  -  d>B ) 
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A;->andix  HI:  PLANE-HAVE  SPECTRAL  Rhf'HCSENrATTtV; 
Of-  IMF  Ffl  I  DS  FROM  AN  APERTURE 


The  spherical  wave  function  (r  —  appearing  in  Lqs.  (i  .l),  (1.3  j,:>) 
and  (I.6.a,b)  has  a  modifier!  plane-wave  spectral  representation  [fro.i:  a  mud 
i  F i cation  of  Stratton,  Chapter  IX,  Eqs.  (23),  (24),  (23) ,  (20),  op.  S77,  07;’.] 


do|,  sin  O'  c 
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ft  f(e,  R) 


O  O 

where  it  is  recalled  that  R  =  (Rj  =|r  -  r'l,  where  (o',..'}  are  the 

1  ■> 1  ■>  •>  1  ft  [S 

coordinate  direction  angles  of  t  referenced  to  a  z-coordinate  in  toe 
of  R,  and  where** 


(nu) 

spheri cal 
direr  Li  on 
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hoting  that  ejkR  cos  -  e'ikk,P-  ,  it  follows  that.  v(ejkv,R)  -  jkf.  ej!:f',R 
and  hence,  taking  the  gradient  of  (III.l)  (differentiating  under  the  integral 
sign)  we  have 
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The  relationship  appearing  in  Stratton  requires  the  use  of  complex  values  of  tno 
components  of  3,  the  propagation  vector,  for  its  implementation.  The  present 
form,  confining*  itself  to  real  values  of  the  3  components,  contains  a  dependence 
on  R*  Hence  it  is  a  modified  plane-wave  spedtral  representation  where  the  com¬ 
pound  ’waves  have  phase  factors  eJ'kf!-R  as  do  plane  waves,  but  are  weighted  with  a 

factor  that  is  dependent  on  the  space  coordinates.  . " 

-The  derivation  of  (III.l)  is  easily  accomplished  by  recognizing  that  * do!, 
sin  o'  ejkRcos  Qp  =  j  sinkR.  Differentiation  of  sinkR  then  produces  the°result 

H  .  fir  -•  i.  r, - ......  n 

and  j  sin  kR  , 
implying  that 


sin  o'  e^kRcos  03  =  j  sinkR.  Differentiation  of  sinkR  then  produces 

cos  k.R  =  X  ['T  do'  sin  Ojj.  ejkKc0S  °K{l+jkR  cos  op.  Adding  cos  kR  j 
dividing  ^by°R  ahd  noting  that  the  integrand  is  independent  of  ,  m 

f?J'  d4> *  =  2ir,  we  obtain  III.l. 
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Revising  (Iil.l)  and  (111.2)  to  include  the  separation  into  upward 
and  downward  vmv<>  .  and  transforming  from  the  spherical  coordinate  system  with 


2 -ax is  in  the  I*  direction  to  our  ground  frame,  v/e  have 
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where  b  =  p.  -t  z  ]  r<  1 

We  now  simplify  the  notation  hy  denoting  f  (j^ ,  R)  by  f"  ,  neglect  the 
R 

term  in  (III. 2)’  ,  neglect  the  field  contributions  due  to  line  currents 
around  the  aperture  edge,  assume  that  the  ambient  medium  is  free  space,  and 
then  substitute  (III.l)1  and  (III. 2)'  into  Eqs.  (I.3-a,b),  interchanging  the 
order  of  integration.  The  results  of  these  operations  are: 
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*  Tliis  term  is  not  included  in  the  analysis  in  Sections  3  and  4,  but  it 
is  approximately  accounted  for  in  the  computations.  It  is  negligibly 
small  except  at  positions  extremely  close  to  the  source. 
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represent  upward  and  dav/nwcmi  propagation,  respectively. 
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where  Y„  =  /—  -  J-  -  wave  admittance  of  free- space  and  +  and  -  super serials 
0  vft>0  Zo 

have  the  same  meaning  as  in  Eq.  (I II. 3a). 

We  now  express  E(r)  and  H(r)  in  the  fon.i  (II. 3)  as  follows: 


E(r)  --  E+(r)  +  E'(r) 


(IIT.flti) 


H(r)  -  H  (r)  +  if  (r) 
->  ->  -*  •>  ->• 


(lll.4b) 


wnere 


e-(j)  .  b<£„> 


.  t'r  o  Jk  (B.  -p+(B  |x)  .. 

Hd  *  jj d  a, e 


and  where,  from  Eqs.  (III.3a,b), 

"(j  k  )?  rr  -jk  (Bh+z|B  1)t’ 

W  ■  —fir--  «’ «  *  <-V.V  *  8<fiJ 

AS' 


i 


-  |{i'm  ■'  y\ x  (iV  x  r  ( >■  • ) )  -  (n'  •  l  (.•'))(  \  »  v: .:])'] 

>n  .  /  >  »  *  f  *  *  ■  ’  *  •  * 


(in.-.,) 


-.1 KM-, 


o  ->r, 

to  e 


*  *  (Yn[»i*  X 

O  r  -  > 


[(Bh  +  2]!i2i)  x  (n*  x  M  ( r 1 } )  -  (o'  -  l!(  r ' ) )  { f  |!^|)]}  t' 


An  important  point  arises  a  t  this  juncture.  It  is  necessary  that  E+(3h) 

and  fi  (f,.)  be  interrelated  through  the  Maywoll  equations,  i.e.» 

•  /■  n 


lA>  -  -Vr  * 


(III.Go) 


•uw  =  Vr *  y^'1 


(in.et.) 


,-+  -E(gh)  ^  0 
-  >-  >  ->  n 


(III.6c) 


£r-  •  H  ((••..)  =■  0 

V  -►  >11 


(111.5c 


Since  Eqs.  (III. 6a)  and  (III. 6b)  imply  Eqs.  ( I ! I -6c)  and  (Iil.ftd), 
respectively,  and  Eq.  (III. 6a)  [with  the  aid  of  (UI.Gd)J  implies  Eq.  (III. 6b), 
it  follows  that  if  Fqs.  (III. 6a, b)  are  consistent  with  either  Eq.  (III. 6a)  or 
(III. fib),  then  Eqs.  (IH.F»a,b)  are  consistent  wi th  the  entire  system  of  Eqs. 

( 1 1 1 .  6a ,  b ,  c. ,  d ) . 

To  test  this  consistency,  we  apply  the  (-Z  pr-  x  ...)  operation  to  Eqs. 

( III .!>)-(III .6),  resulting,  after  some  manipulation  (using  the  vector  identity 


nr-4 


”■  T  i  ■  fff- 


a  x  (l>  x  c)  “  f)(a  •  c)  -  c(<t  •  h)  and  the  fact  that  r;f-  x  -  0,  ft1-  •  (;'t- 

>  '  >  ■>  >  r  ■>  y 

in  the  expression 

k2 


-Z  ft—  X  ft  f  X  ) 
o>  x  .,V  ,ir 


-(>  KX 


w 

AS 1 


-jk  !v  -r‘ 

0>  >  [  l /  [  1 1 1  v  !!/  V-*  i 

0  , 


<«S*  o  u>  >  !«;,r  xtnr'j 


-  [ft-  x  (n1  x  E(r 1 ) )  -  p+ [  n  *  •  E(r’))] 


y  -> 


r  >  > 


+  n+  n*  •  [-Zfft1  X  il(r’))  -  C  ( r ' ) )  ’]  }  f  (III. 7) 

*>  •>■  U  >  >  >  -> 


Expressing  E(r‘)  and  H(r‘)  in  the  last,  square-bracketed  quantity  in  Eq. 
>  >■  ■>  > 

(III. 7)  in  terms  of  its  spectral  representation,  we  have 


-Z0[t5-  x  H ( r * ) ]  -  E(r’) 

U  >  -V  -x  > 


o  jk  ft’— r’  .. 

0>  >■  I..7  y  O  (R«*^ 

1  oV1;  v»-l!?.h" 


d  a e 


(HI  .8) 


When  we  evaluate  the  integral 


I  = 


dS1  e 


jk  (;V  -rf+) *r' 


(III. 9) 


which  arises  in  Eq.  (III. 7),  if  AS1  were  an  infinitely  large  surface,  then 


I  *  6(3-  -  &-) 

•>  -y 


(II  1. 10) 


which  implies  that 


3-’  »  (X 

-y  *> 


(in. n) 


l  n-  :» 


?rssr*s**^ 


If  Lq.  (111. 11)  iio  Id-, ,  then  the  RMS  of  Lq.  (III.  8)  mu'  v.tui.u,  sides;  t!.  ■ 
vanishing  of  ».*i< *  in Log rand  on  the  fills  of  Eg.  (Ill  .f!)  is  ii.»: •'!  led  ny  the  plan"- 
w..ve  Maxwel  I  equation  (111.6;).  Thus,  the  L.l IS  of  lq.  (’II. o)  must  a! so  vanish, 
which  ieu  that  I  lie  III  IS  of  Eq.  (III. 7)  is  equivalent  to  i._|(f'-^)  as  given  by 
Lq.  (111.5a).  This  completes  the  demonstration  that  Eos.  (fit. 5a, b)  are 
consistent  with  Eqs.  ( 1 1 1 .6a,b,c,d) . 

01  course,  the  "infinitely  large  surface"  assumption  cn  which  Eq.  (ill. 10) 
is  based  is  not.  really  valid  here.  However,  for  Egs .  (I.3a,b)  (cn  which  this 
entire  analysis  is  based)  to  be  truly  valid,  it  is  necessary  that  the  surface 
be  closed.  In  the  present:  model,  the  fields  are  assumed  to  be  negligibly 
small  except  on  the  slots;  however,  the  surface  still  includes  the  cable  area 
outside  the  slots  and  the  integral  1  of  Eg.  (III. 9)  still  can  be  considered 
as  an  integral  over  a  closed  surface.  Such  an  integral  will  vanish  unless 
f.1*  is  equal  to  p— ;  hence  Eq.  (III.  11)  holds. 

Also,  as  remarked  earlier,  we  have  neglected  ,  in  the  above  analysis,  the 
R  e.ikR 

>  term  arising  in  the  gradient  of  „  .  This  term  is  only  significant  at 

R7  R 

positions  extremely  close  to  the  source.  It  will  be  included  in  the  computa¬ 
tions. 


Appendix  IV:  ELECTROMAGNETIC  WAVE  PROPAGATION  ALONG  A  COAXIAL  CABLE 


The  theory  of  electromagnetic  wave  propagation  along  a  coaxial  cable  Is 
well-known;  consequently  we  will  Invoke  standard  references^  fn  developing  the 
theory  in  the  context  of  the  specific  problem  that  Is  the  subject  of  this 
report. 

First,  consider  the  modes  of  propagation  down  a  guiding  structure  with 
circular  cross-section.  Using  cylindrical  coordinates  (r",  z")  and  cor¬ 

responding  unit  base  vectors  (r",  5".  z"),  we  have  in  general  a  set  of  TE 
(transverse  electric)  and  TM  (transverse  magnetic)  modes  of  propagation.  The 
fields  between  r"  »  a  and  r"  a  b,  where  a  and  b  are  inner  and  outer  radii 
respectively,  can  be  expressed  as  a  superposition  of  TE  and  TM  modes.  The 
field  vectors  have  the  forms  (where  harmonic  time  dependence  e_Jwt  has  been 
assumed) 


”|(r“ ,z" ) 
H(r"  ,<t>"  ,z" ) 


l 

n»-° 


(C^E)(r»)  e' 
nH 


jno” 


'jko^zn)TEz" 


(IV. 1) 


where  in  general 


cirvi 


jko(k 


1  z" 
zn'TVT  , 


JnE 

nH 


kzr/k 


0 


(kzn^TE  3  ^znR^TE  +  ^Vte 
TM  TM  TM 


and  where 
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( v pn 5 t|t  3  phase  velocity  In  meters/sec  for  TEn  mode 
(vpn)iM  *  phase  velocity  in  meters/sec  for  TMn  mode 
(an)TE  3  attenuation  in  nepers/meter  for  TEn  mode 
(an)im  3  attenuation  in  nepers/meter  for  TMn  mode 

TE  From  Reference  A- 3,  Chapter  3,  we  can  write  the  general  forms  of  the  vectors 

SnE  as  ^ol1ows: 

^nH 

TE  modes  (T^(r")  =  0] 

cj|(r")  3  r"  tj.2)(r")  +  f  ^(r1')  (IV. 2a) 

C JJ(r")  3  r"  H^r”)  +  H^r")  +  z"  fljj^r")  (IV. 2b) 


where 


H^V)  *  /En  Jn(kcn  r")  ♦  BTE°  Nn(kcn  r") 


_ J  ^2n  TE  TE 

*V"V)  *  [A  n  J'{krn  r")  ♦  B  n  N ' ( k  r")] 


n'  cn 


n'  cn 


cn 


_ j  n  k  TE  TE 

V<r") - *-^[A  nOn(kcnr")  +  B  n  Nn(kcn  r">] 

p  kcn 


E^Cr")  =  ~jf  fll?)(r") 


w 


eIh ^(r") 


-  — H^M)(r") 


and  where 


p  =  kcn  r" 


Jn^kcn  r")  *  Vessel  Unction  of  first  kind  and  order  n 


Nn(kcn  r")  *  Bessel  function  of  second  kind  and  nth  order 


TE. 


zn 


Y,.  *  wave  admittance  for  TE„  mode  * 

w  n  a)  y 


kzn  and  k2n  are  given  below  Eq.  (IV. 1) 


IV-3 


A  £  A  j  0 

e  *  -c--c.  s  e  +  ca  _  complex  dielectric  constant 
cac  e0  ca  “  eQ  of  cable  material 

-  £ca 

£ca  =  *  where  £ca  =  Permittivity  of  cable  material 

aca  =  conductivity  of  cable  material 
TEn 

A  *  amplitude  of  Jn(p)  term  of  TEp  mode 
TEn 

B  =  amplitude  of  Nn(p)  term  of  TEn  mode 


(IV. 3a) 


MMMMM 


cJ5J(r")  -  r"  fl^r")  +  J"  H^(r») 


(IV. 3b) 


where 


,ns  TM  TM„ 

F7"V')  *  A  Jn(kcn  r")  +  B  n  N„(k_  r“) 


n'  cn 


/_\  j  k,_  TM  in 

rj^r")  =  -^[A  J'(kn  r")  +  B  n  N'(k  n  r")] 


n'  cn 


/n>  j  n  k  TM  TM 

E<I))(r“)  *  — ^  [A  Jn{kcn  r")  ♦  B  "  Nn(kcn  r-)] 


H^r")  .  Y™n  r(;i(r") 


where  all  quantities  used  In  Eqs.  (IV.3a,b)  were  defined  below  Eqs.  (IV.2a,b) 
with  the  exceptions 


A  n  =  amplitude  of  Jn(p)  term  of  TMn  mode 


n 

B  *  amplitude  of  Np(p)  term  of  TMn  mode 


TM  oj  £ 

Yw  =  wave  admittance  of  TMn  mode  *  ^ -ca 


1: 


The  TEM  (standard  coaxial)  mode  pT^(r")  *  H^(r")  *  0] 

At  low  frequencies  we  can  limit  propagation  in  a  coaxial  line  to  the 
transverse  electromagnetic  (TEM)  mode,  which  has  the  form 


(n) 


E(r\  z" )  =  r“  Er„ 


V  jk  z" 

W  3  r  msfepr  e 


(IV. 4a) 


H(rH,  $\  z”)  *  $'*  Hirj(r“)  *  <P" 


»  u(n )/„„%  _  'ca  ’0  ”ca 


Y  V  ik  z" 

ca  vo  J  *'»z 

in(b/a)r" 


(IV. 4b) 


where  k  ,  is  defined  below  Eqs.  (IV.2a,b)  and  where 

Cd 


a  *  inner  radius  of  cable 


b  *  outer  radius  of  cable 


VQ  =  voltage  between  r"  *  a  and  r"  *  b 


ca 


-  wave  admittance  of  cable  material 


Equations  (IV.4a,b)  constitute  a  complete  solution  for  the  TEM  fields. 
Equations  (IV.2a,b)  and  (IV.3a,b)  require  the  imposition  of  boundary  con¬ 
ditions  at  r"  =  a  and  r"  =  b  for  completion.  In  the  actual  cable  that  is 
used  in  the  configuration  under  study  on  this  project,  the  inner  boundary  at 
r"  =  a  may  be  well  approximated  by  a  perfect  conductor.  The  outer  boundary 
at  r"  *  b,  since  it  contains  the  slots,  may  be  approximated  as  a  perfect 
conductor  except  in  the  slot  region.  An  exact  analysis  would  require  con¬ 
sideration  of  a  boundary,  that  is  not  cylindrlcally  symmetric,  i.e.,  partially 


i 
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free  space.  Such  an  analysis  would  not  be  feasible  within  the  time  limita¬ 
tions  of  the  project;  hence,  the  fields  are  being  approximated  as  if  the 
cable  had  perfectly  conducting  Inner  and  outer  boundaries. 

The  result  of  these  approximations  is  the  condition  that  the  tangential 
electric  fields  must  vanish  at  both  r"  =  a  and  r"  *  b;  hence, 

E^}(a)  =  I^}(b)  =  E^}(a)  =  ^(b)  =  0  (IV. 5) 


for  both  TE  and  TM  modes. 

Applying  the  conditions  (IV. 5)  to  Eqs.  (IV.2a,b)  and  (IV.3a,b),  we  obtain 
the  following  results: 


TE  Modes 


J;<kcn  a> 


TE. 


+  D 


Wkcn 


a) 


(IV. 6a) 


J'(k 
nv  cn 


J  +  D 


TE. 


Wkcn  b> 


where 


TE. 


(IV. 6b) 


TM  Modes 


TM, 


where 


From  Eqs.  (IV. 6a)  and  (IV. 6b),  we  obtain  the  conditions  required  to 
determine  kcn  for  TE  and  TM  modes  respectively.  Dividing  Eq.  (IV. 6a)  by 
Eq.  (IV. 6b)  and  Eq.  (IV. 7a)  by  Eq.  (IV. 7b),  we  obtain: 

For  TE  modes 

J;<kcn  a>  N;<kc„  b>  -  J;<kcn  b>  N;<kcn  a>  *  0  <IV-8a> 

and  for  TM  modes 

J(k„„  a)  N  (k  b)  -  J  (k  b)  N  (k  a)  *  0  (IV. 8b) 

n  cn  n  cn  n  cn  n  cn 

In  the  configuration  under  investigation,  a  and  b  are  extremely  small 
compared  with  wavelength.  Under  the  very  plausible  assumption  that  kcn  is 
of  at  least  the  same  order  of  magnitude  as  kQ  or  at  the  very  least  no  greater 
than  one  order  of  magnitude  greater  than  kQ,  a  study  of  parameter  values 
indicates  that  the  argument  kcn  a  and  kcn  b  appearing  in  Eqs.  (IV.8a,b)  are 
all  very  small  compared  with  unity.  Hence  the  first  few  terms  of  the  power 
series  for  the  Bessel  functions  with  small  arguments  can  be  used  to  put  Eqs. 
(IV.8a,b)  in  a  form  where  they  can  be  easily  solved  for  kcn. 

Also,  it  follows  from  either  Eq.  (IV. 6a)  or  (IV. 6b)  that 
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(IV. 9a) 


J;(krn  a)  J'(k  b) 
n  cn  _  n  cn 

N'(Ca)  "  mrrsj 

n  cn  '  n  cn  ' 


and  from  either  Eqs.  (IV. 7a)  or  (IV. 7b)  that 


JAn  *> 


JAn  b> 


(IV. 9b) 


The  fields  on  the  slots  are  those  corresponding  to  r"  =  b.  Since  those 
are  the  fields  that  must  be  evaluated  in  our  problem,  we  will  focus  on  them 
at  this  point. 

From  Eqs.  (IV.2a,b),  (IV.3a,b)  and  (IV.9a,b),  we  have,  for  the  fields  at 
r"  =  b: 


TE  Modes 


TE. 


=  Nn(kcn  b)  in^cn  ^ 


H^b)  =  E^}(b)  =  E^(b)  -  0 


te  ,  _  r 

u(n)  3  a  J  kzn  ,v  .  N 

V  N;(kcnb)^T~Ln(kcn  b) 


cn 


(IV.lOa) 


(IV. 10b) 


(IV. 10c) 


v  n  Y 


(IV.lOd) 
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where 


Ln(p)  ■  Jn(p)  N^(p)  -  Nn(p)  J^(p) 


TH  Modes 


-l<5>(b)  «  H<:>(b)  -  =  0 


(IV. 11a) 


,n)  *  n  j  k 

Er"  (b)  *  1 - Vkcn  b) 

cn 


(IV. lib) 


T^S*Cb>  -  vw  "  t^'lb) 


(IV. 11c) 


A_4 

The  Bessel  function  series  have  been  evaluated  to  fourth  order  in  the 
argument  p,  based  on  the  observation  that  (k  b)  «  1  in  our  problem.  The 

v*  1 1 

results  are: 


J0(p)  *  1  “  (  f  )2  +  *  (  f  )4 
J-,(p)  -  f  (1  "  \  (  f  )2  +  T7  (  £  )4> 


J2(p)  »  \  0  '  f?  +  > 


(IV. 12a) 


(IV. 12b) 


(IV. 12c) 


For  n  >  3 


Jn(p)  =  -f—  (1 
n  2n  n! 


2n  +  2)(2n  + 


.}  (IV.12d) 
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j  ) 

i 

I  , 


J-j(p)  s  \  (I  -  -j  (  f  )2  +  ■fjr  (  j  )4) 

J2(p)  *  f  {1  _  T  +  } 


For  n  >  3 


Tfn\  ^  pn"^  r n  (n  +  2)p4  .  (n  t  jjp1*  , 

Jn(p)  '  Zn~  {n  “  2(2n  +  2)  2.4(2n  +  2)(2n  +  4)  } 


2"  n! 


n0(p)  91  f  Unf  +  yJ  J0(p)  +  |  {(  f  )2  ‘  |  (  f  )4> 
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» 2  Un  (  £  )[i  -  (  £  )2  (  £  )4] 

+  Cy  +  (  £  )2(1  -  y)  +  £  (  £  )4(2y  -  3)]} 

where  y  3  0.5772156 

N,'»)  *  ■  h 

hM  .  -  A  ( I  m  |  )2) 


(IV. 13a) 

(IV. 13b) 

(IV. 13c) 

(IV. 13d) 


(IV. 14a) 

(IV. 14b) 

(IV. 14c) 
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For  n  >  3 


Nn(P)  *  -  l  (  l  )n((n  -  1):  +  (n  -  2) !(  §  )2  +  (  £  )4} 

(IV.14d) 
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For  n  >  3 
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From  Eqs.  (IV. 12a)  through  (IV.lSd), 


t0(p)  *  J0(»)  «;<p>  -  j;(b)  n0(o)  -  £ 
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(IV. 15a) 


(IV. 15b) 


(IV. 15c) 


(IV. 15d) 


(IV. 16a) 


£-j(p)  B  ^(p)  N|(p)  -  Oj(o)  N^p) 

*  w  0  '  (  ?  }Z  +  7  (  f  )4}  (IV. 16b) 

=  J2^p^  N2^  "  J2^  N2^p^ 

■^d-i(f)4)  ( IV.  16c) 

For  n  >  3 


Ln(P)  -  Jn(p)  n;(p)  -  j;(P)  Nn(P) 


(n2  -  1) 
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(  §  )4)  ( IV.16d) 


The  forms  (IV. 12a)  through  ( I V . 1 6d)  are  those  actually  used  in  the  com¬ 
puter  program.  In  view  of  the  small  values  of  p  applicable  to  our  problem, 
these  forms  are  adequate  approximations  for  these  functions. 

An  aspect  of  this  analysis  we  have  not  yet  discussed  is  the  determination 

of  the  attenuation  an  for  each  mode.  The  usual  analytical  procedure  for  the 

A- 3 

case  of  nearly  perfectly  conducting  walls  is  as  follows:  (a)  Ignore  the 

departure  from  the  perfectly  conducting  wall  case  in  calculating  the  fields 
inside  the  (possibly  lossy)  dielectric;  (b)  Consider  the  actual  conductivity 
of  the  wall  material  and  calculate  the  "surface  Impedance"  along  the  walls; 

(c)  Calculate  the  "surface  current"  J.  5  n  *  H,  where  H  is  the  loss-free 

4  4- 

magnetic  field  at  the  surface  and  n  the  normal  unit  vector;  (d)  Using  the 
surface  impedance  and  surface  current,  calculate  the  power  loss  per  unit  leng*w 
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of  line  due  to  ohmic  losses  along  the  wall  surface;  (e)  Using  the  Poynting 
theorem,  calculate  the  losses  due  to  the  nonzero  conductivity  of  the  lossy 
dielectric  material,  assuming  of  course  that  the  conductivity  i_s_  nonzero; 

(f)  Add  the  power  loss  per  unit  length  due  to  both  mechanisms  (d)  and  (e)  and 
calculate  the  resulting  total  attenuation  an< 

It  was  originally  planned  to  carry  out  the  procedures  (a)  through  (f) 
above  for  each  mode  and  thereby  evaluate  the  attenuation  an  by  purely 
analytical  means.  It  was  decided  later  that  the  degree  of  approximation  of 
those  calculations  was  even  greater  than  that  associated  with  the  evaluation 
of  the  fields  in  the  cable  as  if  the  walls  were  perfectly  conducting.  Also, 
in  the  actual  cable,  the  slots  should  play  a  prominent  role  in  the  introduction 
of  another  source  of  energy  loss.  Obviously  these  slots  will  also  perturb  the 
magnitude  and  direction  of  the  field  vectors,  but  that  effect  is  probably  not 
as  severe  as  the  degree  to  which  they  will  affect  the  attenuation.  All  of 
this  is  speculation  of  course,  but  based  on  such  speculation  and  the  computer 
time  and  other  limitations  on  the  project,  it  was  decided  to  use  empirically 
determine J  values  of  attenuation.  That  is  what  was  done  in  the  computer  pro¬ 
gram. 

Another  issue  that  has  not  yet  been  discussed  in  detail  is  that  of  the 
determination  of  kcn  for  each  mode.  This  can  be  determined  theoretically  by 
solving  the  equations  (IV.8a,b)  for  kcn.  We  invoke  the  assumptions: 


a  |  <<  1 

(IV. 17a) 

b|  «  1 

(IV. 17b) 

thereby  justifying  the  use  of  the  approximations  leading  to  Eqs.  (IV. 12a) 


through  (IV.15d)  in  constructing  the  expressions  on  the  left  hand  side  of 

Eq.  (IV.8a,b).  We  then  construct  these  expressions  and  truncate  the  resulting 

series  beyond  the  fourth  power  in  kcn.  The  results  of  these  operations  are 

2 

two  quadratic  equations  in  kcn,  as  follows: 

TM  Modes 


A  k  +  B 
n  cn  n 


kcn  + 


V 


(IV. 18a) 


TE  Modes 


D  k4  +  E  k2  +  Fn  =  0 
n  cn  n  cn  n 


(IV. 18b) 


where 


A0  «  in  t2  *"<  £  )[a4  +  b4  +  4a2  b2]  +  3[a4  -  b4]} 
B0  -i  ((b2  -  a2)  -  (a2  +  b2)  £n(  |  )} 

c0  •  «(  |  ) 

Dq  =  a2  b2[(a2  -  b2)  +  (a2  +  b2)  £n(  |  )] 

Eo  =  -[2(a4  -  b4)  +  a2  b2  *n(  |  )] 


F 


o 


=  16(a2 


b2) 


IV-15 
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A  *  ■  ■  -  —  -  -  -  aa.sa 

n  32n(n2  -  4) (n2  -  l)(ab)n 

—  { (b2n  -  a2n)  n2(a2  -  b2)2  +  3n(a4  -  b4)(b2n  +  a2n 
+  2(a4  +  b4  +  4a2  b2)(b2n  -  a2n)} 


B  * - — J - {(a2n  +  b2n)  n(a2  -  b2)  +  (b2n  -  a2r 

n  4n(n2  -  1 )(ab)n 


C_  =  — [b2n  -  a2n] 


n  n(ab)n 


3"  32n(n2  -  4) (n2  -  1  )(ab)n 


—  { ( b2n  -  a2n)[-n4(a2  -  b2)2  +  10n2(a4  +  b4  -  a2  b2) 
+  8a2  b2]  +  (b2n  +  a2n)[n3(a4  -  b4)  +  8n(a4  -  b4)]} 


‘n  n 


#h2n  .2nw  2  ,  w2\ 
nB  +  AS - ",  a  ■  b  ? 

4(ab)n 


Fn* 


•n2  C. 


)(a 


Equations  (IV. 18a, b)  can  be  solved  for  k 


The  solutions  are: 


(IV. 19a) 


^cn^TM  mode 
n 


^kcn^TMn  mode  *  + ^ 


-E  +  k 2  -  4D  F 
n  -  V  n  n  n 


to: 


(IV. 19b) 


A  small  program  was  written  to  implement  the  calculations  embodied  in  Eqs. 
(IV. 18a, b)  and  (IV. 19a, b).  From  the  k  values  and  the  equation 

V'l 


k 


zn 


(IV. 20) 


we  can  determine  the  values  of  kzn,  the  complex  propagation  constant  for  a 
particular  mode. 

Since  the  real  part  of  k2n  is  uj/vp,  where  vp  is  the  phase  velocity,  and 
the  imaginary  part  is  a  portion  of  the  attenuation  an  (it  does  not  include 
attenuation  due  to  wall  losses),  it  would  seem  that  the  calculations  repre¬ 
sented  by  Eqs.  (IV. 18a, b),  (IV. 19a)  and  (IV. 20)  would  be  sufficient  to  deter¬ 
mine  the  propagation  and  attenuation  (exclusive  of  wall  losses,  discussed 
elsewhere)  properties  of  a  TEn  or  TMn  mode. 

In  the  numerical  computations  presented  in  this  report,  we  have  not  used 
the  above  formalism  to  obtain  the  propagation  constants  kzn.  Instead,  we 
have  used  empirically  determined  values  of  the  phase  velocity  and  the  attenua¬ 
tion,  which  circumvents  the  need  for  these  computations. 

However,  the  computation  of  both  phase  velocity  and  attenuation  for  each 
mode  directly  from  the  theory  is  quite  feasible  and  simple  to  carry  out.  If 
there  were  justification  for  carrying  out  these  computations  in  a  practical 


IV-18 


problem  (l.e..  If  It  were  known  that  a  number  of  higher  order  modes  were 
propagating  and  If  their  relative  amplitudes  could  be  determined).  It  would 
be  very  easy  to  add  this  capability  to  our  general  computer  program.  At  the 
relatively  low  frequencies  of  Interest  In  this  particular  problem,  most  of  the 
higher  order  modes  would  be  highly  evanescent  and  would  contribute  energy  only 
from  those  slots  near  the  power  source;  hence,  they  would  not  be  very  important 
contributors  to  the  fields  Incident  on  the  scatterer  or  the  fields  at  the 


antenna. 
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